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Abstract 



The dynamical triangulations approach to quantum gravity is investigated in detail for 
the first time in five dimensions. In this case, the most general action that is linear in 
components of the /-vector has three terms. It was suspected that the corresponding 
space of couplings would yield a rich phase structure. This work is primarily motivated 
by the hope that this new viewpoint will lead to a deeper understanding of dynamical 
triangulations in general. Ultimately, this research programme may give a better insight 
into the potential application of dynamical triangulations to quantum gravity. This thesis 
serves as an exploratory study of this uncharted territory. 

The five dimensional (k, I) moves used in the Monte Carlo algorithm are proven to be 
ergodic in the space of combinatorially equivalent simplicial 5-manifolds. A statement is 
reached regarding the possible existence of an exponential upper bound on the number 
of combinatorially equivalent triangulations of the 5-sphere. Monte Carlo simulations 
reveal non-trivial phase structure which is analysed in some detail. Further investiga- 
tions deal with the geometric and fractal nature of triangulations. This is followed by a 
characterisation of the weak coupling limit in terms of stacked spheres. Simple graph the- 
ory arguments are used to reproduce and generalise a well-known result in combinatorial 
topology. Finally, a comprehensive study of singular structures in dynamical triangula- 
tions is presented. It includes a new understanding of their existence, which appears to 
be consistent with the non-existence of singular vertices in three dimensions. The thesis is 
concluded with an overview of results, general discussion and suggestions for future work. 
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Crynodeb 



Mae dull triongliannau dynamegol i ddisgyrchiant cwantwm yn cael ei ymchwilio'n fanwl 
mewn pum dimensiwm am y tro cyntaf. Yn yr achos yma, mae gan yr effaith mwyaf 
cyffredinol, sydd yn llinol yng nghydrannau y /-fector, dri term. Roedd yn debygol y 
byddai'r bwlch cyplysau cyfatebol yn ildio strwythur gwedd cyfoethog. Prif gymhelliant 
y gwaith yw'r gobaith fod y safbwynt newydd yma'n arwain at ddealltwriaeth dyfnach o 
driongliannau dynamegol yn gyffredinol. Yn y pen draw, gall y rhaglen ymchwil lion roi 
gwell mewnwelediad o ddefnydd potensial triongliannau dynamegol mewn disgyrchiant 
cwantwm. Cynigir y thesis yma fel astudiaeth chwilogol i'r maes newydd yma. 

Mae'r symudiadau pum dimensiynol (k, I) a ddefnyddir yn algorithm Monte Carlo 
wedi'u profi i fod yn ergodig yn y bwlch o 5-maniffoldiau symhlygol cyfartal cyfuniadol. 
Gwneir gosodiad ynglyn a bodolaeth posib o ffin uchaf esbonnyddol ar y nifer o drionglian- 
nau o'r 5-sffer sy'n gyfartal gyfuniadol. Mae efelychiad Monte Carlo yn datgelu strwythur 
gwedd annistadl ac fe'u dadansoddir yn fanwl. Mae ymchwiliadau pellach yn ymdrin a 
natur geometrig a ffractal y triongliannau. Dilynnir hyn trwy ddangos nodweddion ffin 
y cyplyn gwan yn nhermau sfferau wedi'u stacio. Defnyddir ymresymiadau damcaniaeth 
graff syml i atgenhedlu a chyffredinoli canlyniad tra chyfarwydd yn nhopoleg cyfuniadol. 
Yn olaf cyflwynir astudiaeth cynhwysfawr o strwythurau hynod yn nhriongliannau dy- 
namegol. Mae'n cynnwys dealltwriaeth newydd o'u bodolaeth sydd yn ymddangos i fod 
yn gyson ag anfodolaeth o fertigau hynod mewn tri dimensiwn. Wrth gloi mae'r thesis yn 
rhoi golwg cyflawn o'r canlyniadau, dadl gyffredinol ac awgrymiadau o waith i ddyfod. 
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Chapter 1 



Introduction 



1.1 Quantum Gravity 

UR current understanding of gravity is based upon Einstein's general 
theory of relativity. General relativity has proved to be one of the great 
achievements in theoretical physics this century. It has predicted a num- 
ber of gravitational phenomena which have subsequently been verified 
experimentally to high degrees of precision. These effects include the advancement 
of Mercury's perihelion, gravitational redshift and gravitational lensing. Another effect 
currently under investigation is gravitational waves. The success of general relativity 
is merited by perfect agreement between theory and experiment at cosmological scales. 
However, general relativity is flawed in that it fails to describe the physics of spacetime 
singularities, such as those found in black holes and at the birth of the universe. In 
other words, we have a theory of gravity which is perfectly consistent with Nature at 
cosmological scales, yet fails for short distance scales. From this, one may cautiously infer 
that general relativity is a low energy limit of a more fundamental theory of gravity. 

Perhaps the most outstanding scientific achievement of the century is the development 
of quantum theory. The effects of quantum mechanics and quantum field theory have 
also been tested in numerous experiments to very high degrees of precision at atomic and 
sub- atomic scales. These tests include the photoelectric effect, electron diffraction, the 
Stern-Gerlach experiment and calculation of the Lamb shift, to name but a few. It is this 
overwhelming success that has led physicists to believe that Nature is ruled by quantum 
theory. Presently, the only fundamental force which still eludes a satisfactory quantum 
description is gravity. 

One of the most important advancements in the understanding of fundamental forces 
was made by Maxwell. He proved that two seemingly disparate forces (electricity and 
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magnetism) may be described by a single theory - electromagentism. More recently, it 
has been shown that the electromagnetic and weak forces may be unified into a single 
electroweak force. This was made possible by realising that such theories can be under- 
stood in terms of symmetry groups. It is hoped that unification may be taken a step 
further by including the strong force, resulting in a grand unified theory (GUT). Un- 
like the electroweak theory, a GUT is characterised by a single coupling constant. Of 
course, the final and ultimate unification would be of all four known fundamental forces, 
incorporating a quantum theory of gravity. Such a theory is commonly referred to as a 
theory-of-everything. 

For these reasons it is generally believed that there exists a quantum theory of gravity. 
Before going any further, we first ask ourselves an important question: what features 
do we expect quantum gravity to possess? First of all, gravitons, the excitations of the 
gravitational field, must be massless - corresponding to long range interactions. Secondly, 
the theory must reduce to general relativity in the low energy limit and to Newtonian 
gravity in the weak coupling limit. Finally, the graviton should be a spin 2 particle since 
the spacetime metric is a rank 2 tensor. 

Quantum gravity has been studied extensively by many research groups over the latter 
half of this century. Despite this fact, no satisfactory theory has yet been formulated. One 
of the main problems facing theorists is the fact that no experimental data from graviton 
scattering is (or will ever be) available, since the quantum effects of gravity are expected 
to emerge near the Planck scale, i.e. ~ 10 19 GeV. This situation is in stark contrast to 
the early days of quantum field theory, when theory was led by experiment. 

On that note, let us briefly mention some of the current avenues of research in the field 
of quantum gravity. For a comprehensive survey of the research programme as a whole, 
the reader is referred to a recent review by Rovelli Perhaps the most promising and 
popular candidate theory to date is superstrings; a theory in which particles (including 
gravitons/gravitinos) are described by the excitations of one dimensional objects in ten 
dimensional spacetime. Loop quantum gravity, which takes the canonical quantisation 
approach, has also been studied intensively over the last decade with encouraging results. 

Theories of discretised spacetime have become increasingly popular of late [0. A 
subset of these, generically called simplicial quantum gravity, was founded upon the Eu- 
clidean quantum gravity programme || and the pioneering work of Regge 0] . Simplicial 
quantum gravity encompasses two closely related and complementary theories, namely, 
quantum Regge calculus and dynamical triangulations - the latter being the subject of 
this thesis. Other candidate theories include non- commutative geometry, topological field 
theory, supergravity and twistors. 

The reconciliation of quantum theory and relativity has proved to be a formidable 
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challenge. Indeed, their seemingly incompatible philosophical underpinnings may only 
be overcome by the formulation of a radical new theory, rather than a quantisation pre- 
scription per se. It is possible (if not likely) that a consistent formulation of quantum 
gravity will force us to abandon our current understanding of the small scale structure 
of spacetime. These difficult issues have been addressed by Isham || |[, amongst many 
others. 

At this point, it is expedient to describe the many problems associated with the (naive) 
conventional field theoretic approach to quantum gravity. This will lead us naturally to 
a review of simplicial quantum gravity. 

1.1.1 Euclidean Quantum Gravity 

We shall now consider the path integral approach to quantum gravity and highlight the 
associated problems. In this formalism one defines a functional integral over field configu- 
rations of a spacetime manifold A4 with metric g^ v . If M. has a Lorentzian signature then 
the path integral is oscillatory and hence divergent, due to complex Boltzmann weights. 
Fortunately, it is possible to overcome this problem by Wick rotating time t clockwise 
by 90° in the complex plane to give the following transformation t — > —it. The resulting 
Euclidean path integral Ze is written as 



The measure includes a factor which prevents the overcounting of diffeomorphic field 
configurations. (Of course, to be precise, Z E is a Riemannian path integral as spacetime 
is not necessarily flat by any means.) In (^-dimensional Euclidean space the Einstein- 
Hilbert action Slg^, A, G] is defined (in natural units) as 



where R is the scalar curvature, A is the cosmological constant and G is Newton's gravi- 
tational constant. 

Having defined the Euclidean path integral Ze, we immediately encounter some serious 
set-backs. Firstly and perhaps most importantly, one can only assume that the signature 
may be rotated back to Lorentzian spacetime. Secondly, one expects the topology of 
spacetime to be dynamical in a quantum theory of gravity. Unfortunately, the topological 
classification of smooth (differentiable) 4-manifolds is still an outstanding problem in 
mathematics. The problems do not end here: the Einstein-Hilbert action Slg^, A, G\ is 
known to be unbounded from below. Perturbation theory, in this framework, leads to 
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yet more complications. The theory is perturbatively non-renormalisable since G has a 
negative mass dimension. It is possible to overcome this problem by introducing higher 
derivative terms into the action, but only at the cost of unitarity |7|. One may strive to 
avoid many of these problems by considering the Euclidean path integral approach in a 
non-perturbative framework - namely simplicial quantum gravity. 

1.2 Simplicial Quantum Gravity 

Simplicial quantum gravity is, in essence, a non-perturbative realisation of the Euclidean 
path integral over spacetime geometries. The key ingredient of this approach is lattice 
discretisation. This scheme brings the theory into the realm of statistical physics, where 
it can be studied numerically. Contact with the continuum may then be made at a 
critical point where the system becomes scale invariant. Traditionally, in such lattice field 
theories, the spacetime lattice serves as an inert background upon which the fields are 
defined. However, in this particular case, the lattice itself is dynamical, which furthermore 
must support a definition of curvature. As we shall see, this compels us to abandon the 
use of hypercubic lattices in favour of simplicial spaces. 

For the time being we shall offer a short resume of simplicial quantum gravity; details and 
precise definitions are to be found later. Simplicial spaces are constructed by connecting 
together a finite set of 'elementary building blocks' called simplices. A d-dimensional 
simplex (or <i-simplex) is, in a certain geometrical sense, the most fundamental bounded 
unit of Euclidean space R d . For example, a 0-simplex is a vertex, a 1-simplex is a link, a 
2-simplex is a triangle and a 3-simplex is a tetrahedron, and so on. In general, a simplex is 
formed by the product of its constituent sub simplices. For example, a triangle is formed 
by the product of three vertices. As a general rule, [d — l)-subsimplices are called faces 
and (d — 2)-subsimplices are called hinges. 

In this thesis, we shall only consider a special subset of simplicial spaces, known 
as simplicial manifolds. As with smooth varieties, simplicial manifolds also support a 
definition of curvature in terms of the rotation of parallel transported vectors around 
hinges; a crucial fact discovered by Regge ||. Using this result, one can formulate a 
discretised counterpart of the Einstein-Hilbert action. These constructs enable us to 
'discretise' the Euclidean path integral. This discretisation scheme, in effect, reduces Ze 
to a partition function Z 

Z=Y, exp(-£[T]), (1.3) 
TeT 

where S[T] is the discretised Einstein-Hilbert action. The partition function is a sum- 
mation over all triangulations T of a space of triangulations T, each weighted by its 
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corresponding Boltzmann factor. 



There are two main variants of simplicial quantum gravity which correspond to the two 
'degrees of freedom' of simplicial manifolds. These are the link lengths and the connec- 
tivity. The connectivity of a simplicial manifold is an abstract notion which refers to 
the different possible ways of interconnecting its constituent simplices. If the geometry is 
governed only by the link lengths, then the resulting theory is called quantum Regge calcu- 
lus, the quantum extension of Regge calculus. Conversely, geometries that are controlled 
solely by the connectivity correspond to a theory known as dynamical triangulations. In 
dynamical triangulations, simplices are treated as being equilateral and hence link lengths 
play no role in the dynamics. Unlike quantum Regge calculus, it is a regularisation since 
the discretisation of spacetime introduces an intrinsic cut-off a, corresponding to the link 
length. 

The goal of simplicial quantum gravity is to make contact with a physical continuum 
limit. This can only be done if the statistical ensemble loses scale dependence; a char- 
acteristic of continuous phase transitions. These critical points are searched for in the 
space of parameters using Monte Carlo simulation methods, which are discussed later in 
section [0| Before looking at Regge calculus, quantum Regge calculus and dynamical 
triangulations in more detail, let us first review some background mathematics. 



1.2.1 Combinatorial Topology 

The study of simplicial spaces and manifolds is known in mathematics as combinato- 
rial topology (or piecewise linear topology) || [|, 10, 11]. This well-established branch of 
mathematics partly forms the theoretical basis of dynamical triangulations. The follow- 
ing sections are devoted to clarifying our definitions and conventions on the subjects of 
topology and combinatorial topology, which will be upheld throughout the remainder of 
this thesis. Fortunately, these basic constructs are easily defined and explained without 
the need for unnecessary mathematical rigour. 



Topology and Manifolds 

Two topological spaces A and B are said to be homeomorphic if there exist a continuous 
map f:A>—>B and a continuous inverse map f~ x : B i— > A between points on A and points 
on B. A continuous map g:C^T> between topological spaces is called an embedding if it 
is a homeomorphism onto a subspace of T>. There are many kinds of topological spaces. 
We are concerned only with those that are locally homeomorphic to Euclidean space, i.e. 
manifolds. Throughout this thesis all manifolds are understood to be closed, compact, 
connected and orientable. 
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One goal of mathematics is to classify smooth manifolds (of a given dimension) ac- 
cording to their topology in terms of a finite set of topological invariants. This has been 
achieved in all dimensions except three and four [|T^]. The topology of 2- manifolds are clas- 
sified by a single integer invariant known as the genus h, which is essentially the number 
of handles. Alternatively, 2-manifolds may be classed in terms of the Euler characteristic 
X, which is related to the genus by x = 2 — 2h. In this thesis we shall only encounter the 
three most fundamental topologies: d-balls B d , d-spheres S d and d-tori T d ,Q. 

Simplices and Triangulations 

In this section, the reader may assume that the subscript i takes integer values from 
to d inclusively. Let R n be an n-dimensional Euclidean space. Formally, a ci-simplex a d 
may be defined as the convex hull formed by a set of d + 1 affinely independent points in 
R n , where n > d. From here on, a 'simplex' is assumed to be <i-dimensional. A simplex 
is composed of d+1 Cj + \ constituent j-dimensional subsimplices, for j = 0, . . . , d — 1. An 
i- (sub) simplex a 1 of a simplex a d is defined as the i-simplex formed from a subset of % + 1 
points of the d + 1 points of a d . As a general convention in this thesis, a '(sub) simplex' 
is taken to mean either a subsimplex or a simplex. 

A simplicial complex K is a topological space constructed from a finite number of 
fc-simplices with the property that any two simplices either intersect at a common sub- 
simplex or are disjoint. The dimension of K is that of its highest dimensional constituent 
fc-simplex. Consider an i- (sub) simplex o % of a simplicial complex K. The star of o % , 
written St(a l ), is defined as the set of fc-simplices in K which intersect o~ % . 

In this thesis, we deal exclusively with simplicial manifolds, which form a special 
subset of simplicial complexes. By this we mean that all manifolds are complexes, but 
not all complexes are manifolds. A rf-dimensional simplicial manifold (or triangulation) 
T is a simplicial complex which satisfies the added contraints that demand all faces to 
be pairwise connected and that every star of T is homeomorphic to B d . For a simplicial 
manifold, the local volume n(a l ) of a 1 is defined as the number of simplices in St(a l ). When 
referring to the local volumes of i- (sub) simplices in general, we shall use the symbol rij. 
To any given simplicial manifold T, one may associate an f -vector 

f=(N ,...,N d ), (1.4) 

where iVj is the total number of i- (sub) simplices of T. In almost all cases the /-vector 
cannot uniquely identify a triangulation. This can only be done with knowledge of the 
connectivity. Given that simplicial manifolds are combinatorial in nature, one may expect 

1 A d-sphere is homeomorphic to the boundary of a (d+ l)-ball, and a d-torus is homeomorphic to the 
product of d 1-spheres. 
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certain relations between the components of an /-vector. This is indeed the case. It 
turns out that there are two sets of equations known as the Euler-Poincare and Dehn- 
Sommerville relations. 

Euler-Poincare and Dehn-Sommerville Relations 

Euler's theorem is a classic result of combinatorial topology. It relates the Euler charac- 
teristic x to the total number of vertices V, links L and faces F of polyhedra. 

X = V - L + F = 2 (1.5) 

The generalisation of this result to rf-dimensional simplicial manifolds is known as the 
Euler-Poincare relation. 

x = EHm (i.6) 

8=0 

The Euler characteristic x is zero f° r °dd dimensional closed manifolds and equal to two 
for even dimensional <i-spheres. 

The Dehn-Sommerville relations are a set of d + 1 equations linear in the number jVj 
of z-(sub)simplices of rf-dimensional simplicial manifolds []I3| . They reduce to (^2) 
linearly independent relations in even (odd) dimensions. Collectively, they are written as 

Ni = E(-l)^' j+1 C i+1 Nj, (1.7) 

j=i 

where i = 0, . . . , d. These equations are necessary but not sufficient for a simplicial 
complex to be a simplicial manifold. In other words, all simplicial manifolds satisfy 



equations (|1.7| ), but not all solutions of (|1.7| ) correspond to simplicial manifolds 



1.2.2 Regge Calculus 

Regge calculus was formulated by Regge in 1961 as a coordinate independent discretisa- 
tion of general relativity M. He showed that one can approximate smooth (spacetime) 
manifolds with piecewise-linear (simplicial) manifolds. His methods were used as an ap- 
proximation of general relativity and have since been used for numerical and computa- 



tional work 14, 15 



As mentioned earlier, curvature has a natural definition on simplicial manifolds in 
terms of parallel transportation of a vector around a hinge. For the moment we shall 
give a heuristic explanation of why curvature is only defined at hinges. Consider the 
embedding of an arbitrary <i-surface S in R n , where n > d. If S cannot be embedded in 
R d , then it may be deemed to be curved. For simplicity, let us apply this notion to two 
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dimensional simplicial manifolds. Clearly, the stars of triangles and links can always be 
embedded in R 2 and are therefore flat. These structures are illustrated in figure |1 . 1| . 






Figure 1.1: Embeddings of two dimensional stars in 



From figure ET , it is clear that a vertex star cannot in general be embedded in R 2 . This 
shows that curvature is defined at vertices (hinges). Of course, the embedding of a vertex 
star in R 2 is possible provided the link lengths satisfy certain conditions. Although difficult 
to visualise, these ideas can be generalised for hinges of higher dimensional simplicial 
manifolds. This brings us to the concept of deficit angles. 

In its simplest form, the Einstein-Hilbert action Sr of a <i-manifold with metric g^ v 
and scalar curvature R is given byQ 



Sr = J d d £(detg^R. (U 



In his seminal paper, Regge showed that the curvature at a hinge h is proportional to its 
associated deficit angle 5h 

5 h = 2n- ]T 6(a d ), (1.9) 
<r d e st(h) 

where 8(a d ) is the dihedral angle subtended by simplex o d at h. The deficit angle 5h is the 
rotation angle of a vector parallel transported around h. By summing over all Nd-2 hinges 
of a simplicial manifold T we obtain an expression for the 'integrated' scalar curvature. 
This gives us the discretised Regge action Sn egge 

Sr > S Regge = ^bVh, (1.10) 



2 Here, for simplicity, the coefficient has been rescaled 'away' and the cosmological constant A is set 
to zero. 
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where Vh is the volume of h. An attractive feature of Sn egge is that it can be shown to 
be equivalent to the continuum action Sr in the limit where the link lengths are reduced 
to zero and their number is taken to infinity. In fact, Sji egge can also be derived from the 
continuum action HTB1. 



Clearly, all the dihedral angles and hinge volumes of a simplicial manifold are determined 
purely by the link lengths, which are subject to the triangle inequality or its higher 
dimensional generalisations. In other words, the curvature of the simplicial manifold is 
governed solely by the link lengths. It is useful at this point to illustrate the connection 
between deficit angles and curvature. Again, for convenience, we consider two dimensional 



simplices. In figure [T7| the top pair of pictures show the stars of hinges hi and hi projected 
onto the plane R 2 . 






Figure 1.2: Curvature and deficit angles of two dimensional hinge stars. 

The top-left picture shows hinge hi with positive deficit angle. By conjoining vertices a 
and g and embedding the structure in R 3 we see that the surface has positive curvature 
(bottom-left). Conversely, the top-right picture shows hinge /12 with negative deficit angle. 
This configuration forms a saddle-type structure with negative curvature when embedded 
in R 3 (bottom-right). If the deficit angle is zero then the structure is flat. 

Classical problems in general relativity may be solved by varying the link lengths of a 
simplicial manifold to find the extremum of the action. This correponds to the solution 
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of Einstein's equation in the limit of infinitesimally fine triangulations. We now move 
on to briefly mention the quantum extension of Regge Calculus. For a comprehensive 
(though dated) bibliography of the subject, the reader is referred to an article by Tuckey 
and Williams [fRj]. 



1.2.3 Quantum Regge Calculus 

Quantum Regge calculus is a realisation of the Euclidean path integral Ze in terms of 
discretised spacetime |[L7| . As with Regge calculus, the dynamical variables are the link 
lengths of a simplicial manifold with fixed connectivity. The partition function is defined 
such that the measure includes a summation over all possible simplicial manifolds which 
satisfy the triangular inequality, or its higher dimensional generalisations. The continuum 
measure transforms as follows. 

Lvfe^fnduft,...,w (lid 

So, the continuum measure of Z E is replaced by a product of Ni integrals over all possible 
link lengths U from to oo. Hence, quantum Regge calculus is not a regularisation scheme; 
the triangulated lattice has no cut-off. Included in the measure is a non-trivial Jacobian 
J(li, ■ ■ ■ , ijvj, whose form is not generally known. The partition function Z Regge is thus 
defined as 

Z E — >Z Regge = Y[dl i J(li,...,l Nl )exTp(-S Re g ge ). (1.12) 

Jo i=i 

Numerical computation of Z Regge is possible only by introducing long and short length 
cut-offs to the integral. Contact with the continuum may be reached at a suitable scaling 
limit where the volume of the simplicial manifold is taken to infinity. One hopes, on the 
grounds of universality, that such a limit would be insensitive to discretisation details, 
such as the choice of triangulation connectivity. 

Quantum Regge calculus is not without its problems. Perhaps the most pressing diffi- 
culty relates to an ambiguity in the choice of measure. In addition to this, the integration 
over link lengths is known to overcount configurations, which may possibly be circum- 
vented by introducing a Fadeev- Popov factor into the measure. In fact, a recent study 
has concluded that the theory 'fails' in four dimensions 



1.2.4 Random Surfaces and Dynamical Triangulations 

The origin of dynamical triangulations can be traced back to the application of random 
surfaces to string theory during the early 1980s and the seminal work of Weingarten |[19|| . 
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String theory may be viewed as being geometrical in nature since propagating strings 
sweep out world-sheets in their embedding space. It was Polyakov who first proposed 
the quantisation of (closed orientable) bosonic strings by defining a path integral over 
surfaces with internal metric g a b, embedded in D-dimensional space One may use 
this formalism as a quantisation of gravity, where the world-sheets are interpreted as 
spacetime manifolds and the D-dimensional embedding space become D scalar fields 0. 
In this framework, the Polyakov action is written as 

S[g ab , ( p} = J2 J d 2 £(det^)V^W + 77 / d 2 £, (1.13) 

i=l 

where rj is a constant. This is equivalent to pure two dimensional quantum gravity when 
D — 0. Here, D is interpreted as the central charge of the conformal matter coupled 
to gravity. This continuum formulation of two dimensional quantum gravity was solved 
analytically by Knizhnik et al. using Liouville theory |pl| . 

One can solve dynamical triangulations analytically using matrix models |22[|. In this 
case, one defines a path integral Zm over Hermitian matrix fields My of dimension N. 

Z M = J dil%exp(-S[il%]) (1.14) 

The action S[Mij,g] is given by 



S[M l3 , g] = \M i3 M* 3 + -l=M l3 M ]k M kl} (1.15) 
^ v JM 



where g is a coupling constant. A perturbative expansion in powers of g yields all pos- 
sible 3 graphs. There exists a direct equivalence between dual 3 graphs and two di- 
mensional triangulations, as triangles have three faces. In other words, there exists a 
direct correspondence between Zm and the partition function of two dimensional dynam- 
ical triangulations. The analytic results of this formalism has been shown to agree with 
those of Liouville theory for c < 1, where c is the central charge 0. It is this success of 
dynamical triangulations which has encouraged physicists to consider higher dimensional 
generalisations. 

Two Dimensional Dynamical Triangulations 

Let us now consider dynamical triangulations in more detail by first deriving the action 



in two dimensions. From section |1.2.2| , we know that the deficit angle of a hinge h is given 
by 

6 h = 2n- £ 6{a 2 ). (1.16) 

<T 2 G St(h) 

3 Liouville theory with c > 1 matter 'breaks down'. 
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For dynamical triangulations, equation ( |1.16| ) requires modification since all simplices are 



equilateral. Specifically, the dihedral angle between two faces of a triangle is a constant 
2 . As a result, the deficit angle of a hinge h depends solely on its local volume n(h). 
Therefore, we have 

6 h = 2n -n(h)& 2 . (1.17) 

Given that Q 2 — 7r/3, it is clear that h is flat if n(h) = 6. As we shall now see, deficit 
angles are not particularly relevant to the two dimensional case. The Gauss-Bonnet 
theorem states that the integrated scalar curvature of a closed two dimensional manifold 
M. is a topological invariant, related to the Euler characteristic \- 



[ d^R = 4n X (1.18) 
Jm 



Using this theorem, one can show that the two dimensional Einstein-Hilbert action is 
given by 

%^,M,G] = M-§, (1-19) 

where A is the area of A4 and \i plays the role of the cosmological constant. This may be 
discretised by realising that the area of a simplicial 2- manifold is proportional to N 2 . 

S[T]=fiN 2 -^ (1.20) 

The partition function of dynamical triangulations is a summation over all triangula- 
tions T of a particular space of triangulations T. 

Z= E exp(-5[T]) (1.21) 

As one might expect, there exist triangulations of T related by symmetry. The measure 
is a sum over physically inequivalent triangulations, therefore those related by symmetry 
must not be overcounted. The issues relating to the meaning of T are discussed in depth 



in section 2.1.2 



1.3 Monte Carlo Simulations 

The numerical study of dynamical triangulations basically entails the detection of critical 
points in the space of couplings. Ultimately, this process involves calculating averages of 
thermodynamic variables. In dynamical triangulations the average value of an observable 
O is given by 

(OJ^E Oexp(-S[T}). (1.22) 
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The (grand canonical) partition function Z obviously cannot be evaluated in practice. 
Instead, we consider the canonical partition function Z c 



where T e T{N ( j) represents the subspace of triangulations T of T with (fixed) volume N^. 
(This sum can only be evaluated by making certain modifications. These are discussed 
later in chapters ||] and y.) Although Z c is a finite sum, the number of terms is, in 
general, very large. For example, even a 16 2 Ising lattice has 2 256 « 10 77 possible spin 
configurations! This scale of calculation is obviously far beyond current computational 
capabilities. For most cases, the exact calculation of the canonical partition function 
is therefore effectively impossible. One can, however, calculate an approximation of Z c 
by sampling only a small proportion of the total number of configurations. Random 
sampling of configurations leads to very poor approximations because, typically, very few 
configurations contribute significantly to Z c . One can overcome this problem using the 
importance sampling methods of Monte Carlo simulations. 

Importance Sampling and Markov Chains 

The sampling of 'important' configurations is achieved by generating a Markov chain. This 
is a sequence of configurations, propagated by performing a series of successive local moves 
on the ensemble. This means that any configuration in the Markov chain depends only on 
the preceding one. Markov chains are therefore, by their very nature, highly correlated. 
This fact must be taken into account when generating statistically independent data. The 
average value of an observable will converge to the thermal average as the Markov chain 
grows in length. 

Local moves refer to algorithms which can transform one configuration into another. 
They are 'local' in the sense that only a small portion of the configuration is altered, 
whose size is independent of the ensemble volume. For example, in the Ising model, the 
local move is a single spin flip. In the case of dynamical triangulations, the local moves 
are much more complicated. Unlike the Ising model, the transformation algorithm is 
dimension dependent. These moves will be reviewed later (see chapter |2| for details). To 
ensure an accurate approximation of a thermal average, the local moves must be ergodic 
in the space of configurations. In other words, all configurations must be accessible to 
the Markov chain. Given that V(T — > T") is the normalised probability of transforming 
configuration T into T', then we insist that 



E 





T £ T(N d ) 



Y^V{T -> T') = 1. 
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This condition ensures that the Markov chain cannot come to an end. The Markov chain 
must also satisfy another condition known as detailed balance. 

exp(-S[T]) V{T -> T') = exp(-S[T'}) ViT' -> T) (1.25) 

This ensures that configurations are chosen with the correct probability distribution of 
exp(— S[T]), provided the local moves are ergodic. 



Metropolis Algorithm 

A configuration may only be added to the Markov chain if accepted by an update accep- 
tance algorithm. The Metropolis algorithm is most often used. Consider the transforma- 
tion of configuration T into T' . If the change in action is given by 

AS = S[V] - S[T], (1.26) 



then the probability V of acceptance is as follows. 

V(T rpi\ ( exp(-AS) if AS>0 

' 1 j \ 1 if AS < 

In other words, the update is automatically accepted providing the change in action is 
not positive. In practice, acceptance is determined in the following manner. If a pseudo- 
random number is found to be greater than exp(— AS) then the update is accepted, 
otherwise it is rejected. 



More often than not, the first few configurations of a Markov chain will be statistically 
unimportant. Therefore, in general, an ensemble requires a finite thermalisation time r 
to evolve towards statistically important configurations. Any measurements taken prior 
to time r will result in inaccurate estimates of thermal averages. Of course, successive 
configurations of Markov chains are highly correlated. Hence, in order to generate sta- 
tistically independent data, one must allow a finite time for the ensemble to decorrelate. 
This is known as the autocorrelation time. The statistical treatment of independent and 
correlated data is explained in appendices BT1 and |B.2j respectively. 

A well-known drawback of the Monte Carlo method is an effect known as critical 
slowing down. Near a critical point, one finds long range correlations in the ensemble. 
Under these circumstances, the number of local moves required to generate a statisti- 
cally independent configuration grows very rapidly. In other words, the algorithm loses 
efficiency 
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1.4 Dynamical Triangulations 

So far, we have reviewed the historical development and certain technical aspects of 
simplicial quantum gravity and dynamical triangulations. The research of dynamical 
triangulations in 'higher dimensions' (d > 2) stems from the success of the theory in 
two dimensions. Let us now consider these extensions by first deriving the discretised 
Einstein-Hilbert action for arbitrary dimension and writing down the partition function. 



1.4.1 Discretisation 

The continuum Einstein-Hilbert action S[g^ V) A, G] for a d- dimensional manifold M. is 
defined as follows. 

A, G] = ^ J m d^ (det (2A - R) (1.27) 

We shall now formulate its discrete analogue S[T], by considering each term of S[g^ u , A, G] 
individually. The cosmological constant term is obviously proportional to the volume of 
M. and discretises to give 

/ d d £(det^ -^V d N d , (1.28) 

JM 

where V d is the volume of a simplex. The volume of an equilateral simplex with link 
length a is given by |23 



a d fd+r 



■2 



Given a d- dimensional simplicial manifold T, the deficit angle of a hinge h in T is 
defined as 

6 h = 2n-n(h)Q d , (1.30) 

where 9^ is the dihedral angle subtended between two faces of a simplex. In arbitrary 
dimension [53|, this is given by 

O d = cos~ 1 (\). (1.31) 



Therefore, the integrated curvature at a single hinge h is 



^d d C (det g^R — > 2V d „ 2 5 h (1.32) 
— > 2V d - 2 (2ix -n(h)Q d ). (1.33) 
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By summing over all N d _ 2 hinges of T and using the fact that each simplex has C^-i 
hinges: 

£n(/i) = d+1 C d ^N d (1.34) 

heT 

= ^d{d+l)N d , (1.35) 

we deduce that 

/ d d £(detg^R — > 2\/ d _ 2 £ (2tt - n(/i)9 d ) (1.36) 

— > 2y d _ 2 ^2 7 riV (i _ 2 -^(d+l)iV d e^. (1.37) 
By combining and expanding ( |1.28| ) and (|1.37|) , we get 

S[g^, A, G\ — > 5[T] = (2AV d iV d - 4vrV^_ 2 iV d _ 2 + d (d + l)V d ^N d ®d) ■ (1.38) 

107TCt 

This may be rewritten simply as 

S[T, K d - 2 ] = K d N d - K d _ 2 N d - 2 , (1.39) 

where the constants n d and K d _ 2 are given by 

2AV d + d(d+l)V d . 2 Q d 
Kd = 16^ (L40) 

= ■ (1.41) 

The exact form of K d can vary slightly depending on how the cosmological constant is 
defined in the continuum Einstein-Hilbert action^. This thesis is primarily concerned 
with dynamical triangulations of five dimensional simplicial manifolds. In this case the 
discretised Einstein-Hilbert action is obviously given by 

S[T,k 3 ] =k 5 N 5 -k 3 N 3 . (1.42) 

Given the discretised Einstein-Hilbert action, one can formally write down the grand 
canonical partition function as 

Z E — > Z(K d „ 2 ,K d ) = ex P(-S[T, K d - 2 \) (1-43) 

TeT 

= exp(K d _ 2 N d _ 2 - n d N d ). (1.44) 



The factor of 2 in (1.32) is occasionally omitted, as in 
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The canonical partition function is defined as 

Z c (K d _ 2 ,N d )= J2 exp(K d _ 2 N d _ 2 - n d N d ). (1.45) 

TeT(N d ) 

Having explained the discretisation scheme in some detail, let us now move on to 
outline the general picture of dynamical triangulations in dimensions greater than two, 



which was first considered by Godfrey and Gross [25 



1.4.2 d > 2 Dynamical Triangulations 

The first extension to three dimensional dynamical triangulations was made in the early 
1990s p6| , |27J. This work was soon followed by preliminary investigations of the four 



dimensional theory E?J. Over the years, dynamical triangulations has become rela- 
tively well-established; a fact reflected by a number of review articles [30, 31, 33] and 
books |12|, |4 . 



Three and four dimensional dynamical triangulations have two distinct phases which 
are fundamentally very similar in nature. The weak coupling phases (small G) correspond 
to elongated triangulations, which resemble branched polymers, and the strong coupling 
phases (large G) are characterised by 'crumpled' triangulations with large (possibly infi- 
nite) fractal dimension. The phase transitions in both cases appear to be of first order, 
which is indicated by bistable behaviour, and hysteresis effects in three dimensions. This 
effectively eliminates the possibility of defining a continuum limit. 

A fundamental distinction between three and four dimensions relates to the phe- 
nomenon of singular vertices. These vertices with unusually large local volumes are 
present in the strong coupling phase of four dimensions, but not in three dimensions. 
Studies have shown that they are significant in terms of the phase transition in four di- 
mensions, though their role is not yet fully understood ||35|| . This subject is dealt with 
in chapter |6|. Other advances include the analytic study of stacked spheres, degenerate 
triangulations and simplicial manifolds with boundaries. 

In recent years, much effort has been devoted to developing 'modified models' of 
dynamical triangulations. Modifications include matter fields or new terms in the action. 
These advances are motivated by the hope that an expanded phase diagram may reveal 
extra phase structure which could lead to a continuum limit. It is in this spirit that we 
consider dynamical triangulations in five dimensions, as the most general action (linear 
in the components of the /-vector) has two independent coupling constants. It is hoped 
that this new perspective will yield a deeper understanding of dynamical triangulations 
in general. Before exploring the model in any detail, one must first show that it is well- 
defined. This involves proving the ergodicity of the five dimensional (k, I) moves and the 
existence of an exponential bound (see chapters and |3| respectively) . 



Chapter 2 

Local Moves and Ergodicity 



2.1 Local Moves 




J ONTE Carlo simulations have a wide scope of application in physics. 
These methods have so far been central to the research of dynamical 
triangulations. This chapter deals with certain technical aspects of 

i the update algorithm in the context of five dimensional dynamical 



triangulations. To be precise, we shall clarify the meaning of the space of triangulations 
T, describe the known local moves and attempt to prove their ergodicity in T. 

The canonical partition function Z c (k, 3 ,N 5 ) cannot be computed in practice, as all 
known local moves change the triangulation volume N 5 (see later). To ensure ergodicity, 
one must therefore allow N 5 to fluctuate around a chosen target volume 7V|. These fluctu- 
ations may be controlled by adding an extra Gaussian term to the action, thus redefining 
Z c (k 3 , N 5 ) as 

Z'(« 3 ,iV|,7)= E exp(/ t3 iV3-/ t5 iV 5 - 7 (iV 5 -iV*) 2 ), (2.1) 
TeT 

where 7 controls the degree of fluctuations in JV5. Let us now consider the subject of 
ergodicity in more detail. 

2.1.1 Ergodicity 

Local moves are ergodic if every configuration of a statistical ensemble is accessible in 
principle. This property is necessary for the accurate approximation of thermal averages. 
The matter of proving ergodicity is trivial in many statistical models. For example, in the 
Ising model the local move is a simple spin flip from | to I or vice versa. It is obvious that 
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all configurations are eventually accessible, providing that spins are chosen at random. 
The flip local move is therefore manifestly ergodic. 

Now consider the issue of ergodicity with regard to dynamical triangulations, beginning 
with the trivial case of one dimension, in which the only closed unbounded simplicial 
manifold is the circle S 1 . The space of triangulations of S* 1 is, of course, identical to 
the set of polygons. For a given volume Ni there is only one possible triangulation, 
up to an automorphic relabelling of vertices (or links). The simplest local move which 
transforms S* 1 triangulations entails exchanging one link for two, or vice versa. It is 
not even necessary to perform the moves on random simplices due to the Z Nl rotational 
symmetry. The moves are obviously ergodic in the space of S* 1 triangulations. The one 
dimensional case is trivial because each simplex (link) has only two faces. Given that all 
faces of simplicial manifolds are pairwise connected, this results in a unique triangulation 
of any given volume. The trivial one dimensional case may be generalised to higher 
dimensions, in which the local moves become much more complicated. In general, a local 
move performed on a (/-dimensional simplicial manifold T involves the retriangulation of a 
small B d submanifold of T. Clearly, proving the ergodicity of local moves is a non-trivial 
matter for dimensions greater than one. 

This brings us to an important aspect of the whole issue of ergodicity in dynamical 
triangulations. That is, over what space of triangulations are the moves required to 
be ergodic? For an Ising lattice of iV 2 spins, the space of possible configurations is an 
unambiguous concept. It is simply the 2 n2 possible combinations of spins. There is no 
overcounting of symmetrical configurations because each configuration corresponds to a 
real physical state of the Ising lattice. 

2.1.2 Combinatorial Equivalence 

For dynamical triangulations, the space of triangulations T requires precise and careful 
definition. Ideally, one hopes that 



is true in the continuum limit. As in the Ising model, there are triangulations which are 
related by automorphisms. These are equivalent by the relabelling of (sub) simplices and 
correspond to the same metric in the continuum. One can prevent the overcounting of 
symmetrical triangulations by including a factor that is equal to the inverse of the order of 
the automorphism group C T (of triangulation T). It is expected that this factor becomes 
less important for larger triangulations, as fewer are related by symmetry. Let us now 
formally define T based on this physical reasoning. 




(2.2) 
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Two simplicial manifolds T and T' are isomorphic if there exists a one-to-one corre- 
spondence between their subsimplices which preserves their connectivity. Such a mapping 
is denoted as T ~ T'. If no such map exists between T and T', then the triangulations 
are said to be distinct. Triangulation T' is a subdivision of T if every simplex of T' is a 
simplex (or part of a simplex) of T. Two simplicial manifolds are combinatorially equiv- 
alent if they have isomorphic subdivisions Pi This is written as T « T". Combinatorial 
equivalence is sometimes called a piecewise linear homeomorphism. Let C represent the 
space of combinatorially equivalent triangulations. 

The partition function is (in our case) defined for fixed topology, therefore the mea- 
sure is a summation over distinct homeomorphic triangulations. From our definitions, 
it is evident that T ~ T' then T « T' . In other words, isomorphic triangulations are 
combinatorially equivalent. However, the opposite is not true in general. Combinato- 
rially equivalent triangulations are homeomorphic. The conjecture that homeomorphic 
triangulations are combinatorially equivalent is known as the Hauptvermutung^. In other 
words, it states that homeomorphic triangulations have isomorphic subdivisions. It has 



been shown to be true for simplicial manifolds of dimension two |36[ and three |37|, but 



false for d > 4 P5fl . Let the space of homeomorphic triangulations be denoted by 
H. When the Hauptvermutung is true, then C C H and the measure is equivalent to 
a summation over combinatorially equivalent triangulations modulo automorphisms. In 
general C is comprised of a number of distinct automorphism classes. However, in the 
five dimensional case, the Hauptvermutung is false, therefore C C H. In other words, the 
measure is not equivalent to a summation over combinatorially equivalent triangulations 
modulo automorphisms. 



2.2 Alexander Moves 

The Alexander moves were developed by Alexander in 1929, who called them simple 
transformations |40|| . He showed that any complex could be transformed into any other 



equivalent complex using these topology preserving transformations. Alexander defined 
two complexes K and L to be equivalent if there exist partitions K' and V respectively 
such that K' and V are congruent. In his now classic paper, Alexander used the 'algebra of 
complexes' to formalise his simple transformations. We shall now review the background 
mathematics leading to explanations of the above italicised terms. 



the mathematical literature the term 'combinatorial equivalent' often has the same meaning as 
'isomorphic'. This unfortunate situation can cause confusion! 
2 German for 'fundamental conjecture'. 
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2.2.1 Algebra of Complexes 

A simplex o~ d can be expressed as the product)^ of its d + 1 constituent vertices of, for 
i — 1, . . . , d + 1. The set of vertices define the point-set of a simplex. Algebraically, a d is 
defined as follows. 

a d = al..a° d+1 (2.3) 

The product of an 2-simplex a % by a j-simplex o" J is an (i + j + l)-simplex a l a^ , providing 
that a 1 and a- 7 do not intersect, i.e. a % R = 0. A complex K is a simplicial complex 
with the added constraint that each simplex in K is of the same dimension. In algebraic 
terms, one can express a complex K of n simplices as the sum of its constituent simplices. 

n 

K = J2af (2.4) 

In general, K can be decomposed into a number of subcomplexes. For example, K may be 
expressed as the sum of two non-intersecting subcomplexes P and Q, such that K = P+Q. 
In general, the boundary of a complex K is denoted by K. Every subsimplex of a complex 
has an associated complement. The complement J of an i-subsimplex a 1 of a complex K 
can be defined as the set of subsimplices such that Ja % is a subcomplex of K. By this 
definition, it is clear that St(a l ) = Jcr\ One may then define the remainder subcomplex 
QoiK as K -Jo' 1 . 

A complex K' is a partition of a complex if their point-sets coincide and if each 
simplex of K' is a simplex or part of a simplex of K. Two complexes K and L (of the same 
dimension) are congruent if there is a one-to-one correspondence between the vertices of K 
and the vertices of L. Following the convention set by Gross and Varsted [III], equivalent 



complexes will be renamed Alexander equivalent. We are now ready to review the simple 
transformations, or Alexander moves as they are now commonly called. 

Simple Transformations 

Consider an i- (sub) simplex a 1 of a complex K with complement such that the remainder 
subcomplex is Q. 

K = a 1 Pi + Q (2.5) 
The set of d + 1 Alexander moves can be written collectively as 

K — > K' = (7° a 1 Pi + Q, (2.6) 



5 Multiplication is associative, commutative and distributive. 
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where i — 0, . . . , d, vertex o° ^ K and o' l Pi is the complement of o°. One may trivially 
define the inverse Alexander moves as 

K' — > K = a l Pi + Q. (2.7) 

Alexander proved that his set of d + 1 transformations are in fact all equivalent to the 
move with i = 1. This corresponds to a link, say the product of vertices o° and o°, with 
complement P. 

K = (o°o 2 °)P + Q <— ► K' = o°(o° + o 2 °)P + Q (2.8) 
In other words, the remaining d moves do not generate any additional complexes. Fig- 



ure [TT] illustrates the two dimensional Alexander move, where the link ac has complement 
b + d and x is a new vertex. For this particular case one may rewrite the transformation 
explicitly as 

ac{b + d)< > x{a + c){b + d). (2.9) 




2.3 (k, I) Moves and Dual Graphs 

The (k,l) moves were discovered more recently by Pachner 03 1, who called them 
elementary shellingsf^. These topology preserving local moves transform k simplices into 
I simplices. The methods used by Pachner to develop his transformations are based on 
advanced principles of combinatorial topology which are beyond the scope of this thesis. 
In this section we merely explain the transformation algorithm in general terms. 

General Description 

The set of (k,l) moves in <i-dimensions may be regarded as the set of all possible decom- 
positions of a minimal rf-sphere into two disconnected rf-balls. By this very definition, it is 



J They are sometimes referred to as the Pachner moves |44| or stellar subdivisions |4q, 46 



CHAPTER 2. LOCAL MOVES AND ERGODICITY 



23 



obvious that the moves preserve topology. In d- dimensions, the (k, I) moves must satisfy 

k + l = d + 2, (2.10) 

since minimal (i-spheres are composed of d + 2 simplices. As an illustrative example, 
the two dimensional (k, I) moves are represented by all the possible decompositions of a 
minimal 2-sphere, which is equivalent to the surface of a tetrahedron, see figure [2.2|. 







(2,2) 






(2,2) 





Figure 2.2: Two dimensional (k, I) moves as the decompositions of a minimal 2-sphere. 



2.3.1 Algebraic Formulation 

Let us now formalise the (k, I) moves in terms of Alexander's algebra of complexes. A 
(i-dimensional (k, I) move may be equivalently considered as the replacement of an i- 
(sub)simplex of a complex K with its dual (d — z)-(sub)simplex, where i = 0, . . . ,d. 

One can define the (k, I) moves in mathematical form for arbitrary k and d. Let us 
first consider the case where k — 1 and the complex K is written in terms of a simplex 
a d and the remainder complex Q. 



K = a d + Q 



(2.11) 



The (1,0?+ 1) move is defined as 



K — > K' = a° m P + Q, 



(2.12) 
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where <r° ^ K and Po is the complement of cr°. The complement Pq is the subcomplex 
surface of a d . This is commonly known as the vertex insertion move. Its inverse, the vertex 
deletion move, is trivially defined as 

K' = a°P + Q^>K. (2.13) 

Now consider the remaining moves with k = 2, . . . , d. Consider a (d — k + l)-subsimplex 
a d-k+i Q £ a com p} ex ^ with complement P^-fc+i- 

= <r d - fc+1 P d _ fc+1 + Q (2.14) 

A (k,l) transformation is given by 

K — > K' = a d - l+1 P d - l+1 + Q, (2.15) 
where <j d ~ l+l is a (d— Z + l)-subsimplex of K with complement P^-z+i- As an example, the 



transformations corresponding to the two dimensional (k, I) moves pictured in figure \2?2 
are given by 

acd < — > (ac + ad + cd)b (2-16) 
ac(b + d) < — > (a + c)bd. (2.17) 

2.3.2 Dual Graphs 

As we have seen, the two dimensional (k, I) moves are easily drawn and visualised. For 
higher dimensions, it is more convenient to transform the moves into their dual counter- 
parts. This allows us to draw moves of dimension greater than two on paper. A dual 
(k, I) graph is created by replacing each simplex with a node and each (d — l)-dimensional 
face by a edge between two nodes. A free edge (see figure |2~3| ) is the dual of a (d — 1)- 
sphere face. We shall now outline the general prescription for constructing dual graphs 
for arbitrary dimension. 

• Determine the initial (k) and final (I) number of simplices. Then draw their repre- 
sentative nodes - say in the shape of a polygon. 

• Draw edges connecting every pair of nodes. 

• Draw extra free edges such that each node has d+1 attached edges. This corresponds 
to each simplex having d+1 faces. 

• Check that no two nodes are connected by more than one edge. This relates to 
the manifold condition that impose simplices to be pairwise connected. The five 



dimensional dual (k,l) moves are shown in figure [T3 
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From this duality, it is clear that an i- (sub) simplex of a (^-dimensional triangulation is 
represented by a set of d — i + 1 interconnected nodes of a dual graph. This relationship 
allows us to count the change in iVj due to a (k, I) move by calculating the change in 
iV*f , for i = 0,...,d. 

ANi = AN^f (2.18) 

jSfduai - g number of different combinations of d — i + 1 nodes in the dual graph. The 
initial (final) dual graphs will have a total k (I) nodes respectively. 

ANi = N^f' 1 - Nf_f k (2.19) 

In other words we want to know the number of combinations of d — i + 1 from k and from 
I. This allows us to write down a very useful equation which gives ANi for any (k, I) move 
in any dimension d. 

AN, = l C d „ i+1 - k C d _ l+1 (2.20) 

Note that the equation is antisymmetric under the exchange of k and I as one would 
expect - the (I, k) move is the inverse of the (k, I) move. Equation ( |2.20|) is used in later 
chapters. Having formally defined the local moves, let us now consider their ergodicity. 



2.3.3 Ergodicity of (k, I) Moves 

In this chapter we are concerned only with proving the ergodicity of the five dimensional 
(k, I) moves. It is felt that a general proof for arbitrary dimension may be rather ambitious 
and unnecessary for our purposes. Ideally, we want to show that the moves are ergodic 
in the space of combinatorially equivalent simplicial 5-manifolds. 

Alexander proved that his simple transformations are ergodic in the space of combi- 
natorially equivalent simplicial manifolds. He achieved this by showing that two simpli- 
cial manifolds are Alexander equivalent if and only if they are combinatorially equiva- 
lent EDI, (TOj. We shall attempt to construct a proof of ergodicity by following the same 



strategy as that taken by Gross and Varsted in their now widely cited paper |[41|| . They 
showed that (k, I) moves are ergodic in the space of combinatorially equivalent simplicial 
manifolds by proving that the (k, I) moves are equivalent to the Alexander move. In other 
words, they proved that one can reproduce the effect of an Alexander move by a finite 
series of (k, I) moves. 

In two dimensions, it is easy to show that the Alexander move may be reproduced by 



a (1, 3) move followed by a (2, 2) move [fil]]. In three and four dimensions the equivalence 
of the moves is much less obvious. This is seen only by first proving that 3-spheres and 
4-spheres are locally constructive. A c?-sphere may be transformed into a <i-ball simply by 
removing a single simplex. For the purpose of this explanation, we call these structures 
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punctured ci-spheres. If a single simplex can be grown into a punctured <i-sphere by the 
successive connection of simplices along i faces (where i — 1, . . . , d) whilst preserving its 
topology, then the corresponding ci-sphere is said to be locally constructive f\. In actual 
fact, we are interested in the locally constructivity of S d ~ 2 submanifolds. The relevance 
of this concept will become apparent later in section 



For simplicity, let us first consider local constructivity of S 1 submanifolds of simplicial 3- 
manifolds. We begin with a single link, and then connect new links to either vertex of the 

1- ball (line). Of course, any 1-ball with volume N x > 2 is a punctured 1-sphere, since one 
can connect the two vertex boundaries with a link to form a 1-sphere (circle). Given that 
a link can be grown into a punctured 1-ball without a change in its topology, the 1-sphere 
is locally constructive. In the case of simplicial 4-manifolds, the submanifold in question 
is the 2-sphere. It has been shown that a single triangle can be grown into a punctured 

2- ball (with volume iV 2 > 3) without a change in topology f4lfl . For five dimensional 



triangulations, we are interested in proving that 3-spheres are locally constructive. 

Local constructivity may seem an obvious property of <i-spheres, but to prove this is 
a non-trivial matter for d > 1. Gross and Varsted attempted to prove that the (k, I) and 
Alexander moves are equivalent in arbitrary dimension d, but were unfortunately unable 



to prove the local constructivity of all ci-spheres. In section |2.4| we prove that 3-spheres 
are locally constructive and use the result to show that the Alexander and (k, I) moves 
are equivalent. 



2.4 Equivalence of Alexander and (k,l) Moves 

Consider a closed five dimensional simplicial manifold T which can be expressed as T = 
(ab)P + Q, where P is the complement of the link ab and Q is the remainder submanifold 
of T. An Alexander move transforms T to give 

T — >T A = x(a + b)P + Q, (2.21) 

where a; is a vertex such that x T. The submanifold P is a 3-sphere comprised of, say, 
n tetrahedra. Let ab be a link in T and pqrs be a tetrahedron in P. One can now express 
T as follows. 

T = abpqrs + ab(P — pqrs) + Q (2.22) 
We now perform a (1,6) vertex insertion move on the simplex abpqrs, where x ^ T. 

abpqrs — > [abpqr + abpqs + abprs + abqrs + apqrs + bpqrs)x (2.23) 
5 A more precise mathematical definition of local constructivity is given by Gross and Varsted 
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This transforms T into a new triangulation T\. 

T — > Ti = abx(pqr + pqs + prs + qrs) + x(a + b)pqrs + ab(P — pqrs) + Q (2.24) 
Now T\ is of the form 

Ti = abxRi + x(a + b)Ri + ab(P - R { ) + Q, (2.25) 

where Pi = pqr + pqs + prs + qrs and R\ = pqrs. One could consider R\ as a minimal 
3-ball with S 2 boundary R\. In general, Ri is a 3-ball submanifold of P, with 2-sphere 
boundary Ri, see figure 




n— 1 



Figure 2.4: This illustration represents the two 3-balls P and P — R{ of the 3-sphere 
P, which share a 2-sphere boundary Pj. In the first (last) picture P — Pi (P n _i) is a 
punctured 3-sphere. 

The strategy of the proof involves subtracting tetrahedra from P — Ri which are then 
added to Pj by performing a series of (k, I) moves. After n moves we expect to have 
P n = P and R n = 0. In this instance, Q2.25|) becomes 



T n = x{a + b)P + Q = T A . 



(2.26) 



The resultant triangulation T n is identical to that generated by the Alexander move, given 
by transformation ( |2.21|) . This process proves that the effect of an Alexander move may 
be reproduced by a finite series of (k, I) moves. In other words, it is a proof of their 
equivalence. 



2.4.1 Local Constructivity 

Every tetrahedron that is removed from P — Ri and added to P« must satisfy the three 
conditions. These conditions ensure that the minimal 3-ball Pi may be grown into a 
punctured 3-sphere P„_i without a change in its topology. The conditions insist that 
each tetrahedron of P — Ri added to Pj must: 

• share at least one triangular face with the 2-sphere boundary Pj. 
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not share exactly one triangular face and four vertices with Ri 
leave Tj in the form of equation ( 2.25Q . 



For the moment, let us concentrate on the first two conditions, which relate to the local 
constructivity of triangulations. Clearly, given that T is a manifold, one can always find 
a tetrahedron of P — Ri which shares at least one triangular face with Ri. Therefore, the 
first condition is always satisfied. The importance of the second condition is illustrated 
by figure 12751. One can see that a change in topology occurs if a tetrahedron of P — Ri 
which shares one triangular face and four vertices with Ri is added to Ri. This ensures 
that the topology of Ri remains spherical, at least for i — 1, . . . ,n — 1. 



P-Ri 




Figure 2.5: This illustration shows a tetrahedron in P — Ri which shares one triangular 
face and four vertices with the 2-sphere Ri. 



Proof of Local Constructivity 

One may show that the second condition can always be satisfied by proving that not all 
tetrahedra of P — Ri can share one triangular face and four vertices with Ri. This is 
done by assuming the opposite and showing that the resulting structure is not a manifold 
(reductio ad absurdum). This would prove the local constructivity of 3-spheres. Consider 
a tetrahedron xaba of P — Ri formed by the product of triangle xab and vertex a in Ri. 
These structures are shown in figure |2.6| . 

Given that T is a manifold, one can always identify three triangular neighbours of xab 
in Ri. One of its neighbours, xbc forms a tetrahedron xbcn with vertex k. Clearly, a and 
k are connected by series of links in Ri. The product of these links with xb form a set 
of tetrahedra in P — Ri. One can then repeat this identification for every triangle that 
is common to vertex x. The corresponding tetrahedra are xaba, xbcn, xcdi, xdy6, xyh^j 
and xha/3. Vertices a, /3, 7, 9, 1 and k form a 1-sphere subcomplex of Ri, which is the 
boundary of a 2-ball X. Again, if T is a manifold, then the product of the triangles of 
X with vertex x must form a set of tetrahedra in P — R4. These also have one face and 
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Y 



Figure 2.6: Tetrahedra of P — Ri which share one face and four vertices with Ri. 



four vertices in common with R^. The tetrahedra defined so far form a solid 3-ball handle 
attached to Ri. 

One may then consider the triangles common to vertex y and apply the same process. 
In this instance the product of triangles yhx, yxd, yde, yef, yfg and ygh with vertices 7, 
6, t], (, e and 5 (respectively) form a set of tetrahedra in P — Ri. In the same way, one 
can identify a second 2-ball Y in R iy which shares a series of links 7-^5 with X. The 
product of the triangles of Y with vertex y also form tetrahedra in P — Ri. 

This process may be thought of as a mapping between two separate 2-balls of Ri. 
One can continue to identify tetrahedra in this manner until almost every triangle of 
Ri is accounted for. Eventually, of course, the 2-balls will approach one another. This 
inevitably leads to collapsed tetrahedra and the breakdown of the manifold structure. In 
this case, one is forced to define a tetrahedron with only three vertices, which is impossible. 



This is illustrated in figure |277], where a tetrahedron pqrip is defined by the product of 
triangle pgr and vertex if). This results in a collapsed tetrahedron pqip. 




Figure 2.7: Example of a collapsed tetrahedron. 



These arguments have shown that T is no longer a manifold if one insists that every 
triangle of Ri forms a tetrahedron in P — Ri from the product of itself and another vertex 
of Ri. 
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Additional Remarks 

After this work was completed, the author became aware that the concept of local con- 
structivity was found to have relevance to the problem of enumerating the distinct tri- 



angulations of the 3-sphere pT| , [48]. It was claimed that the existence of an exponential 
bound of the number of distinct 3-spheres may be proven analytically assuming that 
3-spheres are locally constructive. It would appear that our result completes the proof. 

It was pointed out that the local constructivity of 3-spheres implies their algorithmic 
recognisability f|8|| . A manifold M. is recognisable if there exists an algorithm which can 



prove that a second manifold Ai' is homeomorphic to M.. The algorithmic recognisability 
of the 3-sphere has already been proven |4£| |50j . In fact, it has also been proven that 
spheres of dimension greater than four f\ are unrecognisable [|51]. This result must imply 



that spheres of dimension greater than four cannot be shown to be locally constructive. 
One might expect our proof of local constructivity to generalise to higher dimensions. If 
this is the case, then one is forced to conclude that the proof is incomplete (or perhaps 
flawed) . 

In order to complete the proof of equivalence one must show that every tetrahedron added 
to Ri also satisfies the second condition. A tetrahedron pqrs in P — Ri may share one, 
two, three or four triangular faces with Ri. Let us now treat each case individually to 
show that Ti remains in the form of ( |2.25| ). 



One Triangular Face Shared 

Let pqr be a boundary triangle of Ri. pqs, prs and qrs are internal triangles of P — Ri. 
In this case we perform a (2, 5) move on 5-simplices abpqrs and abpqrx as follows. 

abpqr(s + x) — > (abpq + abpr + abqr + apqr + bpqr)sx (2.27) 

This results in the following transformation: 

Ti — > T i+ i = abx(Ri — pqr + pqs + prs + qrs) 

+x(a + b) (Ri + pqrs) + ab(P - Ri - pqrs) + Q (2.28) 
= abxR i+1 + x(a + b)R i+1 + ab(P - R i+l ) +Q (2.29) 

6 It is not known whether 4-spheres are recognisable. 
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Two Triangular Faces Shared 

Let pqr and pqs be boundary triangles of Ri. prs and qrs are internal triangles of P — R4. 
In this case we perform a (3, 4) move on 5-simplices abpqrx, abpqsx and abpqrs as follows. 

abrs(px + qx + pq) — ► {pqr + pqs + ars + brs)abx (2.30) 

This results in the following transformation: 

Tj — > T i+ i = abx(Ri — pqr — pqs + prs + qrs) 

+x(a + b)(Ri+pqrs) + ab(P - Ri - pqrs) + Q (2.31) 
= abxR i+1 + x(a + b)R i+l + ab(P - R i+1 ) + Q (2.32) 

Three Triangular Faces Shared 

Let pgr, pqs and prs be boundary triangles of Ri. qrs is an internal triangle of P — Ri. 
In this case we perform a (4, 3) move on 5-simplices abpqrx, abpqsx, abprsx and abpqrs 
as follows. 

abp(qrx + qsx + rsx + grs) — > {ab + ap + bp)qrsx (2.33) 
This results in the following transformation: 

Ti — > T i+ i = abx(Ri — pqr — pqs — prs + qrs) 

+x(a + b) (Ri + pqrs) + ab{P - Ri - pqrs) + Q (2.34) 
= abxR i+1 + x(a + b)R i+1 + ab(P - R i+1 ) + Q (2.35) 

All Four Triangular Faces Shared 

Now all four triangular faces of pqrs are boundary triangles of Ri. These correspond to 
the last move which transforms T n -\ — ► T n = Ta- It is a (5, 2) move on the 5-simplices 
abpqrx, abpqsx, abprsx, abqrsx and abpqrs as follows. 

ab{pqrx + pqsx + prsx + qrsx + pqrs) — > (a + b)pqrsx (2.36) 

This results in the following transformation: 

Tj — > T i+ i = abx(Ri — pqr — pqs — prs — qrs) 

+x(a + b) (Ri + pqrs) + ab(P - Ri - pqrs) + Q (2.37) 
= x(a + b)P + Q = T A (2.38) 

since P = R n -i + pqrs and R n -\ is the boundary of a punctured 3-sphere given by 

Rn-l = + PQ S + P rs + 1 rs - 
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We have now shown that both conditions of local constructivity are satisfied. This proves 
that the Alexander and (k, I) moves are equivalent in five dimensions. By showing that 
the Alexander and (k,l) moves are equivalent we have also in effect proven that the (k,l) 
moves are ergodic in the space of combinatorially equivalent simplicial manifolds C. □ 

Gross and Varsted were unable to fully generalise their arguments because local construc- 
tivity could not be proved for arbitrary dimension. Nevertheless, it was proved that the 
Alexander and (k, I) moves are equivalent in all dimensions, assuming that spheres are 
locally constructive. Leading from comments made earlier, it seems that our proof of 
local constructivity certainly cannot be generalised to dimensions greater than four. In 
the context of ergodicity, we are interested in the local constructivity of (d — 2)-spheres. 
Therefore, based on Novikov's result (that spheres of dimension greater than four are 
unrecognisable) ]5T| , one concludes that the proof of ergodicity cannot be generalised to 
spheres of dimensions greater than six. 

During the course of this work it was discovered that (k, I) moves had already been 
proved to be ergodic in the space of combinatorially equivalent simplicial manifolds for 
arbitrary dimension j|2], £43] . One cannot extend the proof of ergodicity to the space of five 
dimensional homeomorphic triangulations, due to the negativity of the Hauptvermutung, 
which implies that C $2 H . In other words, there may exist some triangulations which 
cannot be reached by the Alexander or (k, I) moves. The physical significance of this 
possibility is not clear. In fact, the Hauptvermutung is false in four dimensions and re- 
mains unsolved for the 4-sphereQ. The result could have implications for four dimensional 
dynamical triangulations. 



2.5 Non-computability 

Any initial triangulation T may be transformed into any other final triangulation T" 
homeomorphic to T, via a finite series of (k,l) moves E2L B3|. At a practical level, it 
is hoped that the number of moves J\f(k,i) required to transform T into T' is as small as 
possible. By this, we mean a number which modern computers can manipulate - since a 
finite number can be unimaginably large. If Af(k,i) is n °t a computationally managable 
quantity then certain triangulations may effectively never be reached by a Markov chain. 
This would result in an effective breakdown in ergodicity. In this section we consider this 
concept in the context of five dimensional triangulations. 



7 This fact was pointed out to the author by W. B. R. Lickorish [|5^] . 
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Unrecognisability of Manifolds 

A manifold A is said to be unrecognisable if there exists no algorithm which can prove that 
a second manifold B is homeomorphic to A. Unrecognisability is analogously defined for 
simplicial manifolds. This property of manifolds was shown to have relevance to dynamical 
triangulations by Ben-Av and Nabutovsky |53| . They proved that the unrecognisability of 



a simplicial manifold T implies that the number of moves M(k,i) needed to transform one 
triangulation T into a second T' (both of volume < N) is not bounded by any recursive 
function of N. An example of a very fast growing but nevertheless computable (recursive) 
function / is 

f(N) = N N TT • • • TT N N for N E N, (2.39) 

N N 

using Knuth's up-arrow notation |54| where all powers associate to the right 0. 



mjn = m x m x ■ ■ ■ x m x m = m n (2.40) 

n 

m tt n = rn]m]---]m]m (2-41) 

V v ' 

n 

m ttt n — m TT m TT ■ ■ • TT m TT m (2.42) 



Computational Ergodicity 

If A/(fc n is bounded by a computable function then the moves are said to be computation- 
ally ergodic. Note that computational ergodicity of the (k, I) depends on the topology 
of the simplicial manifold in question. The matter of computational ergodicity was first 
addressed in four dimensions by Nabutovsky and Ben-Av |53] because certain 4-manifolds 
are known to be unrecognisable. The number of distinct triangulations (microcanonical 
partition function) W(N^) of a simplicial 4- manifold appears to be bounded by a com- 
putable function of A^ 4 , see chapter |3|. With this in mind, one may infer that N(k,l) is 
non- computable if T can only be transformed into T' (both of volume < A^ 4 ) via inter- 
mediate triangulations with volumes not bounded by computable functions of A 4 . These 
intermediate triangulations may be interpreted as barriers in the space of triangulations. 

Ambj0rn and Jurkiewicz argued that such barriers should exist for all scales of volume. 
They devised an experiment to test for their existence in the triangulation space of S , 
whose recognisability is not known. This involved growing a minimal S 4 triangulation to 
a finite target volume and then shrinking it by increasing k 4 . If the volume eventually 
reduced to iV 4 = 6 then this is taken as evidence supporting the non-existence of barriers 

8 /(l) = 111 = 1, /(2) = 4 TTTT 4 (extremely large). 
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in S 4 . All independent triangulations tested were successfully shrunk to the minimal S 4 . 
Ambj0rn and Jurkiewicz reached the conclusion that either S A is recognisable, or, there 
are very few or no triangulations separated from the minimal S A by large barriers. 

To test the significance of this result, de Bakker [56| replicated the experiment for the 
5-sphere, which is known to be unrecognisable |HJ . Again, no evidence of barriers were 
found. This proved that the results presented by Ambj0rn and Jurkiewicz tell us nothing 
about the recognisability of S 4 . 

Given that S 5 is known to be unrecognisable, one wonders why no barriers were 
observed. Some possible explanation were suggested by de Bakker |5£|, which are also 
applicable to 4-spheres. Perhaps barriers are very large, even for small triangulations. 
Alternatively, the proportion of triangulations separated by large barriers in triangulations 
space may be very small. Consequently, the chance of observing such triangulations would 
also be very small. 

It seems that the whole issue of computational ergodicity is rather academic. Per- 
haps it is more relevant to consider whether Af(k,t) is computationally managable rather 
than computationally ergodic. Computationally managability refers to whether modern 
compters can perform Af(k,i) transformations over a reasonable time scale. It relates to 
the question of whether triangulations may be reached in practice. This is a more appro- 
priate consideration since, as we have seen, even a recursive function can be extremely 
fast growing. How does this issue affect our study of dynamical triangulations in five 
dimensions? We have no evidence so far that barriers exist in the triangulation space of 
S 5 . Even if 5-spheres were recognisable, there still may be triangulations separated by a 
prohibitively large number of moves. 



Chapter 3 
Exponential Bound 



3.1 Canonical Partition Function 

YNAMICAL triangulations is well-defined as a statistical model provided 
that the canonical partition function is an exponentially bounded function 
of the volume in the thermodynamic limit. The purpose of this chapter is 
to prove that this is indeed the case. We begin by reminding the reader 
that the grand canonical partition function of five dimensional dynamical triangula- 
tions is given by 

Z(k 3 , Kb) = exp(/t 3 iV"3 - n 5 N 5 ). (3.1) 
TeT 

If Z includes a summation over topology, then it is automatically divergent. In this thesis, 
the topology is fixed to that of 5-spheres - though the following treatment may be applied 
to simplicial manifolds of any topology. Equation (|3.1|) can be expressed equivalently as 
a sum over all volumes and a sum over all triangulations T of a given volume N$. 

Z(k 3 ,k 5 ) = ^exp(-K 5 A^ 5 ) exp(n 3 N 3 ) (3.2) 

N 5 TeT(N 5 ) 

= J2 e M-K 5 N 5 )Z c (K 3 ,N 5 ) (3.3) 

The canonical partition function Z c (k 3 , N 5 ) is a sum over all triangulations with fixed 
volume N 5 . It is not possible to calculate Z c in practice because the five dimensional (k, I) 
moves all change N 5 . One must therefore allow the volume to fluctuate. Consequently, 
an extra term is added to the action in the grand canonical partition function. 

Z'(k 3 , K 7 ) = E exp(/€ 3 iV3 - k 5 N 5 - 7 (iV 5 - N*) 2 ) (3.4) 
TeT 
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3.1.1 Entropy Function 

The microcanonical partition (entropy) function W(N§) is obviously equal to Z c (k,3, N$) 
when k 3 = 0. Explicitly, we have: 

W(N 5 ) = Yl 1 ( 3 - 5 ) 

T&T(N 5 ) 

= Z c (0,N 5 ) (3.6) 

Notice that in equation (|3.3| ) the exp(—n 5 N 5 ) factors tend to zero exponentially as N 5 — > 
oo. So, in order for Z to be finite, the canonical partition function must be exponentially 
bounded from above function of volume. 

Z c (k 3 ,N 5 ) = exp( K3 iV3) (3.7) 

TeT(N 5 ) 

< exp(K c 5 N 5 ) as N 5 -> oo (3.8) 
Here Kg is the (finite) critical simplex coupling. If we now substitute 

Z c {k 3 , N 5 ) = exp(^iV 5 ) (3.9) 

into ( |3.3|) , we get 

Z = £exp(- K5 iV 5 )exp(4iV 5 ). (3.10) 

It is now becomes clear that Z is divergent if > and convergent if k\ < k$. From 
equations ( |3.7| ) and ( |3.8| ) we notice that K5 is a function of k%. In general, Kg is also a 
function of iV 5 , i.e. 

«g = ^(« 3 ,iV5). (3.11) 

This dependence on N 5 relates to power law subleading corrections to the exponential 
divergence. 

The partition function Z is thus well-defined (finite) only if k\ is finite in the thermo- 
dynamic limit. The question we therefore must address is: does k% have an asymptotic 
limit as N 5 — > 00? The existence of such an asymptote would imply an exponential bound 
on the number of distinct triangulations of a given topology as N$ — > 00. In the next 
section we describe how this can be answered by numerical means. 



3.2 Numerical Approach 

So far, we have shown that an exponential bound exists providing k\ has an asymptotic 
limit as A^5 — * 00. It is hoped that a numerical study will allow us to make a firm 
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statement regarding this matter. Of course, numerical evidence could never constitute a 
proof of existence, or otherwise. 

From equation (|3.9|) we see that the critical simplex coupling /^(fi^, N 5 ) is defined as 

k 5 {k 3 ,N 5 ) = — . (3.12) 

Our Monte Carlo simulations evaluate the following sum. 

Z'(ks,NI,j) = £ exp(/^3-K 5 Ar 5 - 7 (Ar 5 -jV*) 2 ) (3-13) 
Ter 

= ^exp(lnZ c ( K3 ,iV 5 )- ft : 5 iV 5 - 7 (iV 5 -iV*) 2 )) (3.14) 

One can determine the most important contributions to Z'(k3, iV|, 7) by partially differ- 
entiate with respect to N 5 and setting the result to zero. 

= 5>£ - k b - 2 7 (iV 5 - Nl)) exp(ln Z c (« 3 , ^5) - k 5 N b - 7(^5 - N*) 2 )) (3.15) 

Of course, Z'(ks, iV|, 7) is not zero. This leaves us with the following relation. 

Kl = K B + 2rf(N B -Nl) (3.16) 

The saddle point approximation then relates the critical simplex coupling k\ to the mean 
volume (N 5 ). 

k|(« 3 , (iV 5 » = k 5 + 2 7 ((iV 5 ) - Nl) (3.17) 

Here we describe the method of measuring for a given target volume N§. An 
initial arbitrary k b is chosen at the beginning of the Monte Carlo simulation. Once the 
triangulation is fully grown, the simulation then passes into the thermalisation stage, see 
appendix |A.1.3| . Then, k 5 is tuned using the following iterative procedure. 



k 5 — > «' 5 = « 5 + 27 (N 5 - Nl) (3.18) 

Throughout this thesis 7 is fixed to ~ 0.001. The simplex coupling k 5 is updated after 
each tuning iteration. If N 5 > N§, then k 5 is decreased, see equation QA.lj) . This means 
that the (k, I) moves that increase N 5 have a lower change in action AS, whereas moves 
that decrease N 5 have a higher AS. This has the effect of increasing N 5 . The opposite 
is true when N$ < iVj. The algorithm therefore tunes the system to the preset target 
volume. 
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Clearly, for a simplicial 5-manifold we have the inequality N 3 ^$> N$. In fact, one may 
write N 3 = aN 5 , where a is a function of N 5 and a > 1. Using this relation we may 
rewrite the canonical partition function as follows. 

Z c (k 3 ,N 5 ) = exp( K3 iV3) (3.19) 

TeT(N 5 ) 

= exp(aK 3 N 5 ) ]T 1 (3.20) 

TeT(N 5 ) 

Therefore, in order to show that Z c (n 3 ,N 5 ) is exponentially bounded we only need to 
prove this to be true for one value of n 3 > 0. For simplicity we choose to examine the 
divergence of Z c (0, N 5 ) = W(N 5 ), i.e. with k 3 = 0. 

The pseudo-critical simplex coupling k% was calculated using our measurements of 
(iV 5 ), see equation ( |3.17|) . Calculation of the standard error in the mean volume s m ((N 5 )) 
allows us to evaluate the error in k\ using the following relation. 

a(K c 5 ((N 5 ))) = 2 1Sm ((N 5 )) (3.21) 



3.2.1 Exponential Bound for d < 4 

The enumeration of distinct triangulations is a problem that was first addressed in the 
1960s, though purely for mathematical motivations. It was Tutte who first proved that the 
number of distinct triangulations W(N2) of a 2-manifold is exponentially bounded [57 |. 
In fact, it has been shown that for a 2-manifold of genus g we have 

W(jg,N 2 ) = exp(/i c AT 2 )A^)- 3 (l + O^" 1 )), (3.22) 

where [i c is the critical constant and ■y(g) = \g — \ p8| . 

The generalisation of dynamical triangulations to dimensions greater than two natu- 
rally led to the question of whether such an exponential bound exists in these dimensions. 
This matter was first investigated numerically by Ambj0rn and Varsted, who studied S 3 
triangulations of volumes up to A^3 = 14k |58| . It was found that their results were consis- 
tent with the existence of a bound. A second study using large volumes (up to ^3 = 128k) 
corroborated their findings ||59| . Attempts were subsequently made to construct an ana- 
lytic proof - with varying degrees of success 



With no known formal proof in four dimensions, the exponential bound issue was also 
first considered numerically. Results produced by Catterall et al. (for volumes up to 
N4 = 32k) seemed to indicate that ^(A^) was not exponentially bounded This 
scenario corresponds to a factorial divergence of entropy 

Z c ( Ko ,iV4)~(iV 4 !) 5 , (3.23) 
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where 5 is a positive constant[|, estimated to be 5 = 0.026(5) |6"I|| . Further numerical work 
was presented by Ambj0rn and Jurkiewicz (with iV 4 < 64k), which contradicted the results 
of Catterall et al. ||62|| . The matter was resolved when a third study (with N4 < 128k) 
produced evidence consistent with the existence of an exponential bound |j64"| . More 
recently, certain advances have since been made towards an analytic proof |12|, |65|, |57| . 
We now turn our attention to five dimensions. History has shown that analytic proofs are 
difficult to construct; so let us attack the problem by numerical means. 

3.3 Simulation Results 

One may establish the possible existence of an asymptotic limit, and hence an exponen- 
tial bound, by first measuring k§((N 5 )) over a range of finite volumes N$. The functional 
form of Kg((iV 5 )) may then be determined by fitting curves to the data using chi-squared 
hypothesis testing - see appendix |B.4| for details. Similar studies in four dimensions have 
shown that this approach can be very costly in terms of computational resources. Fur- 
thermore, the results require careful interpretation before any conclusion may be drawn. 
These facts influenced the planning of our strategy - mainly in terms of the chosen target 
volumes. 

The most effective approach clearly entails striking a balance between the range 1Z 
over which volumes are sampled, the number M of volumes we sample and the accuracy^ 
A with which we measure for a given volume. It is reasonable to suppose that for 
a given amount of computational resources there exist optimal values of these variables 
that give the most reliable determination of ^((iVs)). Calculating these optimal values 
is however rather tricky. 

It was first decided that M should be kept relatively small, thus avoiding the problems 
encountered by de Bakker |6~3fl . Obviously, to have any hope of achieving our goal, we 
must have M > 3. The accuracy A typically increases as where M. is the number of 
independent measurements taken. Consequently, gains in accuracy will tend to diminish 
for longer Monte Carlo simulations. Our main concern is that A and 1Z are chosen such 
that logarithmic divergence may be easily distinguished from power law convergence, see 
later. In other words, the accuracy should be sufficiently great that the form of kI({N 5 )) 
is compatible with either divergence or convergence, not both. 

After taking these issues into consideration, the number of volumes M was chosen to 
be five, ranging over an order of magnitude. Measurements of N$ were taken over one 
million Monte Carlo sweeps for each target volume: 5k, 10k, 20k, 40k and 70k. The 

1 exp(xS\nx) — x xS ~ (x\) s 

2 In this sense 'accuracy' is loosely interpreted as the inverse of the error in «§. 
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N 5 


5k 


10k 


20k 


40k 


70k 


(N 5 ) 


5000.1(2) 


10005.0(6) 


20020.0(1.0) 


40002.3(2) 


70011.8(1.6) 


ln(JV 5 ) 


8.5172 


9.2108 


9.9045 


10.5967 


11.1564 




0.8117(5) 


0.8277(6) 


0.8404(6) 


0.8498(4) 


0.8571(6) 



Table 3.1: Pseudo-critical simplex coupling kI({N 5 )) for target volumes 7V|. 



3.3.1 Curve Fitting Analysis 

One may now fit curves to the data and hence determine whether or not /^((A^)) is finite 
as A^5 — > oo. We begin by fitting a power law curve of the form 

K c 5 = a-b(N 5 )- c , (3.24) 

which has an asymptote k% — ► a as A^ — > oo, where a, b and c are positive. This curve is 
therefore consistent with the existence of an exponential bound. Conversely, a logarithmic 
curve of the form 

K c 5 =p + gln(iV 5 ) (3.25) 

does not have an asymptotic limit since k§ — > oo as N5 — > 00, where q > 0. In other 
words, such a curve is incompatible with an exponential bound. Using the chi-squared 
method, the best fit curves were respectively found to be 

K c 5 = 0.890(7) - 1.3(4) (AT 5 )- a33(4) (3.26) 
X 2 ~ 1.38 with 2 d.o.f. 



and 



k% = 0.667(4) + 0.0172000(1) ln(iV 5 ) (3.27) 
X 2 ~ 61.8 with 3 d.o.f., 



• 1 



where the errors in the coefficients are shown in parentheses. These curves are plotted 



along with the numerical data of table |3J] in figure |3.1|P|. 



From figure 3.1 it is obvious that the power law fit is much better than the logarithmic 



fit. The power law curve passes almost exactly through all five data points. This may not 



3 The coefficients were not rounded off, in order to show the very best fit. 
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0.87 



0.86 - 




mean volume 

Figure 3.1: Pseudo-critical simplex coupling k|(x) where x = (AT 5 ), with best fit power 
law (solid) and logarithmic (dashed) curves. (In the key 'log' actually means In.) 

be surprising given that the power law has fewer degrees-of-freedom. Nevertheless, the 
observed goodness-of-fit is proven by comparing their respective values of x 2 P er degree- 
of-freedom. The power law fit has \ 2 ~ 0.7 per d.o.f., compared with y 2 w 21 per d.o.f. 
for the logarithmic fit. (The errors in the measurements of (Ns) are very small, and are 
therefore neglected in the curve fitting analysis.) 

The goodness-of-fit of these curves can be made more apparent by plotting kI({N 5 )) 
against ln(A^). On such a graph a logarithmic divergence would appear as a straight line. 
This is done by making the following transformations. 

x = hi(N 5 ) (3.28) 
K c 5 = a - b (N 5 )~ c — ► k\ = a-6exp(-cx) (3.29) 
K c 5 =p + q\n(N 5 ) — > Kl=p + gx (3.30) 



Curves of the form ( 3.29| ) and (|3.30 ) were fitted to the numerical data. Not surprisingly, 
the coefficients a, b, c, p and q of (|3.29| ) and ( |3.30| ) were the same as those measured earlier 
for equations ( |3.26|) and (|3.27|) . This exercise served as an useful check of the numerics. 
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The best fit curves were found to be 



K c 5 = 0.890(7) - 1.3(3) exp(-0.33(4)x) 



(3.31) 



X 



1.40 with 2 d.o.f. 



and 



K% = 0.662(2) + 0.0173(2)x 



(3.32) 



X 



61.6 with 3 d.o.f.. 



These curves are plotted along with the data in figure 3.2. The best fit curves reinforce 



0.87 




8 8.5 9 9.5 10 10.5 11 11.5 



x 

Figure 3.2: Pseudo-critical simplex coupling «§(x) where x = ln(A^), with best fit power 
law (solid) and logarithmic (dashed) curves. 



the result found earlier that the power law curves fit the data much better than the 
logarithmic curves. This body of evidence appears to be consistent with the existence of 
an exponential bound. 
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Three Parameter Logarithmic Fit 

The logarithmic fit may be improved by introducing a third parameter in the form of a 
scale term r to give 

K c 5 =p + q\n((N 5 )-r). (3.33) 
This functional form results in a much better fit of 

k\ = 0.72(3) + 0.013(2) ln((A^ 5 ) - 3139(12)) (3.34) 
X 2 ~ 2.3 with 2 d.o.L 



This new logarithmic fit is clearly much better than the previous fit of (|3.27|) . From 



figure |3]3|(a), it is clear that fits ( |3.26| ) and ( |3.34| ) are equally good. Again, these curves 



are more easily compared when plotted against a log scale. This is done by making the 
following transformations 

x = ln((JV 5 )-r) (3.35) 
K \ = a - b(N 5 y c — ► k c 5 = a-b(exp(x) + r)- c (3.36) 
k c 5 = p + qln((N 5 ) - r) — ► 4=p + gx (3.37) 

Figure |3"lj|(b) also shows that both curves are very good fits to the data. 

In conclusion, both fits are compatible with the data. Therefore one cannot use these 
results to make a statement regarding the possible existence of an exponential bound 
in five dimensions. From figure |3]3|(b) it is clear that the two curves begin to diverge 
steeply for x > 11.5. This illustrates the fact that as TZ is increased, one fit will become 
progressively better than the other. It will become increasingly less likely that both curves 
give equally good fits. 

Smaller Triangulations 

Unfortunately, it was not really possible to measure for volumes much greater than 
70k; the limiting factor being computer power 0. Further measurements were made of 
for N{: < 5k, in order to check their compatibility with ( |3.26j ), ( |3.25| ) and ( |3.34| ). Again, 



iVs was measured over one million sweeps for target volumes lk, 2k, 3k and 4k. The 
results are presented in figure |3T4". 



It is plain to see that the new data is not even remotely consistent with any of the 
best fit curves. By differentiating (|3.24|) , fl3.25| ) and fl3.33|) twice with respect to volume 



'For example, one million sweeps for iVj = 128k takes at least 3 months, see appendix 
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0.87 



0.86 - 




0.8 1 1 1 1 1 1 1 1 1 1 

7 7.5 8 8.5 9 9.5 10 10.5 11 11.5 

(b) 

Figure 3.3: (a) Pseudo-critical simplex coupling ^(x) where x = (N 5 ), with best fit power 
law (solid) and logarithmic (dashed) curves, (b) Pseudo-critical simplex coupling k^(x) 
where x = ln((JV 5 ) — r), with best fit power law (solid) and logarithmic (dashed) curves. 
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00 

a 
a, 
§ 

CJ 

x 

OJ 



13 

CJ 

■a 

CJ 
I 

o 



0.87 
0.86 
0.85 
0.84 
0.83 
0.82 
0.81 
0.8 
0.79 
0.78 




0.88956- 1.2568*exp(-0.32668*x) 
0.66638+0.017295*x 
0.71543+0.012774*log(exp(x)-3138.5) 

data 



10 



11 



12 



13 



Figure 3.4: Pseudo-critical simplex coupling k§(x) where x = ln(A^), with best fit power 
law (solid), logarithmic (dashed) and scale term logarithmic (dotted) curves. 



we (respectively) get: 



dV 5 

diV 5 2 

d 2 «g 

d 2 «g 
diV? 



-bc(c + l)A^ (c+2) 



-qN 5 ~ 2 

-q (N 5 -r)- 2 



(3.38) 
(3.39) 
(3.40) 



Obviously, these functions are negative for all N 5 . If our numerical data is plotted against 
N$ we see that the second derivative of k$ is positive near N 5 w 3k (x « 8). This means 
that there exist no functions of the form ( |3.24|) , ( |3.25|) or (|3.33|) that can accurately fit 
all the dataQ. From this we can conclude that the data for iV| < 5k is not described by 
the large N 5 limit of nl(N 5 ). In other words, nl(N 5 ) must have higher order corrections, 
thus giving either: 



a-b(N 5 )- c + p(N 5 ) q + ■■ 



(3.41) 



5 By 'accurately' we mean: pass through each error bar. 
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or 



s + t\n(N 5 )+p(N 5 ) q + 



(3.42) 



This shows that we must discard all data for N 5 < 5k from any curve fitting analysis. It 



is interesting that very similar effects were observed in four dimensions [62 



3.4 Concluding Remarks 

In conclusion, it is not possible to determine whether kI(Ns) is finite as N$ — > oo. Hence, 
one cannot make any firm statement regarding the existence of an exponential bound. 
The data is entirely consistent with convergent as well as divergent curves - the power 
law fit is marginally better. The issue may only be resolved by measuring k§ for larger 
triangulations, as was the case in four dimensions p3|. As we increase the range TZ over 
which measurements are taken, it will become easier to differentiate between convergent 
and divergent data. 

This whole approach raises an important question: how can we ever be sure that 
our data corresponds to the true large N5 behaviour? In general, we expect subleading 
corrections to contribute significantly to k^(N 5 ) for small volumes - but what are 'small' 
volumes? This problem may be tackled (in principle) in the following way: 

• First measure k c 5 over large TZ (orders of magnitude). 

• Divide the data into a small number of sets, say five. 

• Then determine the best fit coefficients for each set of data. 

• The final step involves comparing the sets of coefficients. 

If the coefficients are the same (within errors) over a series of consecutive sets, then the 
fit is likely to represent the true large N$ behaviour. Of course, the measurement of k\ 
over orders of magnitude is no easy task in practice. 



Chapter 4 
Phase Structure 



4.1 Action 



=^|j HE action normally used in dynamical triangulations is the discre- 
tised Einstein-Hilbert action S[T,K d -2\- Like its continuum counterpart 
Slg^u, A, G], it has two terms derived from the physical notions of spacetime 
curvature and the cosmological constant A, which are linear in and N d 



respectively. For the following analysis, it is convenient to group the d + 1 elements of 
the /-vector into a set Md- 

M d = {N ,...,N d } (4.1) 

In this thesis we do not restrict ourselves to the study of dynamical triangulations with 
the conventional Einstein-Hilbert action, which is derived for arbitrary dimension in sec- 
tion [TXT1 . 

S[T, K d _ 2 ] = K d N d - K d _ 2 N d _ 2 (4.2) 

Instead, we consider the most general action S d that is linear in the members of Md for a 
given dimension d. Let S d be a linear function of the members of a set where 



Mi c Af d . 

One may then define a second subset M' d of Af d such that 



MuJVicjV, 



(4.3) 



(4.4) 



If the knowledge of M d is sufficient to determine the members of M' d using the Euler- 
Poincare and Dehn-Sommerville relations, then S d is the most general action. Further- 
more, the definition of S d stipulates that the number of members of M d must be minimal 
and that N d e A/"|. 
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In this case, the discretised Einstein-Hilbert action corresponds to a certain limit of 
the most general action. Each term in the action has an associated coupling constant. 
Therefore by taking this approach we, in effect, consider a larger space of parameters. 
However, as we shall see, this is not always the case. The question we must ask ourselves 
is: for a given dimension d, knowledge of which A/"J allows us to calculate the members of 
M' d l This question is now answered for dimensions two to five inclusively []. 

For d = 2, there are two linearly independent Euler-Poincare and Dehn-Sommerville 
relations, (|4.5| ) and (|4.6|) , between three variables N , Aq and N 2 . This means that the 
most general action has only one term - which is linear in N 2 . 

N -Nx+N 2 = 2 (4.5) 
2Aq - 3iV 2 = (4.6) 

For d = 3, the four variables Nq, Aq, N 2 and A^ are related by two linearly independent 
Euler-Poincare and Dehn-Sommerville relations: (4.7) and (4.8). 



N - Aq + N 2 - N 3 = (4.7) 
N 2 -2N 3 = (4.8) 

In this case, the most general action Sf has two terms. It is important to note that Sf 
cannot be a function of any two terms. For example, knowledge of N 2 and A3 does not 
allow us to calculate A" and Aq. 

For d = 4, there are five variables N , Aq, N 2 , N 3 and N± related by three linearly 
independent Euler-Poincare and Dehn-Sommerville relations: Q4.9|) , ( |4.10|) and (|4.11| ). 
Again, this means that the most general action Sf has two terms. Similarly, knowledge 
of A*3 and does not allow us to solve for the remaining variables. 

^0-^1 + ^2-^3 + ^4 = 2 (4.9) 
2Ni - 3N 2 + 4N 3 -5N 4 = (4.10) 
2A3-5A4 = (4.11) 

It is a mere coincidence that the number of terms in Sf and Sf is the same as in the 
Einstein-Hilbert action S[T, Kd-2]- It is important to remember that there is nothing 
physically significant about this fact. 

Usually, in practice, the action used in the study of dynamical triangulations is 
S[T, Kd-2], defined by equation (|4.2|). However, in certain cases a different, but equiv- 
alent, action S[T, kq] is utilised in which the 'curvature term' is linear in A^o- 

S[T, k ] = K d N d - k N (4.12) 



lr This is done formally using Gaussian elimination. Basically, the number of terms of the action is 
d + 1 — 7i, where n is the number of independent equations relating the d + 1 variables. 
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This is done purely for practical reasons; Monte Carlo simulation computer programs may 
be written for arbitrary dimension very efficiently in this way |68| (see appendix |A.f| ). As 
we shall see, this is not generally possible for d > 4. 

Five Dimensions 

The situation becomes a little more complicated in d = 5 as we have six variables N , 
N±, iV~2, A3, A4 and A5 constrained by three linearly independent Euler-Poincare and 
Dehn-Sommerville relations: ( 4.13|) , ( [4.14D and (|4.15| ). 



^0-^1 + ^3-3^5 = (4.13) 
A 2 - 2A 3 + 5iV 5 = (4.14) 
A 4 -3A 5 = o (4.15) 

This means that the most general action S| has three terms 0. We are now faced with a 
situation where Sf has more terms than S[T, /c 3 ]. It is important to note that S[T, k 3 ] 
cannot be rewritten in the form of ( [4.12|) , because it is impossible to express A 3 in terms 
of Ao and A5 alone. Again, the most general action cannot be expressed in terms of any 
set of three variables. Using equations (|4.13| ), (|4.14|) and ([4.15|) we find that knowledge of 



Ml = {A , A 3 , A 5 } allows us to compute the values of Af' d = {Ni, A 2 , A 4 }. We therefore 
choose our most general action to be 

S 9 5 [T, Ko, k 3 ] = k 5 N 5 - k 3 N 3 - k N . (4.16) 

Minor modifications were made to the Metropolis algorithm in order to accommodate this 



new action; details of which are given in appendix [A.l| . Of course, an equally acceptable 
action would be 

S![T, Kl , 4} = - ^N 3 - Kl Nx. (4.17) 

One can relate actions (|4.16|) and ( |4.17|) using equation (|4.13|) , which gives the following 
transformations. 

4 = /€5-3k (4.18) 

K3 = K 3 — K (4.19) 

ki = k (4.20) 

Clearly, the five dimensional Einstein-Hilbert action S[T,k 3 ] corresponds to the kq = 
limit of Sf [T, k,q , k 3 ] . On that note, let us now clarify our reasons for considering these 
general actions. 

2 This fact was first pointed out by de Bakker [M. 
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The principle reason for researching five dimensional dynamical triangulations is based on 
the realisation that the most general space of parameters is three dimensional. Rather than 
attaching any great significance to the discretised Einstein-Hilbert action, we consider the 
most general action that is linear in the components of the /-vector. In doing so, we are 
treating dynamical triangulations as a statistical model of geometry in its own right 
without worrying too much about 'gravity'. 

It is suspected that this larger space of couplings may disclose rich phase structure, 
which could conceivably lead to interesting physics. It is hoped, at the very least, that 
this perspective may uncover new facets of dynamical triangulations, which could teach 
us more about the theory in general and its potential application to quantum gravity. For 
example, this viewpoint may shed new light on the phenomenon of singular vertices (see 
chapter §). These are known to be an important feature of dynamical triangulations, and 
may even drive the phase transition in four dimensions. 



4.1.1 Space of Parameters 



So far, we have established that S§ has three terms, and hence three parameters. Only two 
of these parameters are freely independent; the simplex coupling k§ is effectively fixed by 
the chosen target volume (see chapter ||). This means that the space of coupling constants 
is a two dimensional surface C 2 defined at K5 = ^5(^0, K 3, N5), as shown schematically in 



figure fO| . Points on C 2 are identified by the coordinates (k , k 3 ). 
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Figure 4.1: Two dimensional coupling constant surface C 2 . 



Our aim is to search for non-trivial phase structure in C 2 . As C 2 is two dimensional, 
a thorough exploration is likely to be very costly in terms of computer power and time. 
For this reason, all numerical work requires careful forward planning to make the best 
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possible use of available resources. We shall begin our investigation by considering the 
k 3 = and k,q = limits of Sf [T, k , « 3 ]. It is hoped that these preliminary simulations 
will allow us to identify phase transitions and characterise the associated phases. The 
results are presented and analysed in section fO. 



If any phase transitions are found, then the next natural step is to determine their 
order, either by finite size scaling analysis or other means (see section |4.3| ). Our main con- 
cern is to determine whether a transition is continuous or otherwise, as we are ultimately 
interested in the existence of physical continuum limits. The final stage of the numerical 
work involves exploring regions of C2 where k 7^ and k 3 7^ (see section |4.4| ). These 
simulations will be limited to strategically chosen regions of C2. Before going any fur- 
ther, let us first review some relevant background theory of phase transitions and critical 
phenomena. 



4.1.2 Critical Phenomena 

Dynamical triangulations can only give us a non-perturbative formulation of quantum 
gravity if it has a physical continuum limit. Such a limit is reached by reducing the 
lattice spacing (link length) a to zero. This can only be done if the system is at a critical 
point. 

One can define the connected correlation function between two lattice spin variables 
<7j and o~j (at lattice sites % and j respectively) as 

?ij = (°i°j)-(°i)(oj), (4-21) 

where (<7j) represents the expectation value of a at site i. In general, the correlation 
function is a function of x, the distance separating sites i and j in lattice units. It is a 
measure of the fluctuations in a across the lattice. Away from any critical point we find 
that correlations between lattice variables fall off exponentially for large x as 

Fijix) ~ar A exp^-|^ , (4.22) 

where A is a positive number and £ is the correlation length (in lattice units). At high 
temperatures <Tj become uncorrelated as x — > 00 since the mean value of o~i is zero. 
Typically, at low temperatures we expect the lattice variables to be perfectly correlated. 
However, the fluctuations in <7j are zero. 

A continuum limit exists only if the ensemble is invariant under conformal transfor- 
mations. This occurs if there are fluctuations in the spin variables at all length scales 
(up to the size of the lattice). For this to be the case £ must diverge. Correlation lengths 



CHAPTER 4. PHASE STRUCTURE 



53 



only diverge at continuous phase transitions. At critical points the microscopic details of 
the lattice become irrelevant and T^x) falls off according to a power law function 

ryx)~x- A ', (4.23) 

where A' = d — 2 + rj and d is the dimensionality of the system and rj is a critical exponent. 
In this instance, the renormalised mass m r is related to £ by 

m r = (^a)- 1 . (4.24) 

Therefore, in order for m r to be finite we require £ —>■ oo as a — > 0. 



4.2 Preliminary Monte Carlo Simulations 

Our primary objective is to map the phase structure of Ci- This is done by measuring 
quantities that are sensitive to phase transitions, i.e. derivatives of the free energy F. 

Foe log Z (4.25) 

Past studies have shown that it can also be useful to measure geometrical quantities such 
as geodesic distances. Monte Carlo simulations are run over a wide range of Kq and K3 
for triangulations with target volumes of = 10k. Measurements of N Q , N 3 , N 5 and 
the average geodesic distance g (see section |4.2.2j ) were taken at regular intervals over 10 5 



sweeps. This range was progressively widened until kinematic upper and lower bounds 
were reached. Knowledge of these bounds aid the characterisation of phases and hence 
leads to a more complete picture of C%. Finally, the resolution is increased by running 
more simulations closer to the presumed phase transitions, thus identifying their location 
with ever greater precision. 



4.2.1 Derivatives of the Free Energy 

Let us begin by considering first derivatives of the free energy. Along the Kq axis of C2, the 
first derivative of the free energy with respect to k is the first cumulant of the distribution 
of No. It is proportional to the vertex density po- A discontinuity in p as a function of 
k,q is indicative of a first order phase transition. 

Po = I (4.26) 

Figure |4.2| (a) shows two plots of the mean vertex density {po), in which the solid (dashed) 
curve shows the relationship between (p ) and k (k 3 ) respectively. (Note that the curves 
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Figure 4.2: (a) Mean vertex density (p ) versus the vertex coupling kq (solid) and the 
tetrahedron k 3 (dashed), (b) Mean scalar curvature (R) versus n 3 (solid) and k (dashed). 
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intersect at kq = k 3 = 0, as one would expect.) From this graph, it is immediately 
apparent that both plots highlight two distinct phases. Basically, in each case we find a 
weak coupling phase with large p an d a strong coupling phase with small p . 

One can plainly see that both curves have the same asymptote of po ~ 0.200 for n — > 
00 and k 3 — > 00, which is marked by a dotted line. It appears that this limit represents 
the kinematic upper bound of po- This observation clearly opens up the possibility that 
the K3 = and k,q = limits share a common weak coupling phase. 

In contrast, the strong coupling phases are distinct. The kq — > —00 limit has an 
asymptote of po ~ 0.010, which one might assume to represent the kinematic lower 
bound of po- In comparison, there does not appear to be an obvious asymptotic limit for 
K3 — > —00. Indeed, po fluctuates significantly as a function of K3, but only for k 5 < 3% for 
K3 < — 5. From these results it is evident that po ~ 0.010 is not the true kinematic lower 
bound. This fact is somewhat surprising. A dotted line indicates a possible lower bound 
at po ~ 0.005. Time series plots of Nq have shown that all ensembles are thermalised. 
It seems likely that the ensembles for n 5 < 3k 3 for k 3 < — 5 have become trapped in 
metastable states, and that an asymptotic limit does exist for p ~ 0.005. 

One physically tangible observable is the scalar curvature R, which is effectively the mean 
hinge local volume. It is proportional to the first derivative of the free energy along the 
K3 axis. 

(4.27) 

Figure pT2| (b) shows two plots of the mean scalar curvature (R ), in which the solid (dashed) 
curve shows the relationship between (R) and k 3 (kq) respectively. By comparison with 
figure ^2|(a), we find that both pairs of plots are qualitatively very similar. This illustrates 
the intimate (though indirect) connections between certain members of Md- The mean 
scalar curvature tends to an asymptotic value of R ~ 4.00 in the weak coupling limit, 
which appears to represent the kinematic upper bound of R. One also finds an asymptote 
of R ~ 3.35 for k — > —00. A dotted line highlights a possible kinematic lower bound of 
R w 3.11. Closer inspection reveals a small kink near k 3 pa —1. The possible significance 
of this will be investigated later. 

So far, we have already detected non-trivial phase structure in C 2 . Furthermore, we 
have uncovered interesting similarities and differences between phases along the kq and 
k 3 axes of C 2 - These are now investigated further by looking at second derivatives of the 
free energy. 



Phase transitions are normally identified by divergences in normalised second cumulants, 
which are essentially second derivatives of the free energy. Let us first consider the k 3 = 
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limit. The normalised second cumulant of 02(^0,^5) is a measure of the fluctuations in 
N and may be interpreted as the vertex specific heat. 

N S ) - m Qkq2 m (4.28) 

Figure fO](a) is a plot of c 2 (k ,N 5 ) versus k . The sharp peak at n = 1.5 indicates the 
presence of a phase transition in this vicinity. Analogously, a divergence in the tetrahedron 
specific heat 02(^3,^5) signifies the existence of a phase transition along the K3 axis. 

^ (N S ) dK 3 2 (iV 5 ) * ' 

Indeed, a peak is observed at k 3 = 0.45 (see figure P~3|(b)). It is interesting to see a small 
peak in 02(^3,^5) near k 3 = —1, which coincides with the observed kink in (R(k 3 )) (see 
figure ^72|(b)). This may indicate the presence of a second phase transition. From visual 
inspection, the peaks in 02(^0,^5) and 02(^3,^5) appear to coincide with the points of 
inflexion of (po(^o)) and (R(k3)), as one would expect. Collectively, these results prove 
that five dimensional dynamical triangulations has a non-trivial phase structure. This 
was not entirely unexpected given the situation for d < 4. 



4.2.2 Geodesic Distances 

Triangulated lattices may be viewed as discretisations of continuum manifolds. We have 
already established that the concept of curvature may be extended to triangulations. In 
fact, geodesies also have a natural definition on such discretised spaces. 

A geodesic may be defined as the shortest path connecting two points on a manifold. 
For our purposes, this definition requires some modification since we deal with simplicial 
manifolds. We understand a geodesic between two arbitrary simplices af and af, to mean 
the minimum number of simplices that link af and af along faces. (Alternatively, one 
could define a geodesic as the shortest distance between two vertices along links. Both 
definitions would give the same qualitative results.) Figure [14| illustrates the geodesic 
(dotted line) between simplices a\ and o\ of a two dimensional simplicial manifold^. In 
this particular case the geodesic distance is 7 units. 

By measuring the average geodesic distance ~g one can learn more about the geometric 
of aspects dynamical triangulations. For a given Monte Carlo configuration, one can 
evaluate g by calculating the geodesic distances between many random pairs of simplices. 
Ideally, one ought to average over all pairs. However, this approach can become very costly 

3 Unlike geodesies defined on continuum manifolds, 'discrete' geodesies are not necessarily unique. In 
other words, there may exist more than one path of the same length connecting two simplices (or vertices). 
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Figure 4.3: (a) Vertex specific heat C2(k,q, N 5 ) versus the vertex coupling kq. (b) 
dron specific heat 02(^3,^5) versus the tetrahedron coupling k 3 . 
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Figure 4.4: Geodesic between two 2-simplices of and erf. 

in terms of computer resources since the number of pairs grows factorially with volume. 
Fortunately, reasonable estimates of ~g may be had by sampling just a small proportion of 
the total number of pairs. This is done by randomly selecting a predetermined number 
of 'origin simplices'. Then, for each origin simplex the algorithm calculates the geodesic 
distance to all other simplices, thus giving an average value g. 

Figures [4751(a) and [151(b) show plots of the mean (value of the average) geodesic distance 
(g~) versus kq and K3 respectively. In both cases we find that the weak coupling triangula- 
tions are elongated, with average geodesic distances of about 45 units. This is consistent 
with the existence of a common weak coupling phase. In figure |4.5| (a), the k,q — ► — 00 
asymptote corresponds to triangulations with g ~ 9. For k 3 — > —00, we observe fluctua- 
tions commensurate with those in po and R. In light of these findings, we shall adopt the 
nomenclature prevalent in four dimensions, in which the strong and weak coupling phases 
are commonly referred to as crumpled and elongated phases respectively. 

So far, we have discovered that the mean geodesic distance is as much as four times 
greater in the elongated phase than in the crumpled phases. When this fact is coupled 
with the large fluctuations in g~, we are naturally led to suspect that the dominant tri- 
angulations of the elongated phase may be branched polymers. These tree-like structures 
(with variable branching order) have been observed in two, three and four dimensional 
dynamical triangulations. This matter will be resolved later in chapter [|. 

Previous studies have shown that the fluctuations in g are highly sensitive to phase tran- 



sitions |7T). The normalised second cumulant c 2 (g,N 5 ) of the distribution of g is given 
by 

C2&N 5 )= {92) { -^ 9)2 . (4.30) 
It was found that 02(17, N5) is sensitive to both phase transitions. In the crumpled phases 
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Figure 4.5: Mean geodesic distance (g) versus (a) the vertex coupling k and (b) the 
tetrahedron coupling k 3 . 
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C2(g, N§) is very small, typically of order 10 -6 . Near the transitions we find that C2Q?, N5) 
increases very sharply by orders of magnitude and fluctuates greatly in the elongated 
phase. Hence, c 2 {g,N 5 ) behaves as an order parameter. 

The large fluctuations in g are most probably due to long autocorrelation times (com- 
pared to the measurement intervals). The branched nature of these triangulations means 
that many (k, I) moves must be performed to generate statistically independent configu- 
rations. The local update algorithm can be made more efficient by utilising a technique 
known as baby universe surgery [|69j |70l - In general, this involves identifying a branch B of 
a triangulation T and cutting it off along an S^" 1 submanifold of volume Nd-i- Another 
portion of the manifold V homeomorphic to B d is selected at random. If the S 1 bound- 
ary of V also has volume N^-i, then V and B may be exchanged and reconnected. This 
procedure is a non-local move and generally decreases the autocorrelation time. These 
non-local moves were not employed in our simulations. 



4.3 Finite Size Scaling 

So far we have discovered clear signals indicating the existence of two phase transitions 
along the kq and K3 axes of Ci- Our next goal is to ascertain the order of these phase 
transitions. More specifically, we hope to determine whether the transitions are of first 
order or continuous, using finite size scaling methods. This is achieved in practice by 
looking at how second cumulants scale over a range of target volumes N\. 

As we shall see, this type of numerical work requires substantial amounts of com- 
putational resources, which obviously increase with volume. For this reason, the target 
volumes to be studied were chosen carefully, with a number of factors taken into con- 
sideration. These include the available CPU (central processing unit) time, computer 
performance^ and possible finite size effects. It was first decided that only three target 
volumes would be studied, which is the bare minimum. Past studies of three and four di- 
mensional dynamical triangulations were also influencing factors, for they served as useful 
reference points. The target volumes were finally chosen to be 10k, 20k and 30k. 



In general, the critical couplings are functions of volume, (N5) and K3 = 

k^(N 5 ). Presumably, a prerequisite condition for well-defined thermodynamic limits is 
that Kq(N 5 ) and k^(Ns) converge as N 5 — > 00 |jfl| . Unfortunately, given that we have 



only three data sets, it is not possible to determine whether this is the case. As an 
example, consider the function Kq(Ns). If Kq(Ns) is a monotonic function, then one can 



4 



Computer performance tests show that Monte Carlo simulation times grow faster than linearly with 



volume. Details are given in appendix A. 2 
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always find perfect fits for the two following functional forms (where a, b and c are positive 
constants). 

divergent: /eg = a + bN 5 c (4.31) 
convergent: Kq = a - bN 5 ~ c (4.32) 

In other words, it is impossible to differentiate between the divergent and convergent 
functions. At this stage, one can only assume that Kq(oo) and k^(oo) are finite. To 
produce evidence that supports or contradicts this assumption would require an extensive 
numerical effort, similar in scale to that presented in chapter |3|. It is therefore not in our 
interests to pursue this line of research at the present time. 

4.3.1 Second Cumulant Scaling 

It is assumed that the normalised second cumulant 02(^05-^5) scales with volume at the 
phase transition as 

c 2 (K ,N 5 )~N 5 a /» d ", (4.33) 

where a and v are critical indices and dn is the internal Hausdorff (fractal) dimension Q. 
These quantities are related by Josephson's scaling relation. 

a = 2-vd H (4.34) 

The product vdn has the bound vdu > 1. A first order phase transition is characterised by 
having vdu = 1- Hence, using Josephson's relation we see that in this instance c 2 (k , N 5 ) 
would diverge linearly with volume. The phase transition is continuous if < a < 1 and 
1 < vdu < 2. Relation ( |4.33|) may be rewritten simply as 

c 2 (ko,A 5 ) ~ A 5 A , (4.35) 

where A = a/ vdn- A continuous phase transition therefore has < A < 1. Our task 
is to measure the normalised second cumulants at the phase transitions for each target 
volume and hence determine the scaling behaviour. 

The accuracy with which one can measure second cumulants at phase transitions de- 
pends on two factors. Apart from the usual statistical error, an added error is incurred 
from the uncertainty in the location of the phase transition itself. This drawback is exacer- 
bated if the location is determined from measurements of second cumulants. Fortunately, 
measurements of C2^g, N5), which behaves as an order parameter, allow us to locate tran- 



sitions with good accuracy. Table fO shows the locations of the phase transitions for each 



target volume iV^, with estimates of the uncertainty given in parentheses. 
5 An analogous relation exists for 22(^3,^5). 
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«8(M0 


4(iV 5 ) 


10k 


1.5375(5) 


0.454(1) 


20k 


1.595(1) 


0.469(2) 


30k 


1.640(1) 


0.480(1) 



Table 4.1: Locations of the phase transitions. 

Extensive Monte Carlo simulations were then run near the phase transitions, in an effort 
to minimise statistical errors. These ranged up to 8 million sweeps for N§ = 10k and up 
to 2 million sweeps for iV| = 30k. It was found that c 2 (g, iVs) no longer behaves as an 
order parameter for sufficiently high resolutions of the coupling constants. The scaling of 
C2(g, N5) was found to be similar for both the k% = and k,q = limits. Basically, in the 
crumpled phases c 2 (g,N 5 ) is small and decreases with volume, and the phase transitions 
are characterised by an acute increase in C2(g,N 5 ) which becomes progressively sharper 
with volume. 

The vertex and tetrahedron specific heats were measured close to the phase transitions. 
At this level of resolution, distinct peaks in c 2 (k , iV 5 ) and c 2 (k 3 , N 5 ) are only observed for 
N§ = 30k. In fact, for N§ = 10k, peaks in the specific heats are hardly perceptible. The 
results are shown in table |4.2j . Graphical presentation of the results cannot be justified 
due to the poor statistics. The errors, shown in parentheses, were estimated using the 
binning technique to counter the effect of correlated data (see appendix [B.2j ). It is clear 





C2(K C ,N 5 ) 




10k 


0.09(1) 


1.81(6) 


20k 


0.105(7) 


1.75(10) 


30k 


0.14(2) 


2.0(3) 



Table 4.2: Scaling of the vertex and tetrahedron specific heats. 

that one cannot make any firm deductions regarding A using these results. 

This approach to determining the order of the phase transitions has obviously not proven 
to be very effective. It seems that this type of undertaking is simply too ambitious given 
the computational hardware available. In retrospect, these difficulties should have been 
expected. Fortunately, there exist certain numerical algorithms which can improve the 
efficiency of such studies. For example, histogramming methods have been successfully 
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implemented by other researchers in four dimensions |7T| , |T2| , [H| , did not improve matters 
in our case. This may have been due to an insufficient quantity of raw data. 

Analogous studies in four dimensions have shown that finite size scaling results require 
careful interpretation. A study by Catterall et al. using target volumes of up to N\ = 
8k found A to be 0.259(7), which corresponds to a continuous phase transition [|7TJ. 
Further numerical simulations with target volumes up to N\ = 32k, determined A to be 
0.81(4) ]73]| . This result opened up the possibility that A — > 1 for large iV 4 . It is entirely 
possible that a similar effect could exist in five dimensions. In conclusion, this approach 
is fraught with difficulties mainly relating to the sheer scale of the problem. We shall now 
consider an alternative way of determining the order of phase transitions. 



4.3.2 Metastability 

In this section we exploit the fact that first order phase transitions are characterised by 
metastability. It is possible to observe the co-existence of multiple phases in thermal 
equilibrium from Monte Carlo time series plots. Given sufficient data, one may then plot 
histograms of the data, in which metastability would manifest itself as multiple peaks in 
the distribution. These tests may reveal the existence of first order phase transitions. Of 
course, the absence of metastable effects does not necessarily mean that a transition is not 
of first order. One could envisage cases in which the latent heat is undetectably small Q 
Since latent heat increases with volume, we shall only consider the time series plots of 
our largest target volumes (N§ = 30k). 

Figure |4~6| (a) is a time series plot of Nq taken over 4 million sweeps for kq = 1.640. (It 
becomes immediately apparent that the data is highly correlated - this is evident from the 
vertical streaking effect. However, this is not too important for our current consideration.) 
From closer examination, it appears that the ensemble flips between two bistable states, 
which are highlighted by horizontal lines. This effect is not very clear by any means 
and hence may not be genuine. If this interpretation is correct then the N = 4160 
and iVo = 4290 dashed horizontal lines correspond to the crumpled and elongated phases 
respectively. 

Figure |4.6| (b) is a time series plot of N 3 taken over 4 million sweeps for n 3 = 0.480. In 
this case, we also find evidence of bistability - though much stronger. Here the crumpled 
and elongated phases are highlighted by the N 3 = 113770 and N 3 = 114300 dashed 
horizontal lines respectively. 

In both cases, we find that the ensemble appears more settled in the crumpled phase. 
From this observation, one may cautiously deduce that the exact location of the phase 



6 This was found to be the case for N* 5 = 10k and iV| = 20k. 
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Figure 4.6: Monte Carlo time series plots of (a) the vertex number N for n = 1.640 and 
(b) the tetrahedron number N 3 for k 3 = 0.480. 
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transition are at marginally greater values of kq and K3. Unfortunately, it was not pos- 
sible to extend this study due to time constraints. The CPU (central processing unit) 
time required to complete these time series runs is of the order of 7 weeks - even for our 
most powerful computer. In conclusion, these preliminary results appear to be consistent 
with first order phase transitions. However, the study must be extended to larger vol- 
umes before our claims can be accepted with confidence. At this point we conclude our 
investigation of phase transitions and concentrate on further exploration of C 2 - 



4.4 Two Dimensional Coupling Space C2 

Until now, we have only considered the K3 = and = limits of C 2 , where Sf [T, Ko, K3] 
effectively has only two terms. The next stage in our investigation involves probing regions 
of C 2 where Ko and k 3 are non-zero. Given the difficulties encountered earlier, it is clearly 
not feasible to undertake any form of finite size scaling analysis in C 2 . Our objective is 
simply to detect other non-trivial phase structure. 

It is tempting to speculate that there may exist a line of phase transitions (phase 
boundary) passing through (1.5375, 0) and (0, 0.454) in C%. One may then consider 
whether the phase boundary extends indefinitely or ends in a critical point, as in the 
phase diagram of water. Such question are very difficult, if not impossible, to answer 
predictively. It is not even clear whether such a boundary would be straight or curved. 
The preliminary simulations of section POI showed marked differences between the strong 



coupling phases of the k% = limit and the kq = limit. This may hint at the existence 
of other structure in the region of C 2 where Kq < and K3 < 0, which would explain the 
small peak in c 2 (/t3, N5) found near K3 = —1. Let us now attempt to resolve some of these 
issues. 



We begin by exploring C2 for triangulations with target volumes of N§ = 10k. This 
was done by first fixing k and then tuning k 3 to phase transitions. The first series 
of simulations were run for kq = 1.0, from which a transition was detected at K3 = 
0.157(1). This result effectively proves the existence of a line of phase transitions in C 2 
(see figure fO](a)). The equation of the best fit line was found to be 

k 3 = -0.2955(1)« + 0.4536(1), (4.36) 

with x 2 ~ 1-9 d.o.f. -1 . The good fit indicates that the phase boundary is a straight line. 
The phase boundary was then traced in the same manner at regular intervals ranging 
from kq = —20 to k,q = 20. In each case, k 3 was tuned until the location of the phase 
transitions were known to within at least ± 0.001. The results are shown in figure |4.7| (b). 
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Figure 4.7: (a) Possible line of phase transitions, (b) Phase boundary in C 2 . 
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4.4.1 Phase Boundary 

From figure [4.71 (b), it is evident that the phase boundary does not form a straight line 
in C2. Analysis has shown that the boundary is straight in the approximate range of 
—2 < k < 2. However, significant deviations are observed for large positive values of Kq. 
Careful observation also reveals minor deviations for large negative values of k . If we 
ignore these deviant points then the best fit line is given by 

n 3 = -0.2943(3)k + 0.4564(4), (4.37) 

with x 2 ~ 23 d.o.f. -1 . This line is drawn in figure f4.7| (b). The possibility of an ending to 
the phase boundary remains, as yet, an open question. The curved nature of the phase 
boundary is not currently understood. 

Interestingly, it was noticed that the phase transition appears 'sharper' for large positive 
values of kq. By this we mean that for a given resolution of n 3 , the first derivatives of the 
free energy appear to be more like discontinuous functions of k 3 . We illustrate this point 
using the k = —20 and k = +20 simulations as examples (see table |4~3| ). 



K 


^3 


(Po) 


(R) 


(9) 


+20.0 


-4.425 


0.19185(3) 


3.97006(8) 


40.2(13) 


+20.0 


-4.426 


0.01285(1) 


3.6867(7) 


10.44(1) 


-20.0 


+6.405 


0.1443(6) 


3.812(2) 


21.7(11) 


-20.0 


+6.404 


0.1440(6) 


3.811(2) 


21.6(14) 


-20.0 


+6.403 


0.1431(6) 


3.808(2) 


19.2(7) 



Table 4.3: Measurements taken near phase transitions. 



This may indicate that the phase boundary is more strongly first order for large positive 
values of n . Further exploration of C 2 is temporarily postponed. This chapter is ended 
with another look at the strong coupling phases. 

4.5 Strong Coupling Regimes 

Triangulations in the crumpled phase of the k 3 = limit rapidly thermalise to a particu- 
lar domain of the space of triangulations T, which is asymptotically fixed as — > —00. 
Time series plots of iVo have confirmed this fact. Naturally, one assumes that this do- 
main corresponds to the kinematic lower bound of Nq. Let us explain the basis for this 
assumption. 
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In the k 3 — limit, the action is a function of iVo and N$. When kq is large, tri- 
angulations with large N will dominate the partition function - due to their increased 
statistical weight. Obviously, for a given N 5 , one can identify a subset of T in which N 
and po are maximal. As — ► 00 we expect this subset to saturate the partition function. 
Conversely, one can identify another subset of T in which Nq and po are minimal. One 
expects this subset to saturate the partition function as Ko — > —00. However, certain 
measurements of po for k 3 < K3 have shown that this is not the case (see figure |4.2| (a)). 
This observation raises the question: why does the k — > — 00 limit of p n °t correspond 
to the true kinematic lower bound? 

In constrast, the crumpled phase of the kq = limit is more complicated - in the sense 
that there does not appear to be a single asymptotic domain of T to which the ensembles 
equilibrate. In actual fact, there is evidence of a kinematic lower bound for N 3 and R, but 
triangulations for k 3 ^ -00 do not coincide with the bound as one would expect. What 
could be the reason for this? This effect may be due to ensembles becoming trapped in 
metastable states. 

One expects triangulations with small iV 3 to dominate for k 3 < k^, and those with 
minimal iV 3 to saturate the partition function as K3 — > —00. Our observations may 
be explained by the existence of a number of domains in T having minimal iV 3 . In 
other words, the free energy surface has a number of vacua. Time series plots of N 3 
seem to be consistent with this theory - we find that during thermalisation iV 3 often 
reaches apparently stable equilibria followed by sudden drops. These may translate to 
the ensemble migrating from one minimum of the free energy to another via random 
fluctuations. It is possible to test for the existence of such metastable states by simple 
experiments. This is done by thermalising a set of different starting configurations. If the 
final state varies, then this is indicative of metastable states. This avenue of research was 
not explored, due to time constraints. 

Using equation ( j2.20| ) one can write down an expression for the change in S§[T, n , k 3 ] for 
any five dimensional (k, I) move. In the k 3 = limit, it is written as 

AS 9 5 [T, kq, 0] m = k 5 (1 — k) — K ( l C 6 - fc C 6 ) = -ASf [T, «„, 0] (I)fc) . (4.38) 

Numerical simulations have shown that A5f [T, k , 0]tk,i) is positive for k = 1,2,3 and 
hence negative for k = 4,5, 6. This is true for all Kq. In the kq = limit, the analogous 
relation reads 

AS 9 5 [T,0,K 3 ] m = k 5 (1 - A;) - k 3 ( 1 C 3 - k C 3 ) = -ASi[T,0,K 3 ] m . (4.39) 
In this case we find that for k 3 = —10, the change in action ASf-[T, 0, ft 3 ](fc,z) is positive for 
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k — 1,2,4 and hence negative for k = 3, 5, 6. Of course, when Kn = k 3 = the change in 
action is positive for k = 1,2,3. This means that ASf [T, 0, «3](3,4) changes sign at some 
point in the range —10 < k 3 < 0. By substituting k = 3 and I = 4 into ( [4.39[ ) we get 



AS![T, 0, «3] (3)4) = k b - 3k 3 = -AS! [T, 0, k 3 ]( 4j 3)- (4.40) 

In other words, the sign change occurs when = 3k 3 . Figure |4.8| plots the tuned simplex 
coupling k 5 versus K3 (for kq = 0). The straight line corresponds to k 5 = 3k 3 . Clearly, 
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Figure 4.8: Tuned simplex coupling k 5 versus tetrahedron coupling K3. 

we see that K5 < 3K3 for K3 < — 5. It is exactly this region where we observe the strange 
metastable effects. In conclusion, it appears that the sign change of ASf [T, 0, K3] ( 3) 4) 
results in a destabilisation effect, which seems to prevent the ensemble from reaching true 
equilibrium. 

Notice from figure |4.8| that k 5 = 3% near k 3 = —5. Recall that k 5 is calculated 
via an iterative tuning procedure during the thermalisation stage. The change in action 
due to the (3,4) move for K3 = —5 is zero when Kb = 15. In this case, k 5 was tuned 
to —14.992139. Clearly, these values satisfy K5 > 3^3. This accounts for the apparent 
complete thermalisation. This simulation was repeated as an experiment. This time 
k 5 was tuned to —15.138613. As expected, the ensemble became trapped in metastable 
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states. The critical simplex coupling k 5 obviously plays a crucial role in determining 
whether the ensemble reaches true equilibrium. This type of effect will always emerge 
whenever the overall change in action of a (k, I) move is close to zero. 



Chapter 5 

Fractal Geometry and Stacked 
Spheres 

5.1 Spatial Distributions 

far, our extensive series of Monte Carlo simulations have revealed the existence 
of a phase boundary and at least two phases in C 2 - The weak coupling phase 
corresponds to elongated triangulations with large scalar curvature and vertex 
density. In contrast, the strong coupling regimes of the K3 = and kq = limits 
appear to be distinct. The main purpose of this section is to characterise these phases 
in geometrical terms. This goal may be achieved by visualising the the S 5 triangulations 
at a global scale; in other words, as structures embedded in M 6 . Specifically, this entails 
measuring the number of simplices N at a given geodesic distance r from an arbitrary 
origin simplex^]. The resulting spatial distribution N(r) is normalised and averaged over 
all possible origin simplices and over many Monte Carlo configurations. Therefore, the 
mean of the distribution is, by definition, located at r = (g), the mean geodesic distance. 
From such information one can gain knowledge of the spatial distribution of simplices 
and, as we shall see later, learn more about the geometrical and fractal nature of the 
triangulations. This will give us a more complete picture of the phase structure as a 
whole. 



The first step of this project is to establish how N(r) varies with kq and K3. Measurements 
of N(r) were taken at regular intervals over 10 5 Monte Carlo sweeps, for target volumes 
of iVg = 10k. Simulations were run over a wide range of couplings, which included points 
near the phase transitions: = 1.5375 and K3 = 0.454. Figures |5.1| (a) and |5.1| (b) show 



1 Dctails of the algorithm are provided in appendix A. 3 
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plots of the mean spatial distribution (N(r)) for the k 3 = and kq = limits respectively, 
where each curve is labelled by its corresponding coupling. Both graphs are plotted with 
the same scale to allow comparisons to be made. Measurements of N(r) for k = k 3 = 
serve as a point of reference. The error estimates^ (indicated by dots) are considered 
reliable as they are compatible with the fluctuations in N(r). 



5.1.1 Results and Interpretation 

From figures |5.1] (a) and |5.1| (b), we find that the weak coupling spatial distributions are 
virtually indistinguishable. This reaffirms the fact that the weak coupling regimes of the 
K3 = and /to = limits correspond to the same elongated phase. (Some simplices of 
these triangulations are separated by geodesic distances of up to 140 units.) One can 
plainly see that the spatial distributions near the phase transitions are also very similar. 
This could be consistent with the phase transitions being of the same order, as originally 
suspected. Perhaps one should not read too much into this observation as there are hidden 
errors derived from the uncertainty in the location of the phase transitions. 

Let us now consider the spatial distributions in the strong coupling regimes. It is 
here that we find striking differences between the two sets of distributions. The k,q *C 
plots of figure |P|(a) are consistent with the existence of an asymptotic set of crumpled 
triangulations, and correlate perfectly with measurements of (po( K o)) and (R(kq)) from 
section [4.2. 1|. These results highlight the distinction between the two strong coupling 



regimes. From figure |5.1| (b) it is clear that k 3 = —5.0 plot corresponds to highly crumpled 
triangulations. This is also in agreement with earlier results. The most surprising result of 
this study relates to the k 3 = —10.0 plot, which has two peaks! This spatial distribution 
of simplices accounts for the unexpected increases observed in p , R and g for k 3 < — 5. 
These represent the possible metastable configurations discussed in section |4~5| . 

One may wonder what the corresponding triangulations look like. Unfortunately, it is 
not possible to uniquely determine the geometry of the triangulations from knowledge of 
N(r) alone. Despite this fact, one can deduce certain features from basic arguments. First 
of all, the double peak tells us that the geometry is not even remotely spherical. One can 
only conclude that the geometry is comprised of a number of 'baby universes' connected 
by thin necks. The last and perhaps most powerful deduction is that the number of these 
necks is small - probably just one or two. This becomes clear when we consider the spatial 
distribution of the elongated phase which are believed to be branched polymers. As these 
triangulations have a large number of necks, the mean spatial distribution is smoothed 
out, resulting in a single peak. 



The errors were estimated using the Jackknife method, which is explained in appendix B.3 
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Figure 5.1: Mean spatial distribution (N(r)) for (a) the k% — limit and (b) the kq = 
limit, for triangulations with target volumes = 10k. 
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5.2 Fractal Structure 

This section deals with the fractal nature of five dimensional triangulations. Our primary 
aim is to prove that the elongated phase is characterised by branched polymer triangula- 
tions. This is done by evaluating the internal Hausdorff dimension dn as a function of r. 
In general, the internal Hausdorff dimension of a geometrical structure is defined in terms 
of the volume V contained within a radius r. 

Vocr dH (5.1) 

This definition may be naturally extended to triangulations, by considering the number 
of simplices contained within a given geodesic distance r of an arbitrary origin simplex. 
Now, the internal Hausdorff dimension of a triangulation is a function of r, which is given 
by 

V{r) = ar dH{r \ (5.2) 

where a is a constant and V(r) is the volume of a <i-ball with geodesic radius r. This 
expression is valid only for a limited range of r, since we are dealing with finite triangu- 
lations. Of course, N(r) is the volume of a shell at a geodesic radius r. Therefore, V(r) 
is given by the sum over N(r) from to r. By averaging over all origin simplices, one 
obtains the following expression. 

(V(r))=j2(N(i)) (5.3) 

i=0 

This is calculated in practice using the results of section |5.1| . By taking logs and differen- 
tiating (|5.1|) with respect to lnr, we obtain a formal definition of the internal Hausdorff 



dimension. 

d H = — 5.4 

dmr 

This may also be extended to (discrete) triangulations by replacing derivatives with differ- 
ences. Using (|5.3|) one can define an expression for the mean internal Hausdorff dimension 
function. 

(d H (r)) = HV(r))-HV(r-l)) ^ 
In r — ln(r — 1) 

Using equation ( |5.3|) one can then easily express <i#(r) in terms of N(r). 



5.2.1 Branched Polymers 

Figures |5T2|(a) and [53(b) show plots of {dn{r)) for the K3 = and k q = limits respec- 
tivelyQ. The internal Hausdorff dimension of the triangulations is given by the maximum 
3 The method of error estimation is discussed in appendix 
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of dn{r), which widens into a plateau as the target volume is increased. From these fig- 
ures one can plainly see that the fractal dimension of the weak coupling triangulations is 
du ~ 2. This result proves that the weak coupling triangulations are branched polymers. 
The strong coupling triangulations are characterised by large internal Hausdorff dimen- 
sion. It is important to realise that these results are subject to finite size effects; one can 
learn more by measuring <i#(r) as a function of volume. 

Further Monte Carlo simulations were run in order to establish the scaling behaviour 
of dn(r). The study was repeated for triangulations of target volumes of = 20k and 
iVg = 30k. It was found that du m 2 in the weak coupling phase, regardless of volume. 
From now on we shall refer to the common weak coupling phase as the branched polymer 
phase. In the strong coupling regime of the k$ = limit, we find that du increases 
approximately linearly with volume. This would be consistent with du — > oo in the 
thermodynamic limit, as conjectured by Ambj0rn |54 . Double peaks in N(r) were also 
detected for K3 <C k%. Finally, the internal Hausdorff dimensions measured near the phase 
transitions also increase with volume. However, these results are not considered reliable 
due to relatively poor statistics and hidden errors due to the uncertainty in the location 
of the phase transitions. 



5.3 Stacked Spheres 

For a given space of triangulations T with fixed volume, it is clear that one can identify 
a subspace of T in which Nq is minimal. These triangulations correspond to the lower 
kinematic bound of Nq. Conversely, there also exists a subspace of T in which Nq is max- 
imal. These triangulations of the upper kinematic bound are known as stacked spheres. 
This section deals with the analytic study of these structures. The name given to these 
structures is quite apt when one considers how they are constructed. Minimal d-spheres 
are composed of d + 2 vertices and d + 2 simplices pairwise connected along their faces. 
Stacked ^-spheres are generated from minimal <i-spheres by performing only vertex inser- 
tion moves. In a sense, new simplices are 'stacked onto' the sphere. It is the fact that only 
(1, d + 1) moves are used that allows us to study stacked spheres by analytical means. 

It is a trivial matter to calculate the scalar curvature R of a stacked sphere. Consider a 
minimal 5-sphere, which has N$ = 35 and N$ = 7, 4 . Using equation ( |2.20|) one can show 
that AN$ = +10 and AN 5 = +5 for the (1,6) move. These results give us the following 
relations for N% and N$ after n vertex insertion moves, where n£N. 



jV 3 = 35 + 20n 
4 The number of i-(sub)simpliccs of a minimal d-sphere is given by d+2 Ci+i. 



(5.6) 
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N 5 = 7 + 5n (5.7) 

By eliminating n and dividing by N 5 we get the following equation for the scalar curvature, 
expressed in terms of N 5 . 

iV, 7 



Clearly, we find that R — > 4 in the thermodynamic limit. One must note that (|5.8|) does 
not hold for all volumes - only those that satisfy (|5.7|). Substitution of N 5 = 10,002 
into ( |5.8| ) gives R ~ 4.0007. This value of the scalar curvature matches exactly the 
weak coupling data presented in figure P~2|(b). This proves that the partition function 
in this limit is saturated by stacked spheres. It is important to bear in mind that these 
structures are branched polymer stacked spheres. These are constructed by performing 
vertex insertion moves on randomly chosen simplices. One can generate singular stacked 
spheres if the moves are performed on simplices in a specific pattern. This resulting 
triangulation has a vertex with a very large local volume. Such structures clearly do not 
resemble branched polymers. 

5.4 Walkup's Theorem 

It is the fact that stacked spheres represent a kinematic bound that allows for the existence 
of inequalities relating Ni and No, at least in three and four dimensions. Collectively 
the inequalities are called Walkup's theorem ||74|| . These results have been central to 



certain advances made in the understanding of the weak coupling limit of dynamical 
triangulations []75|, [7(| . 

d = 3 : Ni > AN — 10 (5.9) 
d = A: Ni > 5iVo-15 (5.10) 



Inequalities ( |5.9| ) and ( |5.10| ) are true for all spheres and hold with equality only for stacked 
spheres 

5.4.1 Generalisation of Walkup's Theorem 

This chapter is concluded with a simple proof of Walkup's theorem by considering the 
effect of (k, I) moves on stacked spheres. In fact, the proof is a generalisation of Walkup's 
It is worth noting that there exist four dimensional triangulations which satisfy the lower bounds 



of inequality (5.10) that are not stacked 4-spheres. This is because in even dimensions there are (k, I) 
moves that transform triangulations without changing the /-vector. In other words, it is a necessary 
but insufficient condition for stacked 4-spheres to satisfy the lower bounds of Walkup's theorem in four 
dimensions. This does not happen in odd dimensions. 
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theorem to arbitrary dimension. In one and two dimensions Nx is uniquely denned for 
any given N . 

d = l: Ni = N (5.11) 
d = 2 : Ni = 3N - 6 (5.12) 

The following analysis applies to all dimensions d > 0. From chapter [l] we know that a 
minimal <i-sphere is composed of d+2 C i+1 i-(sub)simplices. 

N = d+2 d = d + 2 (5.13) 
Nx = d+2 C 2 (5.14) 

Starting from a minimal <i-sphere one can construct a stacked <i-sphere by performing a 
series of n vertex insertion moves. Again, using equation ( p. 20 ) one can calculate AN 
and ANx for a (1, d + 1) move. 

AN = +1 (5.15) 
ANx = +(rf+l- 1 C d ) (5.16) 

Using these equations one can write very simple expressions for iV and Nx in terms of 
just d and n. 

No = d + 2 + nANo = d + 2 + n (5.17) 
Nx = d+2 C 2 + nANx = d+2 C 2 + n(d+l- l C d ) (5.18) 

By eliminating n one can express Nx in terms of only iVo and d. The following equation 
is true for any stacked c/-sphere. 

Nx = d+2 C 2 + (N -d-2)(d+l- ^d) (5.19) 



Equation ( |5.19| ) clearly reduces to ( |5.11| ) and ( |5.12|) when d = 1 and d = 2 respectively. 



It is conjectured that the following inequality holds for any rf-sphere. 

Nx > d+2 C 2 + (N - d - 2)(d + 1 - ^d) (5.20) 



This inequality is offered as a generalised form of Walkup's theorem. It clearly reduces 
to Walkup's theorem when d = 3 and d — 4. This is taken as evidence supporting the 
validity of our conjecture. The following sections form a proof of this result. 
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5.4.2 Proof Strategy 

The strategy for proving our generalisation of Walkup's theorem is as follows. We begin 
with a triangulation that satisfies (|5.19|) , i.e. a stacked sphere. It is known that a stacked 
sphere may be transformed into any other combinatorially equivalent (piecewise linear 



homeomorphic) triangulation using a finite series of (k, I) moves [42, f|3). If one can show 



that the effect of each (k,l) move is consistent with ( |5.20| ), then one may be able to 
conclude that the lower bound holds for all <i-spheres. 



From now on we only consider d > 1. Our conjecture reads 

Ni > d+2 C 2 + (N -d-2)(d+l). (5.21) 



Inequality (|5.21 ) may be shown to be false if there exist (k, I) moves that have 8 < d + 1, 

(5.22) 



where 5 is defined as follows. 

= AN, d+2 - k C d - k C d 
AN d + 2 ~ k C d+1 - k C d+l 



For k = 1 (vertex insertion moves) and k — d + 1 (vertex deletion moves) we see that 
S = +{d + 1) which means that the moves are consistent with inequality (|5.21|) . Now 
AN = for all moves with 1 < k < d+1. Therefore the only moves that are inconsistent 
with the inequality are those with 1 < k < d+1 andANi < 0. Using AN% = d+2 ~ k C d - k C d 
we see that AA^ < only when k = d, d + 1. Therefore 8 < d + 1 only for the (d, 2) 
move. Such a move could potentially be used to generate triangulations that do not satisfy 
( |5.21| ). This situation presents us with two possibilities. Either 



• (a) our conjecture is wrong, or 

• (b) W( 2 ,d) > Af(d,2), 

where J\f(k,i) is the total number of (k, I) moves required to generate an arbitrary triangu- 
lation from a stacked sphere. Option (a) seems unlikely because (|5.21|) is known to reduce 
to Walkup's theorem for d = 3 and d = 4. In the next section we hope to convince the 
reader that (b) is true. 



5.4.3 Clusters and (d,2) Moves 

From here on, it is necessary to visualise simplices in their dual form - details of which are 
to be found in chapter ^|. The reader is reminded that simplices and faces are represented 
by nodes and edges respectively. In <i-dimensions each node has d+1 edges connecting 
it with d+1 other nodes such that pairs of nodes are connected by only one edge. This 
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condition imposes that simplices are pairwise connected. In this chapter, these ideas are 
extended to construct entire triangulations. Their dual versions are called graphs. Let us 
begin by considering the graphs of minimal 3-spheres and minimal 4-spheres. These are 
shown in figure |5.3| , f\. 




Figure 5.3: Graph of a minimal 3-sphere (left) and a minimal 4-sphere (right). 



Before going any further, we must introduce the concept of a cluster. An n-cluster in 
(i-dimensions is a set of n nodes pairwise interconnected by n C2 edges. Each node in the 
cluster has d — n + 2 free edges linked to neighbouring clusters. In general, an n-cluster 
contains n Cj z-subclusters, for i — 2, . . . , n— 1. An z-cluster may be regarded as the dual of 
a (d+1 — z)-(sub)simplex star. The action of a (k, I) move in the dual picture is therefore 
the replacement of a fc-cluster with an Z-cluster. For stacked ci-spheres, (^-clusters can only 
exist as subclusters of (d + l)-clusters, because only (l,d+ 1) moves are allowed. This 
is evident from figure |5.4j which shows graphs of stacked 3-spheres and stacked 4-spheres, 



cf. figure 5.3. Their constituent 4-clusters and 5-clusters are labelled. 



In the dual picture, a (2, d) transformation may be considered as the replacement of a 
2-cluster with a d-cluster. A 2-cluster is the structure formed by two connected nodes. 
Each node has d free edges since the total number of edges is d + 1. The 2d free edges 
are pairwise connected to a = 2d neighbouring nodes. Now, each node of a <i-cluster has 
(3 = 2 free edges. The fact that a (3d means that the nodes of the d-cluster can be 
pairwise connected with the neighbouring nodes []. In other words, it is always possible 



6 The graphs are drawn symmetrically merely for clarity. 

7 This argument only applies to d > 2 since (2, d) moves are not denned for d < 3. 
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Figure 5.4: Graph of a stacked 3-sphere (left) and a stacked 4-sphere (right). 



to perform a (2, d) move whilst maintaining the manifold structure of the triangulation. 
Figure |5.5| illustrates this argument schematically for the four dimensional case. 



>f >f >f >f >f >f >f ^neighbouring ^^^^^^^^ 



nodes 



a = 8 




(2,4) 



2-cluster 



4-cluster 




(3 = 2 



Figure 5.5: Dual representation of the (2,4) move for a stacked 4-sphere. 

Let us now consider (d, 2) transformations. (These moves cannot be performed on 
minimal <i-spheres as > d + 2 for <i-spheres.) Each node of a (d + l)-cluster has 
one free edge connected to another node outside the cluster. By removing d nodes of 
a (d + l)-cluster we are left with one neighbouring node with a' = d free edges and 
d other free edges from d neighbouring nodes. These are to be connected with the 2 



nodes of the replacement 2-cluster. Again, figure |5.6| is a schematic illustration of the 
four dimensional case. Clearly, one cannot pairwise connect the nodes with this edge 
configuration. Forced connection results in multiple edges, which signifies the breakdown 
of the manifold structure. This is because a < (3d. In other words, (d, 2) moves cannot 
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Figure 5.6: Dual representation of the (4, 2) move for a stacked 4-sphere. 

be performed on stacked <i-spheres[| Logically, we conclude that a (d, 2) move can only 
be executed as the inverse of a (2, d) move. Therefore <i-clusters only exist as subclusters 
of (d + l)-clusters or as the end product of (2, d) transformations. For these reasons 
we conclude that N(2,d) > ^{d,i)- This fact proves that our generalisation of Walkup's 
theorem is correct. □ 



The dual formulation of (A;, I) moves were (apparently) first employed by Godfrey and 
Gross p5f to aid their visualisation in dimensions greater than two. It seems, however, 
that this is first time the dual graphs have been used to visualise entire triangulations and 
solve problems in this manner - although their use was recognised by Gabrielli |75] . This 
work has demonstrated the power of simple graph theory concepts to produce analytic 
results in dynamical triangulations. The actual generalisation of Walkup's theorem is not 
in itself a significant result. It is hoped that more complex applications of graph theory 
may lead to new analytic results. 



In fact, it can be shown that this is the case for all moves with I < k < d + 1. 



Chapter 6 



Singular Structures 



6 . 1 Background 




URING 1995 a new feature of dynamical triangulations was discovered by 



Hotta et al. [j3q1 . This development centred on measurements of the vertex 
local volume distribution J\f(no) of four dimensional triangulations. They 
found that triangulations in the strong coupling phase are characterised 
by a small number of vertices with unusually large local volumes. Here we present a 
short historical review of the discovery and subsequent developments. 

Figure is a schematic normalised log-log graph of A/"(n ) for four dimensional 
triangulations with target volumes of N\ = 32k, drawn from results presented by Hotta 



et al. ||35|| . It shows three plots which correspond to the crumpled phase (k 2 = 0.0), the 
branched polymer phase (k,2 = 2.0) and to a point near the phase transition («2 = 1.267). 
For k 2 > the distribution is as one may naively expect: the vast majority of vertices 
have smallQ local volumes (n ~ 10 1 ) and J\f(n ) rapidly tails off to zero for n Q ~ 10 3 . 
In contrast, the distribution of the strong coupling phase is markedly different; the main 
distinction being the existence of a small number of vertices with very large vertex local 
volumes (about a third of the triangulation volume). It was confirmed that this peak in 
A/"(too) is representative of two vertices with approximately equal local volumes n™ ax |35 



Hotta et al. also looked at how A/"(n ) scales with volume [|3l|. In the weak coupling 
phase, there was found to be no significant dependence of M{uq) on volume. However, 
for k 2 < K\ y it was found that n™ ax increases linearly with volume. As a consequence, 
these vertices were dubbed singular vertices since their local volumes seemed to diverge 
in the thermodynamic limit. This was perceived to be a potential problem in reaching a 
continuum limit. For this reason, Hotta et al. constructed a new variant of the discretised 



x It is easy to show that no > d + 1 for simplicial d-manifolds. 
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Figure 6.1: Vertex local volume distribution jV(no) for four dimensional triangulations 
with target volumes N\ = 32k. The singular vertices have local volume n™ ax of order 10 4 . 



Einstein-Hilbert action S[T, k 2 \ by introducing an additional term which inhibited the 
formation of singular vertices. The new action is written as follows. 



S[T,k 2 ,\)=S[T,k 2 ) + ^ £ (n(a°)-k) 2 

1 CT o eT 



(6.1) 



The extra term is a summation over all Nq vertices a of triangulation T. Here, A is a 
coupling parameter and k is a constant, arbitrarily chosen to be 5, since n > 5 in four 
dimensions. Monte Carlo simulations were performed in order to establish the effect of A 
on the singular vertices. It was found that the singular vertices vanish when A > 10~ 3 . In 
fact, the phase transition also disappears and the dominant triangulations are branched 
polymers for all values of k 2 . For A < 10" 6 , the additional term has negligible effect and 
there is no suppression of singular vertices. From this investigation, it became evident that 
singular vertices are an important feature of four dimensional dynamical triangulations, 
and may even drive the phase transition fl35| . 

In a later paper by Hotta et al. (with the same title), it was claimed that no ob- 



servations were made of links with unusually large local volumes |77]. Furthermore, it 
was found that (^-dimensional triangulations (for Kd- 2 = 0) have d — 2 singular vertices, 
at least for 3 < d < 7. It was thought that perhaps these singular triangulations may 
be the result of the ensemble becoming trapped in metastable states. This possibility 
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was explored by running Monte Carlo simulations with three very different starting con- 
figurations. It was found that in each case the ensemble thermalised to the same set 
of singular triangulations. This was taken as evidence to suggest that singular vertices 
are not a result of incomplete thermalisation, but a genuine feature of the model [[77]. 
It was also discovered that toroidal triangulations in the strong coupling phase are also 
characterised by singular vertices. In fact, the Af(no) distribution for T 4 was found to be 
virtually identical to that for S 4 (at least for large volumes). On the basis of this fact, one 
could tentatively speculate that the existence of singular vertices might be independent 
of topology. 

6.1.1 Further Developments 

Further research was carried out by Catterall et al. on singular vertices , which we shall 
now review. Monte Carlo simulations were run in four and five dimensions using S[T, kq] 
with kq = 0. This choice of coupling ensured that all triangulations carried equal weight in 
the partition function^. On the basis of their numerical results it was conjectured that such 
(i-dimensional triangulations are characterised by a single singular (d— 3)-subsimplex. This 
primary singular subsimplex is the product of d ~ 2 Cj + i secondary singular i-subsimplices, 
for i = 0, . . . , d — 4. These results appear to be in contradiction with those of Hotta et 
al, who claim that no singular links were observed in four dimensions 1)771 . However, 
this is probably not the case because Hotta et al. describe a link as singular if it has an 
unusually large local volume. In this thesis - following the convention set by Catterall et 
al. - a singular link is defined as the product of two singular vertices, regardless of its 
local volume. Driven by numerical results, it was conjectured that the local volume of 
primary singular subsimplices increase with volume as 

rC™<xN d l (6.2) 

In addition, it was conjectured that all secondary singular i-subsimplex local volumes 
scale linearly with volume, for i = 0, . . . ,d — 2. All multiple singular structures of a given 
dimension were confirmed to have approximately equal local volumes. 

Dual-Spheres and Local Entropy 

These observations were explained heuristically in terms of entropy. Consider an i- 
subsimplex a 1 of a ^-dimensional triangulation T. The star St{a l ) of a % is homeomorphic 
to B d , because T is a simplicial d-manifold. Since a % has % + 1 vertices, each simplex of 

2 In five dimensions, S[T, kq] is only equivalent to the discretised Einstein-Hilbert action S[T, K3] when 
kq = and K3 = 0. 
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St(a l ) has d — i vertices that are not in common with cr\ These structures are a set of 
(d — i — l)-subsimplices which form an S^ 1 ^ 1 submanifold of T. Moreover, this structure 
is a submanifold of the S^ 1 boundary of St^a 1 ). The S ,d ~*~ 1 submanifold is called the 
dual-sphere of cr\ The local volume n(a % ) of o % is therefore equal to the volume of its 
dual-sphere. To a given i-subsimplex a 1 , one may associate a local entropy u)(o~l). This is 
equal to the entropy of its dual-sphere, which is denoted by n(S 

Let us now consider some examples, beginning with i = d — 2. In this case, the 
dual-sphere of a (d — 2)-subsimplex a d ~ 2 is a 1-sphere. Clearly, the S 1 dual-sphere has a 
unique triangulation for a given volume (up to automorphisms). This means that the local 
entropy u(a d ~ 2 ) does not increase with volume. Consequently, the formation of singular 
(d — 2)-subsimplices are not entropically favoured. Now consider the S 2 dual-sphere of a 
(d — 3)-subsimplex a d ~ 3 . The entropy of a 2-sphere is known to grow exponentially with 
volume fl57|. This fact implies that the local entropy of o~ d ~ 3 also increases exponentially. 
This accounts for the existence of singular (d — 3)-subsimplices. For singular subsimplices 
of lower dimension (i < d — 3), the corresponding dual spheres have dimension d > 2. The 
entropy of these structures are known to increase (at least) exponentially with volume. 
Hence, for the same reasons, this explains the existence of singular i-subsimplices, for 
i — 0, . . . , d — 4. The problem with this picture is that it cannot account for the apparent 
non-existence of singular vertices in three dimensions. According to the reasoning of 
Catterall et al, vertices (in three dimensional triangulations) with large local volumes 
should be entropically favoured since their dual-sphere are 2-spheres. This point seems 
to be a problem and is recognised as being so ]7jJ. 

This brings us to the question: why should there be only one singular primary sub- 
simplex? Catterall et al. answer this important question with the following argument. 
Consider a triangulation with a number of randomly distributed singular vertices. Each 
vertex strives to maximise its local volume. This can be done if the vertices tend to share 
their local volume. In other words, the dual-sphere of a singular vertex will itself contain 
other singular vertices. Eventually some singular vertices will dominate by 'absorbing' 
others. This coalescence continues until only d — 2 vertices are left. These maximise their 
local volumes by sharing simplices and hence form a primary singular subsimplex. 

It is quite easy to determine the intersection of these dual-spheres. Consider the four 
dimensional case, where two S 3 dual-spheres whose intersection forms a 2-sphere. This 
can be shown using simple geometry. Of course the situation becomes more complicated 
in higher dimensions. In five dimensions we have three 4-spheres dual to the singular 
vertices. Again, simple algebra shows that each pair of 4-spheres intersects at a 3-sphere 
dual to a singular link. In fact, one can show that all three 4-spheres intersect at a 2-sphere 
dual to the singular triangle. In higher dimensions, one can follow the same pattern to find 
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that the intersection submanifold is always S 2 , the dual-sphere of the primary singular 
subsimplex. 



6.2 Three Dimensions 

According to the entropy arguments proposed by Catterall et al., one would expect to 
observe singular vertices in three dimensions; since the dual-sphere of a vertex in three 
dimensions is a 2-sphere, whose entropy is known to grow exponentially with volume [5?|. 



None have ever been detected JTBj . The purpose of this section is to resolve this apparent 
discordance. Let us begin by repeating the work of Catterall et al, as their numerical 
results were not published in explicit form. 

A new function incorporated into the standard dynamical triangulations code allowed 
for the measurement of the vertex local volume distribution AA(no) and hence detection of 
any singular vertices that may exist (for d < 6). The algorithm is described in appendix p9|. 
Monte Carlo simulations were run for target volumes up to N§ = 30k and over a wide 
range of kq. In each case, the maximum vertex local volume n™ ax was measured 10 3 
times over a period of 10 5 sweeps, yielding a mean value (n™ ax ). The results presented in 
figure BT2| clearly show that singular vertices are not found in three dimensions. 



6.2.1 Non-existence of Singular Vertices 

In this section, our aim is to construct a plausible explanation of why singular vertices are 
not prevalent in three dimensional triangulations. The following arguments are developed 
from those put forward by Catterall et al. WE . 



Consider a simplicial 3-manifold T. Now identify an arbitrary B 3 submanifold of T with 
volume N$. The 3-ball submanifold has a 2-sphere boundary with volume N 2 . Let the 
space of triangulations of the 3-ball be symbolised by T^-baii- There exists a special 
subspace of T^-baii whose constituent triangulations have N 2 = N 3 , at least for N 3 > 4. 
The triangulations of this subspace, denoted by T sing fl, are described as vertex- singular, 
since every tetrahedron of the 3-ball is common to a single vertex. In other words, we 
have T sing fi C T 3 - ba u- Furthermore, it is clear that the 2-sphere boundary is the dual- 
sphere of a vertex iff N 2 = N 3 . One can now define the subspace of non-vertex-singular 
triangulations T non as follows. 

Tsingfi U T norlj o ^ 1~3-b a u (6-3) 



This space of 3-balls are obviously triangulated such that N 2 ^ N 3 . 
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Figure 6.2: Maximum vertex local volume (n™ ax ) versus kq for three dimensional trian- 
gulations. (The lines serve only to guide the eye.) 



Now one must consider how the entropy of these structures grow with volume. The 
entropy of 3-balls is given by f^T^aiz) — Q{B 3 ). The entropy of vertex-singular tri- 
angulations effectively grows as that of 2-spheres because N 2 = iV 3 ; therefore we have 
Q(%3ing,o) = Q(S 2 ). The entropy associated with these spaces of triangulations can be 
equated to give 

^(Tsing.o) + ^(^non,o) = ^ {T^-ball) > (6-4) 

which may be rewritten as 

Q(S 2 ) + Q(T nont0 ) = Q(B 3 ). (6.5) 
The observed non-existence of singular vertices for k,q = would therefore be due to 

^(Xion,o) > tt(T singt0 ). (6.6) 
Using equation ( |6.5|) one can re-express this inequality as 

Q(B 3 ) >2fi(S 2 ). (6.7) 
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For a given simplicial <i-manifold of fixed volume Nd, there exists an upper kinematic 
bound to the local volume of a vertex, which is a function of Nj. We shall call this upper 
bound riQ tn . The non-existence of singular vertices may then be explained by comparing 
the entropy of 3-balls and 2-spheres with volume n^ n according to inequality (|6.7|). Let 
us illustrate these ideas with an example. Consider a three dimensional triangulation of 
volume N 3 = 1000. In this case n($ m could be, say, 400. That is to say: it is geometrically 
impossible for any vertex to be common to more than 400 simplices. If the entropy of an 
arbitrary 3-ball of volume N$ = 400 is much greater than twice the entropy of a 2-sphere 
of volume N2 = 400, then vertex-singular triangulations (with volume N$ = 1000) will 
not contribute significantly to the partition function. Obviously, in order to extend this 
idea to the thermodynamic limit, one must consider the growth rates of entropy. 

Although it may be possible to prove inequality (|6.7| ) in principle, it is likely to be 
very difficult We therefore offer the result as a conjecture. One can, however, show that 

Q(B 3 )^Q(S 2 ). (6.8) 

This is seen by realising that the boundary of a 3-ball is a 2-sphere. In this case the space 
of triangulations of S 2 may be regarded as a special subspace of the space of triangulations 
of B 3 , cf. Tgingfi C T 3 - ba u. From this fact one can deduce (|6.8| ). This inequality may be 
trivially extended to arbitrary dimension d giving 

fi(5 d ) > ft^ 1 ). (6.9) 

By combining inequalities ( |6.7| ) and ( |6.8| ), we find that singular vertices can only exist in 
three dimensions if 

2Q{S 2 ) > Q(B 3 ) > Q(S 2 ). (6.10) 

In order for (|6.10|) to hold in the thermodynamic limit, the entropy functions of S 2 and 
B 3 must grow at the same rate in this limit. This seems unlikely. 



Our new way of looking at the problem involves comparing the entropy of the singular 
vertices with those of 'ordinary' triangulations of B 3 , and determining which is the greater. 
This is in contrast to the approach taken by Catterall et al. who only considered the 
entropy of the dual-spheres. These arguments are valid for k,q = 0, since all triangulations 
are equally weighted in the partition function. For k 7^ 0, one must consider the effect of 
the action. Figure |6.2| shows that vertices in the weak coupling phase have slightly lower 
local volumes, because the action forces extra vertices into the triangulations. 



3 Depending on the growth rates of fl(B 3 ) and fl(S 2 ); inequality (6.7) may not be true for all N3. As 
physicists, however, we are only really interested in the thermodynamic limit. 
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6.3 Four Dimensions 

The next step in our investigation of singular structures takes us to four dimensions. 
All studies published to date have concentrated on k = or regions near the phase 
transition. For a more complete picture, it was decided to measure the vertex local 
volume distribution J\f(no) across a much wider range of couplings, including negative 
values of k . For the time being we shall focus on the special case of k = 0, and attempt 
to reproduce the findings of previous studies. Measurements were recorded of the two 
largest vertex local volumes and the corresponding link local volume. If the two vertices 
did not form a link in the triangulation, then the link local volume was zero by definition. 
For each Monte Carlo simulation measurements were taken at 10 2 sweep intervals over a 
period of 10 5 sweeps, for target volumes up to N\ = 40k. 



Singular Vertices 

As expected, two singular vertices were detected for each target volume with k = 0. 
These are labelled as and cr° 2 . They were found to have approximately equal local 
volumes, which consistently differed slightly. For this reason, the two quantities were 
distinguished by calculating two separate averages. The largest and second largest vertex 
local volumes are denoted by ^(cr^) and n(a° s2 ) respectively. Figure |0|(a) plots the local 
volume of the two singular vertices versus Ko for each target volume. The k,q = data 
was fitted to the following curves using the chi-squared method, which is described in 
appendix [B.4. 



n(a° sl ) = ai + 6iA^ 4 Cl (6.11) 
n{a° s2 ) = a 2 + b 2 N 4 c * (6.12) 

The best fit curves (both with x 2 ~ 5 d.o.f. -1 ) have coefficients c\ = 1.017(2) and c 2 = 
1.006(6). Again, as expected, these results are consistent with the conjecture that the local 
volumes of secondary singular subsimplices diverge linearly with triangulation volume. 
Furthermore, the two singular vertices were found to be common to approximately one 
third of the triangulation volume, since b\ = 0.307(6) and b 2 = 0.34(2). This is also in 



perfect agreement with the work of Hotta et al. p5| , [77fl and Catterall et al. |78 



One should bear in mind that ( |6.11D and fl6.12|) are valid only for large volumes, since 



the relationship between vertex local volumes and N± is discrete in nature. In fact, no 
is not even defined for N 4 < Q. Consequently, one cannot assume that the large iV 4 
behaviour is consistent with a curve which intercepts the origin. This fact was accounted 
for by including constants ai and a 2 in the fitted curves. 
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Figure 6.3: (a) Local volume of singular vertices versus vertex coupling k . (b) Local 
volume of corresponding singular link versus k . (The lines serve only to guide the eye.) 
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Singular Link 

Let us now consider the second conjecture made by Catterall et al. which states that the 
local volume of primary singular subsimplices diverge with triangulation volume according 
to a two-thirds power law |f78| . Figure |6.3| (b) plots the local volume of the singular link 



(formed by the product of singular vertices a° sl and a° s2 ) versus k q for each target volume. 
If the link local volume n\ is zero, then this implies that the two singular vertices are not 
connected. From the Monte Carlo data it became apparent that o° sl and a° s2 always form 
a link a] with local volume n{a\). In other words, n{a\) was always non-zero. By fitting 
the kq = results to the following curve 

n(al) = p + qN 4 r , (6.13) 

we find that r = 0.66(1) with x 2 ~ 17 d.o.f. - . This is perfectly consistent with the 
conjectured two-thirds power law divergence. 

It was claimed by Catterall et al. that it is possible to understand the origin of this 
power law divergence in terms of simple geometry [ 78[] . In four dimensions the dual-sphere 



of a vertex is a 3-sphere, and the dual-sphere of a link is the 2-sphere intersection of two 
3-spheres. By assuming that simplices are evenly distributed one may equate the volumes 
of simplicial and smooth spheres. It is then claimed that the two-thirds power is derived 
from the fact that the volume of smooth 2-spheres and 3-spheres scale with the square 
and the cube of the radius respectively J78| . When considering singular structures in 
five dimensions, we indeed encounter three overlapping 3-spheres which intersect at a 2- 
sphere. However, each 3-sphere is the intersection of two 4-spheres. Presumably one could 
apply the same arguments to the 4-spheres. Given that the local volume of the 4-spheres 
diverge linealy with triangulations volume, one would conclude by the same arguments 
that the local volume of the 3-sphere diverges according to a three-quarters power law. It 
seems that this leads to a contradiction since singular links in five dimensions are known 



to diverge linearly with volume (see figure 5 in [78 



Additional Remarks 

So far we have detected the presence of singular links in four dimensional triangulations for 
kq = 0, and reproduced known results. Close scrutiny of Monte Carlo time series plots has 
revealed that, occasionally, the local volume of one of the singular vertices dramatically 
reduces to 'normal levels' (but always remains the second largest local volume). However, 
the local volume of the link never reduces to zero. In other words, the two vertices with 
the largest local volumes are always connected by a link. The most astonishing thing 
about this phenomenon is that the transformation from singular link to singular vertex 
and back again occurs very quickly (typically within about 200 Monte Carlo sweeps). 
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Furthermore, these events are very rare - with an occurrence rate of only about 1 in 10 3 
measurements. This type of phenomenon was apparently first reported by Catterall et 
al. following an investigation into the possible existence of 'volume barriers' in the space 
of triangulations J78[] . It was found that the ensemble 'tunnels' between triangulations 
characterised by a singular link (so-called super-crumpled states), a single singular ver- 
tex or those with no singular vertices. These occurrences were only observed for small 
triangulations N\ = 4k. 

Our time series plots show that the ensemble is generally very stable in the space 
of triangulations with singular links. These rare tunneling events may be a finite size 
effect. However, our results are consistent with a tunneling rate that is independent of 
the volume, at least up to N\ = 40k. This effect has also been observed for kq < 0. In 
fact, the tunneling rate is marginally higher deeper into the crumpled phase. The reason 
behind these events is not currently well-understood. This brings us to our last comment 
of the section. 

These tunneling events may be related to another phenomenon concerning singular 
structures discovered deep in the crumpled phase of four dimensional dynamical trian- 
gulations. From figure |673|(a) it is evident that triangulations deep within the crumpled 
phase have only one singular vertex, and hence no singular link. Figure |6.3| (b) also il- 
lustrates this fact - though less convincingly. Interestingly, this effect was predicted by 
Bialas et al. [BO]. Perhaps these effects deserve closer attention. 



6.3.1 Generalisation of Entropy Arguments 

The non-existence of singular vertices in three dimensions was earlier accounted for in 
terms of entropy. In the simple three dimensional case one only considers the entropy 
contributions of a single subsimplex. In higher dimensions, primary singular subsimplices 
are formed by the product of secondary singular subsimplices. This complication requires 
some modification of our (new) arguments. In this section, we generalise our arguments 
to arbitrary dimension. 

Consider an arbitrary B d submanifold of a <i-dimensional triangulation. Let the space 
of triangulations of this structure be represented by Td-baii- For simplicial (i-manifolds, 
the star of a {d — 3)-subsimplex is homeomorphic to B d , which has an S^ 1 boundary. 
Now, the dual-sphere of a (d — 3)-subsimplex is an S 2 submanifold of the S^ 1 bound- 
ary. Such d-balls are described as (d — 3) -sub simplex- singular because each simplex of 
the (i-ball is common to the (d — 3)-subsimplex. The corresponding space of triangula- 
tions is symbolised by %i n g,d-z- From our definitions of Td-baii and T sing ^-3 it is clear 
that T s ing,d-z C Td-baii- One can define the subspace of non-(d — ?>)- sub simplex- singular 
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triangulations T n on,d-3 such that the following expression is true. 

T S ing,d-3 U ^non,d-3 Q ^d-ball- (6.14) 

One must now consider the entropy of these structures as a function of volume. The 
entropy associated with these spaces of triangulations can be equated to give 

Q\T s i ng> d-z) + ^(7"non,d-3) = Q\Td-ball)- (6.15) 

It was claimed by Catterall et al. that the existence of singular subsimplices can be 
accounted for by considering the entropy of the dual-spheres individually. In this thesis, 
we propose that it makes more sense to consider the entropy associated with the entire 
singular structure, whose neighbourhood is homeomorphic to a (i-ball. One then compares 
the entropy (or growth rates of entropy) of the singular structure <i-ball and an arbitrary 
<i-ball of the same volume. These ideas may be clarified with an example. 

Consider the purely hypothetical case of two singular vertices in three dimensions whose 
product forms a singular link. The vertices have S 2 dual-spheres which intersect at a 1- 
sphere dual to the link. Our previous understanding of the existence of singular structures 
accounted for their existence by quantifying the entropy of the various singular subsim- 
plices individually. Instead, one ought to consider the entropy of the entire singular 
structure - that is to be treated as a single object. Previous thought would have it that 
the entropy associated with the link is effectively that of the 1-sphere. In fact, we propose 
that the entropy fl associated with the structure is given by 

Q = 2Q{S 2 ) -QiS 1 ), (6.16) 

since the entropy contribution of the 1-sphere intersection is overcounted. 

By generalising this argument to higher dimensions, one can show that the entropy 
associated with a primary singular subsimplex in ^-dimensions is given by 

ftCWd-g) = E^ 2 Q+i (-lynts^- 1 ). (6.17) 

8=0 

The existence of the singular structure depends on 

^(%ing,d-3) 3> ^(%ion,d-3) (6.18) 



By substituting (|6.15|) and fl6.17|) in fl6.18j ) and using the fact that Q(Td-baii) = Q(B d ), we 



get the following expression. 



d-3 

2^ d - 2 C !+ i(-l) i fi(5 rf - i - 1 ) >0(5 d ) (6.19) 

i=0 
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Clearly, this inequality applies to d-balls of the same volume Nd- Using the same argu- 
ments explained earlier, one can show that this is true when 

N d {B d ) = Y J d - 2 C i+1 {-l) i N i {S i ), (6.20) 

where N^S 1 ) is the volume of an S l dual-sphere. (One assumes that multiple dual-spheres 
of the same dimension have equal volumes.) In other words, inequality ( |6.19 ) is valid only 
for structures which satisfy ( |6.20|) . 

For example, by substituting d — 3 in ( |6.19| ) we see that singular vertices exist if the 
following inequality is satisfied, cf. inequality (|6.7|). 



2tt(S 2 ) > tt(B 3 ) (6.21) 

In the same way one can show that the existence of singular structures in four dimensions 
would be due to 

AQ{S 3 ) - 2Q{S 2 ) > fi(5 4 ), (6.22) 
providing the volumes of the dual-spheres satisfy 

N A (B A ) = 2N 3 {S 2 ) - N 2 {S 2 ). (6.23) 



Unfortunately, it is very difficult to show that ( |6.22| ) is true. Obviously, in the thermody- 
namic limit one considers the growth rates of entropy and so ( |6.20| ) does not apply. 



6.4 Five Dimensions 

Finally, let us now turn our attention to singular structures in five dimensional triangu- 
lations. These have already been investigated to a certain extent by Catterall et al. [|78|| . 
Here, our objective is to add to their work in the light of the newly discovered phase 
structure of five dimensional dynamical triangulations. When k = k 3 = 0, one expects 
to observe three singular vertices and three singular links which collectively form a sin- 
gular triangle. Our algorithm is therefore designed to measure the three largest vertex 
local volumes of a given triangulation. It then calculates the local volumes of the three 
links formed by pairs of the three vertices. Lastly, it determines the local volume of the 
triangle formed by the product of all three vertices. Monte Carlo simulations were run 
for the K3 = and k,q = limits over a wide range of couplings, with target volumes up 
to iV* = 40k. 
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6.4.1 Vertex Local Volume Distributions 

Let us begin by taking snapshots of the vertex local volume distribution jV(no) at various 
points in C 2 , which are identified by coordinates (k ,k 3 ). These distributions (taken 
from triangulations with = 30k) are neither normalised or averaged over a number of 
triangulations. Their features are more discernible when plotted on a logarithmic scale. 
In each plot a horizontal line indicates the triangulation volume and a vertical line shows 
the mean vertex local volume no, defined as follows. 



= ^ E n(a°) = ^ (6.24) 



The distribution plotted in figure |BT4|(a) was recorded at point (0, 5) in C 2 , which is 
deep within the branched polymer phase. It is immediately evident that the vast ma- 
jority of vertices have a local volume of six. This is the minimum for five dimensional 
triangulations, i.e. uq > 6. The distribution has no ~ 30 and tails off inversely propor- 
tional to no- The largest vertex local volumes are of order 10 2 . Not surprisingly, this is 
qualitatively very similar to the analogous distribution of the four dimensional branched 
polymer phase, schematically drawn in figure |6.1| . The absence of singular vertices in this 
phase is accounted for by the fact that large k 3 (and/or large Kq) forces extra vertices into 
the triangulations. 

As expected, the distribution for (0, 0) is radically different. The main distinction 
being the existence of a small number of vertices with very large local volumes (in excess 
of 10 4 ). These are plainly seen in figure |6]4](b). Analysis of the raw data confirms that 
these points represent three singular vertices. Other differences between the (0, 0) and 
(0, 5) distributions include far fewer vertices with n ~ 10 1 . This feature also mirrors the 
four dimensional distribution (see figure |Q|) . 

Figure |0|(a) corresponds to a triangulation with (—10,0) that is representative of 
the strong coupling limit of /to — > — oo (see figure |4.2| (a)). By comparing this figure with 
the previous two, one can see how the transition from the branched polymer phase to 
the crumpled phase affects the distribution of vertex local volumes. There are two main 
changes that occur. The first is obviously the appearance of three singular vertices, and 
the second is the disappearance of vertices with small local volumes of order 6 < no < 10 1 . 
This results in an overall increase in no, shown by the vertical lines. In fact, using the 



ideas presented in section |5]3| one can show that (in the thermodynamic limit) the mean 
vertex local volume is bounded from below by no > 30, or in general no > d(d + 1). 

Our last distribution snapshot, shown in figure |6.5| (b), was taken for (0,-5). This 



point might be located in a new phase, whose possible existence was bought to light in 



section 4.2. The distribution is quite unlike the others. As with the previous three cases, 
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we find that virtually all vertices have uq = 6. In this phase there are no vertices with very 
large local volumes, though a small number have n Q ~ 10 3 . More surprisingly perhaps, 
there are no vertices with 6 < n < 100. Clearly, this distribution is fundamentally 
different to those of the branched polymer phase and the strong coupling phase of the 
K3 — limit. The origin of this distribution is not currently understood, though the key 
to its understanding may be the realisation that the five dimensional Einstein-Hilbert 
action cannot be expressed in terms of N$ and Nq. In this strong coupling regime, the 
ensemble seeks to minimise N 3 irrespectively of N . 



6.4.2 Primary and Secondary Singular Subsimplices 



The local volumes of secondary singular subsimplices are conjectured to diverge linearly 
with volume [|78| . In five dimensions, this translates to the linear divergence of vertex and 
link local volumes with N$. As in four dimensions, the local volumes of the three vertices 



are only approximately equal. These are denoted by n{a Q sXj 



n(a° s2 ) 



and n(cr° 3 ), such that 



n«) > n(a° s2 ) 



> n(a® 3 ). The link local volumes are given by 71(0-^), n(al 2 ) and n(o~\ 



These structures form a triangle with local volume n(o~ 2 s ) in the arrangement shown in 



figure 5.6 




Figure 6.6: Labelling of primary and secondary singular subsimplices. 



Figures |BT7|(a) and |6[7](b) plot the mean values of n^a^), n(a° s2 ) and n{a° s3 ) versus k 
and k 3 respectively, for each target volume. On each graph, a key labels the vertex local 
volumes of each target volume in a logical manner. Figures |6.8| (a) and |6.8| (b) plot the 
mean values of n^a^), n{a\ 2 ) and n(al 3 ) versus k and k 3 respectively, for each target 
volume. Again, a key labels the link local volumes in each graph. The k,q = K3 = data 
was fitted to the following sets of curves. 

n(a° sl ) = ai + hN^ n( ( x s 1 1 ) = f 1 + g x N^ (6.25) 

n(a° s2 ) = a 2 + b 2 N 5 c > n(a l s2 ) = f 2 + g 2 N 5 h * (6.26) 

n(o-%) = a 3 + b 3 N B c ° n(al 3 ) = f 3 + g s N 5 h3 (6.27) 
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Figure 6.8: Scaling plots of the three largest link local volumes versus (a) the vertex 
coupling k and (b) the tetrahedron coupling k 3 . (Lines serve only to guide the eye.) 
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The exponent coefficients of the best fit curves are given in table |6.1| along with an 
indication of the 'goodness-of-fit' in terms of \ 2 per degree of freedom. As expected, 
these results are consistent with linear divergence. 
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Table 6.1: Best fit exponent coefficients for secondary singular subsimplices. 



Let us now consider the primary singular subsimplex, which in five dimensions is a 
triangle. The algorithm measures the local volume of the product of the three singular 
vertices, whether they form a triangle or not. Figures |6.9| (a) and |6.9| (b) plot the mean 
value of n(a 2 ) versus kq and K3 respectively, for each target volume. In the same manner, 
we fit the Ko = K3 = data to the following curve. 

n(a 2 ) =p + qN 5 r (6.28) 

The best fit curve was found to have r = 0.63(2) with x 2 ~ 3.2 d.o.f. -1 . These results are 
reasonably consistent with the conjectured two-thirds power law divergence. 

New Results 

Our study has shown that the vertex local volume distribution Af(no) of five dimensional 
triangulations is independent of volume. There is a direct correspondence between the 
branched polymer phases of four and five dimensional dynamical triangulations. This is 
not so surprising. For large kq and/or k$, the action promotes vertex insertion moves. 
Consequently, the number of vertices of the triangulations grows which in effect bars the 
existence of singular structures. 

One can draw parallels between singular structures in four dimensions and the ^3 = 
limit of five dimensions. It was interesting to find that triangulations in the crumpled 
phase are not always characterised by singular triangles. Monte Carlo time series plots 
of 

^(^si)) n (° r s2) an d M "^) f° r K o — showed that occasionally the ensemble tunnels 
between states with singular triangles and singular links. Like the four dimensional case, 
the tunneling events occured very quickly (within 200 sweeps). However, in five dimen- 
sions the tunneling rate was of order 1 in 100. It was found that the three largest vertex 
local volumes always form a triangle. 
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Figure 6.9: Scaling plots of the largest triangle local volume versus (a) the vertex coupling 
k and (b) the tetrahedron coupling k 3 . (Lines serve only to guide the eye.) 
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Deeper into the crumpled phase ensembles were found to thermalise towards trian- 
gulations with singular links for n <C 0. Analogous effects have been observed in four 
dimensions and may represent metastable states in the space of triangulations [fT8[. A 



simple experiment was devised in order to test this was the case in five dimensions. The 
ensembles were originally thermalised from cold configurations, i.e. branched polymers. 
The simulations were repeated by thermalising from the k,q = triangulations. It was 
found that the singular triangles remained intact and stable for all kq < 0. Despite this, 
the ensemble continued to tunnel between singular triangle and singular link states in the 
same manner. These are the results which are shown in figures |0](a), |6.8| (a) and |6.9| (a). 
In conclusion, this evidence suggests that triangulations with singular links (in five di- 
mensions) represent metastable configurations. 

Close examination of the results reveals that 72(0"^), n(a Q s2 ) and n((T° 3 ) increase slightly 
for hq < 0. This is presumably attributable to the action which tends to reduce N . It 
appears that the local volumes of secondary singular subsimplices diverge linearly with 
volume throughout the crumpled phase. It is difficult to confirm the two-thirds power 
law scaling in this limit, since we have data for only three target volumes. Interestingly, 
the singular triangle local volume decreases (possibly to an asymptote) as k — > —00. 
Triangulations with k 3 < in the k,q = limit have no singular vertices. This fact was 
discovered earlier in section |6.4.1| . Unfortunately, it is not clear from our limited study 
how (or whether) the vertex local volume distribution depends on triangulation volume. 



6.4.3 Entropy 

As an explicit example of our generalised entropy arguments, we shall now apply these 
ideas to singular structures in five dimensional triangulations. Consider an arbitrary B 5 
submanifold of a five dimensional triangulation. The boundary of a 5-ball is a 4-sphere. 
Let the space of triangulations of this structure be represented by T^-baii- For simplicial 
5-manifolds, the star of a triangle is homeomorphic to a 5-ball, and its dual-sphere is a 
2-sphere. Clearly, the dual-sphere is a submanifold of the 4-sphere boundary of the star. 



This is illustrated schematically in figure |6.10| . Such 5-balls are described as triangle- 
singular. The corresponding space of triangulations is symbolised by T sing ^- Clearly, 
from our definitions of T 5 - ba u and T sing 2 we find that T S ing,2 C T^-baii- One can then define 
the subset of non-triangle- singular triangulations T non ^ such that 

%ing,2 U T non ^ C %-hall- (6.29) 



The entropy of an arbitrary 5-ball may be represented by Q(B 5 ). By substituting d = 5 
into equation ( |6.17|) one can show that the entropy associated with the singular structure 
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Figure 6.10: Schematic diagram of a triangle-singular B 5 with S A boundary. The S 2 
dual-sphere of the triangle in the B 5 is a submanifold of the S 4 boundary. 



as a whole is given by 



sing, 2 ) 



3tt(S 4 ) -3tt(S 3 ) + tt(S 2 ). 



The entropy of these spaces of triangulations may be equated as follows. 

-ball I 



(6.30) 



(6.31) 



Therefore, the existence of the observed singular structures in five dimensions would be 
due to 

^(T sing>2 ) >fi(T„ onj2 ). (6.32) 

By substituting equations (|6.30| ) and ( |6.31| ) into inequality ( |6.32| ) and using the fact that 
Q(T 5 - ba ii) = Q(B 5 ) one can derive the following expression. 



6Q(S 4 ) - 6Q(S 3 ) + 2Q(S 2 ) > Q(B Z 



(6.33) 



Again, such an expression is very difficult, if not impossible, to verify. One cannot even 
qualitatively substantiate this relation. Obviously, we are interested in comparing the 
entropy of two 5-balls of the same volume. Now, the triangle-singular 5-ball is comprised 
of three 4-spheres, three 3-spheres and one 2-sphere. Let their volumes be given by N^B 5 ), 
N^S 4 ), N 3 (S 3 ) and N 2 (S 2 ) respectively. Using ( |6.20|) one can show that inequality ( |6.33| ) 
is only valid when the volumes satisfy the following equation (assuming that the 4-spheres 
and 3-spheres have equal volume). 



N 5 (B 5 ) = 3N 4 (S 4 ) - 3N 3 (S 3 ) + N 2 (S 2 



(6.34) 



6.4.4 New Phase Structure 



In this chapter, it has become clear that the distinction between the two strong coupling 
regimes (first brought to light in chapter |j) is directly related to the phenomenon of 
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singular structures. Our study has shown that the crumpled phase of the K3 = limit is 
characterised by the existence of singular triangles. This phase is known to include the 
point (0, 0) in C 2 . One may then trivially deduce that the k,q = limit also corresponds to 
the same phase at k 3 = and perhaps its vicinity. Having observed vertex local volume 
distributions for k 3 < that are fundamentally distinct from those of k 3 = 0, one may 
wonder whether they represent different phases in C%. The main purpose of this section 
is to explore this possibility. 

Recall that evidence of a phase transition near k 3 = —1.0 was discovered in section [4.2| . 
It seems very likely that this point in C 2 separates two phases highlighted clearly in 
figures |6.71 (b) and |6.8| (b). To avoid confusion, the location of the new phase transtion is 
denoted by k.^. For clarity and convenience, the phases are named as follows: the crumpled 
phase is located in the region k 3 < k 3 < k 3 and the compact^ phase is found at k 3 < k^. 

Our specific goal is therefore to confirm that a phase transition does separate these 
phases. This is done by measuring the tetrahedron specific heat X3( K 3) an d the local vertex 
distribution for higher resolutions of k 3 . These simulations were run for triangulations 
with target volume = 10k. Figures |6.11| (a) and |6.11| (b) (on page |107|) plot the three 



largest vertex local volumes and %3 versus k 3 respectively. At k 3 = 0, triangulations 
are characterised by three singular vertices. As k 3 is increased, the singular vertices 
are progressively suppressed. By comparing both figures, their disappearance is obviously 
marked by the phase transition k\ ~ 0.454(1). As k 3 is decreased we find that the singular 
vertices vanish near k 3 w 0.951. Again, this correlates with a small peak in x 3 . It would 
appear that the transition is strongly first order. 

The new phase transition was traced for kq > and for kq < 0. Again, this was done 
by fixing k and tuning k 3 to the transition. The transition was traced from k = — 2 
to Ko = 1 (see figure |6.12j on page |108|). The known phase are labelled accordingly. The 



question marks indicate regions in C 2 which are yet to be explored. These results confirm 
the existence of a second line of phase transitions. Unfortunately, it was not possible to 
reconnoitre C 2 any further due to time constraints. 



4 This name reflects the fact that triangulations of the compact phase have smaller mean geodesic 
distances than those of the crumpled phase, as seen in figure 4.5(b). 
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Figure 6.11: Vertex local volume versus the tetrahedron coupling k 3 for k = 0. (b) 
Tetrahedron specific heat X3 versus k 3 . 
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Chapter 7 
Summary 



LTHOUGH dynamical triangulations has been actively researched for over 
ten years, physicists are yet to fully understand its possible links with a 
consistent theory of quantum gravity in four dimensions. It appears that 
the naive generalisation of two dimensional dynamical triangulations does 
not yield a physical continuum limit. Despite this fact, research into dynamical trian- 
gulations for d > 2 has continued, albeit in new directions. One such direction relates to 
the implementation of modified actions. This is usually done either by including matter 
fields or by adding measure terms. The net effect is the addition of an extra coupling 
constant. One then hopes that a physical continuum limit may be reached at a point in 
the expanded space of couplings. 

In this thesis, the space of couplings is expanded by considering the most general action 
that is linear in components of the /-vector. These actions happen to be equivalent to 
the discretised Einstein-Hilbert actions for d < 5. The purpose of this work is to discover 
new aspects of dynamical triangulations by viewing the model from a fresh slant. This 
approach will teach us more about dynamical triangulations in general. Let us now 
summarise the results of the preceding five chapters. 

This thesis is naturally divided into two parts. The first part aims to give credibility to 
the second. A fundamental requirement of Monte Carlo simulations is that local moves 
are ergodic in the space of configurations. In chapter the five dimensional (k, I) moves 
were proven to be ergodic in the space of combinatorially equivalent triangulations. Un- 
fortunately, the negativity of the Hauptvermutung in five dimensions means that the proof 
cannot be extended (in general) to the space of homeomorphic simplicial 5-manifolds. It 
is not known (to the author) whether the Hauptvermutung is true for the 5-sphere. The 
physical implications of its negativity are not well-understood. 
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Beyond this formal proof, one must contend with the issues of computational ergodicity 
and the recognisability of manifolds. In the context of triangulations, the unrecognisabil- 
ity of simplicial manifolds T and T' implies that the number of (k, I) moves needed to 
transform T into T" is not bounded by a recursive function of their volume. Of course, 
even recursive functions can grow extremely quickly. Therefore, it is argued that, on a 
practical level, it is more relevant to consider whether the number of (k, I) moves required 
to transform T into T' is computationally manageable. This relates to whether a series of 
(k, I) moves can be performed by modern computers in a reasonable time scale (i.e. days 
or weeks). 

Chapter ^| aimed to show that the canonical partition function of five dimensional 
dynamical triangulations is exponentially bounded from above as a function of volume. 
Of course, numerical evidence can never constitute a proof. No attempt was made to 
construct an analytic proof - mainly because of the obvious complexity of the problem. 
Unfortunately, the outcome was inconclusive since the results were found to be equally 
compatible with exponential and factorial divergence. The chapter was concluded by high- 
lighting the problems associated with the numerical approach, and suggesting a possible 
solution. 

The second part of the thesis deals with the more interesting aspects of dynamical trian- 
gulations, at least from a physicist's point of view. In chapter |], Monte Carlo simulations 
reveal the existence of a phase boundary in C2, the two dimensional space of coupling 
constants. The weak coupling phase is characterised by elongated triangulations of high 
scalar curvature and vertex density. The triangulations of the strong coupling regimes 
are generally crumpled with low scalar curvature and vertex density. The phase transi- 
tions were identified by sharp peaks in the specific heats. The results also hinted at the 
existence of a third phase in the strong coupling regime of the kq = limit. 

Traditional finite size scaling methods failed to determine the order of the phase tran- 
sition at points (kq, 0) and (0, Kg). This was probably due to a combination of insufficient 
statistics and strong finite size effects. However, time series plots showed some evidence of 
bistability, which may be indicative of first order phase transitions: the signal was some- 
what clearer for (0, k%). Further simulations at large volumes are required for verification. 
The phase boundary was traced over a large range of couplings. It forms a straight line 
for small kq and k 3 , but is slightly curved in regions away from the origin. The reason for 
this is unknown. Interestingly, it was found that the phase boundary appears to be more 
strongly first order as «o ~* °o- 

In chapter the weak coupling phase triangulations were found to have du ~ 2, 
irrespective of volume. This discovery confirmed their branched polymer nature. In 
fact, it was found that the partition function is saturated by (branched polymer) stacked 
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spheres for large k,q and/or K3. Finally, a generalisation of Walkup's theorem to arbitrary 
dimension was proved from elementary arguments and simple graph theory concepts. This 
was achieved by constructing dual graphs of whole triangulations. The real result of this 
work is not so much the generalised inequality itself, but the methods used to prove it. 

The thesis was brought to a close with a broad study of singular structures in chap- 
ter |6|. Their existence is accounted for in terms of new entropy arguments, based on those 
proposed by Catterall et al. [78]. They appear to explain the non-existence of singular 
vertices in three dimensions. Monte Carlo simulations were run for four dimensions over 
a range of volumes and couplings. Time series plots revealed that the singular link occa- 
sionally transforms into a singular vertex and back again very quickly (within about 200 
sweeps). This phenomenon was observed for various values of kq in the crumpled phase, 
including k = 0. It was found that triangulations deep within the crumpled phase are 
characterised by a single singular vertex. This confirmed the prediction made by Bialas 



et al. g. 

The very last section dealt with singular structures in five dimensions. It was found 
that singular structures play an important role in the K3 = limit for k < k . Perhaps 
the most interesting outcome of this research was the discovery of a new 'compact' phase 
in the strong coupling regime of the k,q = limit, where the action is the discretised 
Einstein-Hilbert action. It was found that triangulations endowed with singular structures 
dominate the partition function in the range K3 < K3 < where k\ is the location of 
the newly discovered phase transition. Its existence was hinted at by measurements of 
c 2 {k3,N 5 ) in section fL2| . The compact phase was found to be characterised by a new 
distribution of vertex local volumes, which is shown in figure |6.5| (b). Finally, the account 
ended with an in-depth study of the new phase transition, which appears to be strongly 
first order. The phase transition was then traced in C2 for k,q ^ 0. 



Let us now look at how these results relate to dynamical triangulations in general. Our 
study of stacked spheres in chapter || allows us to make the following generalisation. 
Consider <i-dimensional dynamical triangulations, in which the action is linear in any 
components of the /-vector - one of which must be Nd- One can always identify a 
'weak coupling phase' which is characterised by branched polymer triangulations. Let us 
substantiate this claim. One can generalise the ideas of section |5.4j to show that branched 
polymer stacked spheres of a given volume Nd correspond to the upper kinematic bound 
of Ni, for i = 0, . . . ,d — 1. In other words, regardless of the form of the action, one can 
always set the couplings either to +00 or —00 such that the partition function is saturated 
by branched polymer stacked spheres. Using equation (|2.20|) , it is quite simple to prove 
the following upper kinematic bound, which holds in the thermodynamic limit for d > 1 
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and i = 0, . . . , d. 

Ni (d+l)(d 



- (7.1) 
N d i\(d-i + l)\ { } 

The substitution of d = 4 and i = 2 yields i? — > 2.5 as expected. Of course, when i — 
we are left with the upper kinematic bound of the vertex density. 

* - i (7.2, 

Singular structures are known to play an important role in dynamical triangulations, 
at least for d > 3. Based on our knowledge of dynamical triangulations in four and five 
dimensions, one may speculate that a phase transition separates the weak coupling limit 
and zero coupling in all cases where the action is linear in components of the /-vector. 
For an action with n independent couplings, one may expect a (n — l)-dimensional surface 
of phase transitions to exist in the space of couplings. This is exactly what we have found 
in five dimensions. 

One can make further speculations as a result of this work. Naively, dynamical trian- 
gulations could give us 'quantum gravity' in any dimension we like. If this was the case, 
then one might expect the phase structure of four dimensions to be similar to that of the 
Ko = limit of five dimensions, assuming that the Einstein-Hilbert action is special in 
some sense. Our research has shown that this is not the case. In fact, the phase structure 
of the K3 = limit is more similar to the four dimensional case - in terms of singular 
structures. From this, one may speculate that the four dimensional phase structure is as 
it is because the action is proportional to N d and N , not because it is the discretised 
Einstein-Hilbert action. 



Finally, let us end with some suggestions for future research. First of all, as mentioned 
earlier, the negativity of the Hauptvermutung implies that the (k, I) moves may not be 
ergodic in the space of homeomorphic simplicial manifolds. It is not clear whether any 
unreachable triangulations are physically significant; that is, whether they should be 
included in the measure. This matter seems to be quite important and ought to be 
resolved. 

The existence of an exponential bound remains an outstanding problem in five dimen- 
sions. In the absence of an analytic proof, one must rely on numerical evidence. This 
issue might be 'resolved' simply by extending the investigation of chapter § to larger 
triangulations. However, there is no guarantee that the outcome would be conclusive. 

Although much of C2 is now mapped, there are a number of issues which are yet 
to be addressed. Perhaps the first task would be to establish the order of the phase 
transitions at (kq, 0) and (0, k^) with certainty. This may be achieved by seeking evidence 
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of metastability, particularly over a range of volumes. One may then consider whether 
the whole boundary is of the same order. Other prospects include the possible existence 
of critical points. Perhaps the best course of action is to further trace the second phase 
boundary, which separates the crumpled and compact phases. The identification of a 
critical point in C 2 is a real possibility, and ought to be pursued. If contact can be made 
with the continuum, then one must consider whether the theory has anything to do with 
'gravity'. The answer would be very interesting, regardless of the outcome. 

Our generalisation of Walkup's theorem demonstrates the power of simple graph theory 
concepts when applied to the dual construction of triangulations. Perhaps, this connec- 
tion between dynamical triangulations and regular non-planar graphs can be exploited 
further. It would also be interesting to see if (or how) our results relate to recent ana- 
lytic advances jS^l . Lastly, it would be nice to explain the nature of the vertex local 
volume distribution of the compact phase. This may help us understand more about the 
transition between the crumpled and compact phases and may reveal connections with 
other work. Perhaps this could be achieved in terms of the ' Balls- in-Boxes' model ||S3" |. 
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Computer Programs 



A.l Dynamical Triangulations Computer Program 



N this section we present a condensed overview of the main computer program 
used in our research of dynamical triangulations. In particular, we outline the 
general algorithm and mechanics of the Monte Carlo program, which was written 
in C by Simon Catterall for arbitrary dimension [|B8[1 . This program, which we 



refer to as main, was written in C primarily to utilise its dynamic memory allocation 
capability^. It also allowed for the effective use of pointers. These features are by no 
means mandatory as equivalent programs have been written in FORTRAN ||85|| . One of 
the most difficult problems in creating such a program is finding a good way of abstractly 
representing the simplices and their connectivity. Let us now outline the framework of 
main with this issue in mind. 

Each simplex of a <i-dimensional simplicial manifold may be represented by d + 1 
labels for each of its vertices and some form of identification of its d + 1 neighbouring 
simplices. This knowledge is sufficient to reconstruct any triangulation. In our program, 
each simplex is abstractly represented by a set of variables known as a structure, which 
has an associated pointer. These pointers are grouped together into an array called 
simplex_point. Each structure contains (amongst others) two main variables: an array 
of d + 1 integers which label the vertices of the corresponding simplex and an array of 
d + 1 pointers which point to the d + 1 neighbouring structures (simplices). In main 
these arrays are called vertex and neighbour respectively. These constructs enable the 



1 From here on, any technical C terms will be printed in italics. Programs, functions and variables will 
be printed in lower case typed text and global constants are shown in UPPER CASE TYPED TEXT. For 
full details of the C programming language the reader is referred to a popular standard text by Kelley 
and Pohl [H. 
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program to simulate triangulations of any dimension. This use of pointers allows us to 
'connect' neighbouring simplices in a very intuitive manner. 

A. 1.1 Update Algorithm 

Recall that a (k, I) move in ci-dimensions may be regarded as the replacement of an i- 



(sub)simplex with its dual (d — i)-( sub) simplex (see section 2.3.1 ). In this context we have 
% = d + 1 — k and hence d — i = d + 1 — I because k + I = d + 2. This realisation of the 
(k,l) moves allows for their definition for arbitrary dimension since a move is uniquely 
defined for a given i and d, where i — 0, . . . , d. Here we describe the four main stages in 
the update procedure as they occur in order. 

• In order to maintain ergodicity it is imperative that the moves are performed on 
random (sub) simplices. This is done by first selecting, at random, a pointer from 
simplex_point which corresponds to a simplex, followed by a move of type sub, 
where sub = i = d + 1 — k. The (sub)simplex to be transformed is chosen by 
randomly selecting d + 2 — k of the d + 1 vertices of the simplex. 

• Before executing the transformation the algorithm must first check that the chosen 
(sub)simplex is 'legal'. An i- (sub) simplex is considered legal if its local volume is 
d — i + 1. For example, one cannot perform a (3, 1) move on a vertex if its local 
volume is four. If the (sub) simplex is illegal then the algorithm restarts. Notice 
that this check is unnecessary if i = d since vertex insertion moves are always legal. 

• Once a legal (sub) simplex has been identified, the corresponding transformation will 
be accepted by the Metropolis algorithm with a given probability - determined by 
the change in the action AS. If AS is zero or negative then the move is accepted. 
If AS is positive then the probability of acceptance is exp(— AS). The original 
Metropolis algorithm of main calculated AS[T,Kq]. Some minor alterations were 
made in order to study five dimensional dynamical triangulations with our new 
action Sf[T, k , k 3 ]. 

• If the move is accepted by Metropolis then the algorithm checks that the move 
is geometrically 'allowed'. This condition insists that the star of each vertex is 
homeomorphic to a (i-ball and hence the triangulation retains its manifold structure. 
Transformations that do not satisfy the condition are therefore discarded. The 



importance of this point was stressed by Schleich and Witt [S6[ . 



If a candidate move passes through each step then the transformation will be executed. 
Having explained the update process of main, let us now describe each phase of the Monte 
Carlo simulation life cycle. 
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A. 1.2 Growth 

The C program written by Catterall is capable of generating simplicial manifolds of spher- 
ical topology for arbitrary dimension Jj3]. This is made possible by realising that such 
structures can be grown from minimal d-spheres using (k, I) moves. The initial configura- 
tion is equivalent to the surface of a (d+ l)-simplex and is generated by a function named 
initial_conf ig which is called by main. This can then be grown to any chosen target 
volume N l d within technological bounds. This is done most efficiently by only perform- 
ing vertex insertion moves that bypass the Metropolis algorithm. This process generates 
stacked ci-spheres (see section |5.3| ). 

Clearly, the triangulation target volume is limited by RAM (random access memory) f\ 
However, the numerical performance of Monte Carlo simulations is limited more by CPU 



power (see section |A.2|) . Once the required target volume has been attained, the algorithm 



then passes into the thermalisation stage. 



A. 1.3 Thermalisation 

In general, the initial fully grown triangulation will not be in thermal equilibrium. It 
is known from our investigations in chapter ^] that stacked spheres saturate the weak 
coupling limit. If order to study the model at stronger coupling then one must allow a 
finite Monte Carlo time for the initial triangulations to thermlise. In these simulations 
Monte Carlo time is measured in units of sweeps. A sweep is defined as N\ legal attempted 
moves. (Alternatively, one could equally define a sweep as N l d completed updates.) The 
number of thermalisation sweeps THERM is defined as a global constant. 

During this stage the simplex coupling is tuned at regular intervals of TUNE_COUPLING 
sweeps using the following iterative procedure 

= + 27(i\fr- A®, (A.l) 

where is the updated simplex coupling and 7 is related to the volume fluctuation 
parameter DV by 

DV = 7"i (A.2) 

It has been verified that thermal averages are independent of 7 for small 7, but have 
errors that grow as 7 — > ||68|| . In our simulations 7 was chosen to be 7 ~ 0.0005, 
corresponding to DV ~ 32. It is important that DV is not set too low, as one may lose 

2 According to Catterall, the memory required for a d-dimensional triangulation of volume is of 
the order of 8(d + 6)N c i bytes |6q] . For example, a five dimensional triangulation of volume N 5 = 10k 
requires approximately 0.9MB of memory, which is a modest amount for current technology. 
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ergodicity. Ergodicity demands that all of the d + 1 (k, I) moves are implemented, which 
means that DV must be greater than d. 

At intervals of CHECKPOINT sweeps an output file called dump is created. It contains 
all the necessary information needed to regenerate a triangulation of a given volume and 
coupling. One can set main to read-in a dump file so that triangulations do not have to be 
grown and thermalised for each simulation. Once the simplex coupling is tuned and the 
ensemble is fully thermalised main then passes into the measurement stage. 



A. 1.4 Measurement 

At this point the main begins a series of SWEEPS sweeps during which measurements of 
certain physical observables are made at regular intervals of GAP sweeps. These include 
the triangulation volume Nd, the total number of vertices Nq and the average geodesic 



distance g (see section |4.2.2|) . The number of origin simplices is controlled by HITS. All 
measurements are sent to an output file called data. As with the thermalisation stage, 
the triangulation connectivity is sent to dump at regular intervals of CHECKPOINT sweeps. 

It was necessary to measure A3 in order to calculate the scalar curvature R and 
tetrahedron specific heat Xs- This minor problem was solved by introducing a mechanism 
of monitoring the /-vector. A d+1 dimensional array of integers called Tot_up was created, 
whose elements correspond to the components of the /-vector. Prior to the growing stage 
the elements of the array were set to those of a minimal 5-sphere. 

Tot_up[i] = d+2 C l+1 = 7 C m (A.3) 

Then for each (k, I) move that was accepted, the elements of Tot_up were updated ac- 
cordingly using equation ( |2.20|) . Table [A. 1| shows AiVj for each of the six (k, I) moves. On 



completion of SWEEPS sweeps, a final dump file is created. 



(M) 


sub 


AiV 


AAi 


AiV 2 


AiV 3 


AiV 4 


AiV 5 


(1,6) 


5 


+1 


+6 


+15 


+20 


+ 15 


+5 


(2,5) 


4 





+1 


+5 


+10 


+9 


+3 


(3,4) 


3 








+1 


+3 


+3 


+1 


(4,3) 


2 








-1 


-3 


-3 


-1 


(5,2) 


1 





-1 


-5 


-10 


-9 


-3 


(6,1) 





-1 


-6 


-15 


-20 


-15 


-5 



Table A.l: AAj for the five dimensional (k,l) moves. 
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A. 2 Computer Performance Tests 

In this section we give an indication of the scale of computational resources and time 
required for this type of work. During the course of this research Monte Carlo simulations 
were run on a variety of computers, which included: 

• DEC Alpha 3000/133 MHz and 500/333 MHz workstations (Digital UNIX). 

• PCs: Pentium 200 MHz, 233 MHz and Pentium II 333 MHz (Linux). 

• SGI Origin2000 supercomputer - non-parallel (IRIX). 

Tests were run on each computer to quantify their power when applied to our particular 
needs. The test involve measuring the CPU time needed to complete 10 3 sweeps, with 
K 3 — K o — over a range of target volumes^. Figure |A.1| shows the best fit curves of the 
resulting data. This is a plot of time per sweep (in seconds) versus triangulation volume 
N5. A key identifies each curve with its corresponding computer. 



<D 

E 

CL 
<D 
CD 
g 



DEC Alpha 3000/1 33 MHz 
Pentium 200 MHz 
SGI Origin2000 
DEC Alpha 500/333 MHz 
Pentium II 333 MHz 




5000 10000 15000 20000 25000 30000 35000 40000 45000 50000 

volume 

Figure A.l: Sweep time (in seconds) versus triangulation volume N 5 . 



3 In each case the code was compiled using cc -03 . . . . 
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It is clear that DEC Alpha 3000 machines are considerably slower than the other machines 
tested - in fact, Pentium 200 MHz PCs are over twice as fast. The newer DEC Alpha 
500/333 MHz and Pentium II 333 MHz are faster still. Interestingly, we see that the single 
processor Origin2000 is the fastest for small volumes, but slows very rapidly for larger 
volumes. It seems that this system is best suited for parallelised code. From these results 
we see that, for example, a Pentium 200 MHz PC can complete 10 6 sweeps in about 100 
hours. In conclusion, the most efficient means (in terms of cost as well as time) of running 
these simulations is by using a set of the most powerful PCs available (currently Pentium 
III processors). 

A. 3 Additional Code 

Some additional code was incorporated into main in order to study spatial distributions, 
and singular structures in various dimensions. These function were written by the author. 

Spatial Distributions 

The spatial distribution of simplices n(r) was evaluated by a C function called spatial. 
It was called by main at regular intervals of SPATIAL_GAP sweeps. Memory is allocated 
in main for a two dimensional array of floating point numbers called simplex_data. The 
first dimension corresponds to the geodesic distance r. Its size is set to the maximum r 
expected, typically 10 2 . The second dimension relates to the number of measurements of 
n(r) taken during the Monte Carlo simulation. Its size will therefore be 

SWEEPS .. A . 

• ( A - 4 ) 

SPATIAL_GAP v ' 

Each time spatial is called by main it allocates memory for an array of integers 
spatial_tally whose elements are initialised to zero. The r th element counts the number 
of simplices separated by a geodesic distance of r. Then for each element of simplex_point 
acting as an origin simplex O s , spatial calls geodesic which returns the geodesic distance 
between O s and every other simplex. By using every simplex as O s one obtains an 
'exact' measurement of n(r) rather than an estimate^. Once this loop is completed, 
spatial_tally is then normalised and passed back to simplex_data in main. As the 
volume is not fixed, spatial_tally was normalised for each configuration then averaged 
over configurations, rather than vice versa. 

When the measurement stage is finished, main calls mean_spatial which calculates 
the mean spatial distribution {n(r)) and estimates the errors a((n(r))) using the Jackknife 

4 This can become costly (in computational terms) for large volumes. Fortunately, it is not a problem 
for the scale of volumes studied here. 
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method (see appendix [B.3| ). This is the reason why the distributions must be stored in 
simplex_data. Clearly, such storage of data can become a problem for intensive simula- 
tions. This algorithm is not very efficient. 

A separate C program used measurements of (n(r)) and a((n(r))) to calculate the nor- 
malised volume (V(r)) contained within a geodesic radius r of an arbitrary origin simplex. 
The error in (V(r)) was calculated (rather crudely) by defining the upper and lower limits 
as follows. 

(V(r)) + = E«n(i))+a((n(<)») (A.5) 

i=0 

(V(r))- = E«n(<)>-rT«n(<)») (A.6) 

i=0 

This method may not give accurate estimates of uncertainty, but does have the virtue of 
being conservative. The mean fractal dimension function (dn{r)) was calculated using 
equation (|5.5|) and the uncertainty in this quantity was estimated using equations ( |A.7f ) 
and ( |Op . 

{Mr)) ~ ln(r + 1) - ln(r) (A ' 7) 

Id (r))~ - HV(r + l))--HV(r))- . 
{Mr)) ~ ln(r + 1) - ln(r) (A ' 8) 

Singular Structure 

A new function of main named vertex_order was used to study singular structures for 
dimensions d < 6. It is called at regular intervals of SPATIAL_GAP sweeps. vertex_order 
measures the local volume of each vertex by exploiting the fact that each structure stores 
an array of d + 1 integers called vertex, which label the vertices of the corresponding 
simplex. The algorithm may be divided into the following stages. 

• Memory is allocated for an array of integers called vertex_tally, whose elements 
are vertex integer labels. A loop is then run over all simplices and vertex_tally 
counts the vertex labels. This process reveals the number of simplices common to 
a given vertex, i.e. its local volume. A check then confirms that the number of 
non-zero elements of vertex_tally equals N Q . 

• A second array of integers called singular is then generated which counts vertices 
according to their local volume. For example, the tenth element gives the number 
of vertices with a local volume of 10. The elements of singular are then sent to an 
output file. 
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• From chapter |6], we know that the primary singular subsimplex is the product of d— 2 
vertices for d > 3. The function then identifies the labels of the d — 2 largest vertex 
local volumes and calculates the local volume of the corresponding subsimplex. 

• Given that a primary singular subsimplex has d — 2 vertices, the function then 
calculates the local volume of each of the d ~ 2 Cj + i secondary singular i-subsimplices 
for i — 1, . . . , d — 3. vertex_order works out every combination of subsimplices and 
finds the local volume for each case. 

Once the measurement stage is completed, the output files are then analysed. Specifically, 
a separate program calculates the mean local volume of each structure and the standard 
error in the mean (see appendix [B]). 



Appendix B 



Error Analysis and Curve Fitting 



B.l Independent Samples 



FFECTIVE error analysis is a vital part of Monte Carlo simulation methods. 
It allows us to establish the uncertainty associated with a given measurement. 
Let us first consider the simple case where we have n independent (random) 



Monte Carlo simulations, independent samples are those that are not correlated. Given 
a set of n measurements of X, the mean X is, of course, simply defined as 

X= l -±X t . (B.l) 

"' i=i 

(Since we are dealing with finite systems, it is important to remember that X is an 
estimate of the expectation value of X.) The standard deviation a and standard error in 
the mean s m of the sample are given by 

a 2 = -jZ{X % -Xf (B.2) 

and i i 

s m = . a ps —= a for large n. (B.3) 

One may envisage a scenario in which measurements are not independent, i.e. cor- 
related. Under these circumstances, one cannot obtain reliable estimates of errors using 
equations ( B.2|) and (|B.3|) . Such problems may be avoided by allowing the Monte Carlo 
ensemble sufficient 'time' to decorrelate between successive measurements. This charac- 
teristic time is known as the autocorrelation time r and can vary for different observables. 
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Unfortunately, calculation of r is a non-trivial matter and can be a clumsy process. In 
this thesis, the problem of correlated data was dealt with in another way - by using a 
technique known as binning. 



B.2 Binning 



Suppose that our n measurements of X are correlated. Binning involves dividing the n 
data sets into b 'bins' of equal size and calculating the average of each. If the size of each 
bin is greater than the autocorrelation time, then the averages can be considered as being 
independent. Therefore, in essence, binning is a way of generating independent data from 
correlated data P7| . 

The mean Yj of bin i is given by 

^=-E^%-i)+i fori = 1,..., b. (B.4) 

n j=i 

One can then calculate Y, which is equal to X. 

F = (B.5) 



8=1 

1 



6 J 



„ EE (B.6) 
X (B.7) 



By treating the set of 'measurements' Y{ as independent data, one can define the standard 
deviation and standard error in the mean s m as 



a 2 ^W{Yi-Yf (B.8) 
i=i 

and | | 

s m = ; a ps — ^ cr for large 6. (B.9) 

Clearly, when b = n, binning has no effect and we expect inaccurate estimates of the 
errors in X. As b is decreased, F 4 become progressively decorrelated. Consequently, the 
error in Y (and hence X) increases and thus becomes more accurate. When b ps r, Yi 
are effectively independent and the resulting error estimates are accurate. This method 
is effective providing the sample size is much greater than the autocorrelation time. 
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B.3 Jackknife 



As the name suggests, the jackknife is a rather crude (yet reliable) general purpose method 
of estimating errors |T7|]. It is most commonly applied to samples with unknown distri- 
butions. The technique involves dividing a sample of n independent measurements of X 
into, say, c blocks. One then calculates the c averages Yi of the sample by leaving out 
each block in turn. 



Yi 



c — 



Y E *i+E 



X, 



for i — 1, 



The mean of the jackknife estimators Y is given by 

1 



Y 



-E — 



E E ^ 



i=l 



7 = — 2 



X, 



and the standard error in the mean s m is 



\ 



c —±±{ Yi -Yy. 



(B.10) 

(B.H) 

(B.12) 
(B.13) 

(B.14) 



i=i 



B.4 Curve Fitting 



The testing of observations against theoretical predictions is a fundamental aspect of scien- 
tific research. Here we are concerned with 'goodness-of-fit' hypothesis testing - specifically 
chi-squared tests |p8 |. For example, given a set of observed data, chi-squared tests allows 
us to quantify how well the data fits a hypothesised theoretical prediction. To be more 
precise, chi-squared tests tell us how confidently we may accept or reject our hypothesis. 

Consider a set of n observed data sets in the form (xi, yi.Oi) and a theoretical prediction 
y* = f(x). Chi-squared tests involve calculating the \ 2 statistic 

2 



X 



E 



V -Hi 



(B.15) 



. =i v a > 

where y* is the predicted value of y for a given x. Clearly, if we have a perfect fit then 
X 2 = 0. In order to judge the fit correctly, one must specify the number of degrees of 
freedom k (or d.o.f.), which is defined as 



k 



n 



P, 



(B.16) 
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where p is the number of parameters (coefficients) of y* = f(x). Once the functional form 
of y* is chosen, a computer program can then find the best fit by tuning the coefficients 
such that x 2 is minimised. (Typically the program outputs the best fit coefficients, their 
errors and the resulting x 2 -) This fit is usually called the null hypothesis Hq. 

For a given number of degrees of freedom, the chi-squared distribution is defined as 



where x represents x 2 ■ The final step involves using the appropriate distribution to 
determine whether the null hypothesis can be accepted or rejected. Let us assume that 
for a certain null hypothesis H the chi-squared statistic is x 2 — u - m this case, the 
probability Q{x 2 \k) of getting x 2 > u if H is true is given by 



F k (x) 



1 



for x > 



(B.17) 



2tr(fc/2) 



i 2 e 2 , 




(B.18) 



For example, if Q{x 2 \k) is small, then it is considered unlikely that the discrepancies be- 
tween observation and prediction are chance fluctuations and hence, H may be rejected. 
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